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i We review recent progress in 2D gravity coupled to d < 1 conformal matter, based on 

Os ■ 

a representation of discrete gravity in terms of random matrices. We discuss the saddle 
point approximation for these models, including a class of related 0{n) matrix models. 
D ! For d < 1 matter, the matrix problem can be completely solved in many cases by the 



introduction of suitable orthogonal polynomials. Alternatively, in the continuum limit 
the orthogonal polynomial method can be shown to be equivalent to the construction of 



H ' representations of the canonical commutation relations in terms of differential operators. In 
the case of pure gravity or discrete Ising-like matter, the sum over topologies is reduced to 
the solution of non-linear differential equations (the Painleve equation in the pure gravity 
case) which can be shown to follow from an action principle. In the case of pure gravity 
and more generally all unitary models, the perturbation theory is not Borel summable and 
therefore alone does not define a unique solution. In the non-Borel summable case, the 
matrix model does not define the sum over topologies beyond perturbation theory. We 
also review the computation of correlation functions directly in the continuum formulation 
of matter coupled to 2D gravity, and compare with the matrix model results. Finally, we 
review the relation between matrix models and topological gravity, and as well the relation 
to intersection theory of the moduli space of punctured Riemann surfaces. 

6/93, submitted to Physics Reports 
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0. Introduction 



It was proposed some time ago [|1| tliat tlie integral over the internal geometry of a 
2D surface can be discretized as a sum over randomly triangulated surfaces. The use of 
such a lattice regularization allows the partition function of 2D quantum gravity coupled to 
certain matter systems to be expressed as the free energy of an associated hermitian matrix 
model. This matrix realization can frequently be solved by means of large techniques 
(see also @ and references therein), and the solutions restricted to fixed topology of 
the two dimensional spacetime are found to be in agreement with the continuum Liouville 
results of Q (as we shall review in detail later on here). 

More recently a continuum limit that includes the sum over topologies of two 

dimensional surfaces was defined for certain matter systems coupled to 2D quantum gravity. 
The continuum limit specific heat for these models was moreover found to satisfy an 
ordinary differential equation, in principle allowing a full non-perturbative solution. 

The above progress suggests hope for extracting nonperturbative information 

from string theory, at least in some simple contexts. A prime obstacle to our understand- 
ing of string theory has been an inability to penetrate beyond its perturbative expansion. 
Our understanding of gauge theory is enormously enhanced by having a fundamental for- 
mulation based on the principle of local gauge invariance from which the perturbative 
expansion can be derived. Symmetry breaking and nonperturbative effects such as instan- 
tons admit a clean and intuitive presentation. In string theory, our lack of a fundamental 
formulation is compounded by our ignorance of the true ground state of the theory. 

String theory is an attempt to overcome the difficulties encountered in the quantization 
of 4D gravity by replacing particles by string-like one dimensional objects, which describe 
some two dimensional "worldsheet" E as they evolve in time (interactions are encoded in 
the genus of the surface E). Polyakov showed that such theories could be interpreted as 
theories of two dimensional quantum gravity, in which the world-sheet is exchanged with 
the space-time, and the string coordinate X^{a), cr G E, is considered as a D-dimensional 
"matter" field defined on E. 

In string theory we thus wish to perform an integral over two dimensional geometries 
and a sum over two dimensional topologies, 

topologies 
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where the spacetime physics (in the case of the bosonic string) resides in the conformaUy 
invariant action 



Here ^, v run from 1, . . . , iD where D is the number of spacetime dimensions, G^J^^{X.^ is the 
spacetime metric, and the integral T>g is over worldsheet metrics. TypicaUy we "gauge-fix" 
the worldsheet metric to g^^ = Q'^bah-, where 99 is known as the Liouville field. Following 
the formulation of string theory in this form (and in particular following the appearance 
of [§), there was much work to develop the quantum Liouville theory (some of which is 
reviewed in sec. ^ here) . 

The method of using a discretization of the string worldsheet to incorporate 

in the continuum limit simultaneously the contribution of 2d surfaces with any number 
of handles, makes it possible not only to integrate over all possible deformations of a 
given genus surface (the analog of the integral over Feynman parameters for a given loop 
diagram), but also to sum over all genus (the analog of the sum over all loop diagrams). 
This progress, however, is limited in the sense that these methods only apply currently 
for non-critical strings embedded in dimensions D < \ (or equivalently critical strings 
embedded in D < 2), and the nonperturbative information even in this restricted context 
has proven incomplete. Due to familiar problems with lattice realizations of supersymmetry 
and chiral fermions, these methods have also resisted extension to the supersymmetric case. 

In addition an investigation of the large order behavior of the perturbative (topo- 
logical) expansion shows that terms at large orders k have a typical (2/c)! behavior. The 
perturbation series are divergent, and for the most interesting models (pure gravity, unitary 
models) the series are non-Borel summable because all terms are positive. Perturbation 
theory does not define unique functions, and in the pure gravity case it can be shown that 
the real solution of the Painleve equation has unphysical properties. It is thus conjectured 
that in the non-Borel summable case, the matrix model does not define the sum over 
topologies beyond perturbation theory. 

The developments we shall describe here nonetheless provide at least a half-step in 
the correct direction, if only to organize the perturbative expansion in a most concise 
way. They have also prompted much useful evolution of related continuum methods. Our 
point of view here is that string theories embedded m. D < \ dimensions provide a simple 
context for testing ideas and methods of calculation f^. Just as we would encounter 
much difficulty calculating infinite dimensional functional integrals without some prior 





4 



experience with their finite dimensional analogues, progress in string theory should be aided 
by experimentation with systems possessing a restricted number of degrees of freedom. 

Other review references on the same general subject can be found in []T^. In this 
review we concentrate mainly on the properties of D < 1 systems coupled to 2D gravity. 
Many interesting recent developments in the field, including study of issues of principle 
such as topology change in 2D quantum gravity, and as the relation to recent work on 
d — 2 black holes in string theory, are based on D = 1 matter coupled to gravityS (see e.g. 

El)- 



1. Discretized surfaces, matrix models, and the continuum limit 

1.1. Discretized surfaces 

We begin here by considering a "D = dimensional string theory", i.e. a pure the- 
ory of surfaces with no coupling to additional "matter" degrees of freedom on the string 
worldsheet. This is equivalent to the propagation of strings in a non-existent embedding 
space. For partition function we take 

Z = 5:/D,e-^^+™, (1.1) 

h 

where the sum over topologies is represented by the summation over h, the number of 
handles of the surface, and the action consists of couplings to the area A = f y^^, and to 
the Euler character x = J \f9 R = 2 — 2h. 

(Recall that the Einstein action for pure gravity with cosmological term reads S{g) = 
J d'^Xy/g{KR + A), in which is the metric tensor, R the scalar curvature and K,A 
are two coupling constants. The cosmological constant A multiplies the volume element. 
In two dimensions classical gravity is trivial because the scalar curvature term J ^/gR 
is topological (Gauss-Bonnet theorem) and thus does not contribute to the equations of 
motion. In the quantum case, however, even two dimensional gravity is non-trivial because 
large quantum fiuctuations may change the genus of the surface and the partition function 
hence involves a sum over surfaces of all genus. In addition, on higher genus surfaces there 
are non-trivial topological sectors.) 

^ which, together with the Liouville field, results in a d = 2 dimensional target space in the 
critical string interpretation 
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Fig. 1: A piece of a random triangulation of a surface. Each of the triangular 
faces is dual to a three point vertex of a quantum mechanical matrix model. 

The integral J T>g over the metric on the surface in ( p. . 1|) is difficult to calculate in 
general. The most progress in the continuum has been made via the Liouville approach 
which we briefly review in sec. |^. If we discretize the surface, on the other hand, it turns out 
that ( p. . 1| ) is much easier to calculate, even before removing the finite cutoff. We consider in 
particular a "random triangulation" of the surface ^ , in which the surface is constructed 
from triangles, as in fig. 1. The triangles are designated to be equilateral,!! so that there 
is negative (positive) curvature at vertices i where the number Ni of incident triangles is 
more (less) than six, and zero curvature when A^^ = 6. Indeed if we call V , E, and F 
the total number of vertices, edges, and faces respectively, of the triangulation, due to the 
topological relations 2E = J2i 3F = 2E (a relation obeyed by triangulations of 

surfaces, since each face has three edges each of which is shared by two faces) , all quantities 
can be expressed in terms of the A^i's. The discrete counterpart to the Ricci scalar R at 
vertex i is Ri = 2n{6 — Ni)/Ni, so that 

j ^ R -^^Ati{1 - Ni/Q) =At:{V -\F) = Ati{V - E + F) = Atix , 

i 

^ We point out that this constitutes a basic difference from the Regge calculus, in which the 
link lengths are geometric degrees of freedom. Here the geometry is encoded entirely into the 
coordination numbers of the vertices. This restriction of degrees of freedom roughly corresponds 
to fixing a coordinate gauge, hence we integrate only over the gauge-invariant moduli of the 
surfaces. 
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coincides with the simphcial definition which gives the Euler character x- The discrete 
counterpart to the infinitesimal volume element is cXi = Ni/3, so that the total area 
IS"! = just counts the total number of triangles, each designated to have unit area. 

(The factor of 1/3 in the definition of is because each triangle has three vertices and 
is counted three times.) The summation over all such random triangulations is thus the 
discrete analog to the integral / Vg over all possible geometries, 



genus h 



Vg 



(1.2) 



random 
triangulations 



In the above, triangles do not play an essential role and may be replaced by any set 
of polygons. General random polygonulations of surfaces with appropriate fine tuning of 
couplings may, as we shall see, have more general critical behavior, but can in particular 
always reproduce the pure gravity behavior of triangulations in the continuum limit. 



1.2. Matrix models 



We now demonstrate how the integral over geometry in (|1 . 1|) may be performed in 
its discretized form as a sum over random triangulations. The trick is to use a certain 
matrix integral as a generating functional for random triangulations. The essential idea 



goes back to work on the large limit of QCD, followed by work on the saddle point 
approximation 0. 

We first recall the (Feynman) diagrammatic expansion of the (0-dimensional) field 
theory integral 

;i.3) 



d<^ -(/?V2 + Av?V4! 
: e 



-oo v27r 

where (p is an ordinary real number.i In a formal perturbation series in A, we would need 
to evaluate integrals such as 



4! 



:i.4) 



Up to overall normalization we can write 



d 



2k 



e 



j=o 



;i.5) 



j=o 



The integral is understood to be defined by analytic continuation to negative A. 
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Since -§jQ'^ = J^^ applications of d/dJ in the above need to be paired so that 
any factors of J are removed before finaUy setting J = 0. Therefore if we represent each 
"vertex" Xtp'^ diagrammaticaUy as a point with four emerging hnes (see fig. 2b), then ( |1.4| ) 
simply counts the number of ways to group such objects in pairs. DiagrammaticaUy we 
represent the possible pairings by connecting lines between paired vertices. The connecting 
line is known as the propagator (see fig. 2a) and the diagrammatic rule we have 

described for connecting vertices in pairs is known in field theory as the Wick expansion. 



(a) (b) 

Fig. 2: (a) the scalar propagator, (b) the scalar four-point vertex. 

When the number of vertices n becomes large, the allowed diagrams begin to form 
a mesh reminiscent of a 2-dimensional surface. Such diagrams do not yet have enough 
structure to specify a Riemann surface. The additional structure is given by widening 
the propagators to ribbons (to give so-called ribbon graphs or "fat graphs" ) . From the 
standpoint of (|1.3|), the required extra structure is given by replacing the scalar ip by an 
N X N hermitian matrix M*j. The analog of ( |1 . 5| ) is given by adding indices and traces: 



trMV2 ■ d d -trM^/2 + trJM 



M 



e ' M'S- •■■M*",- = — — e 

Jl Jn ^ ^ - ^ ^ _ 



d 


d 






d 


d 







J=0 



tr JV2 

e 



;i.6) 



j=o 



where the source is as well now a matrix. The measure in ( |1.6|) is the invariant dM = 
rii ni<j dReM^j dlmM^j, and the normahzation is such that ^ ^ 

calculate a quantity such as 



^! Jm 



e-''^'"/^{trMY , (1.7) 



we again lay down n vertices (now of the type depicted in fig. 3b), and connect the legs 
with propagators (M'^jM^i) = Sj (fig. 3a). The presence of upper and lower matrix 
indices is represented in fig, 3 by the double linej and it is understood that the sense of 



This is the same notation employed in the large expansion of QCD [|12|. 
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the arrows is to be preserved when hnking together vertices. The resulting diagrams are 
similar to those of the scalar theory, except that each external line has an associated index 
z, and each internal closed line corresponds to a summation over an index j = 1, . . . , N. 
The "thickened" structure is now sufficient to associate a Riemann surface to each diagram, 
because the closed internal loops uniquely specify locations and orientations of faces. 



(a) (b) 
Fig. 3: (a) the hermitian matrix propagator, (b) the hermitian matrix four-point vertex. 

To make contact with the random triangulations discussed earlier, we consider the 
diagrammatic expansion of the matrix integral 



(with M an N X N hermitian matrix, and the integral again understood to be defined 
by analytic continuation in the coupling g.) The term of order g"^ in a power series 
expansion counts the number of diagrams constructed with n 3-point vertices. The dual 
to such a diagram (in which each face, edge, and vertex is associated respectively to a dual 
vertex, edge, and face) is identically a random triangulation inscribed on some orientable 
Riemann surface (fig. 1). We see that the matrix integral ( p..8|) automatically generates all 
such random triangulations.i Since each triangle has unit area, the area of the surface is 
just n. We can thus make formal identification with ( |1 . 1|) by setting g = e~^. Actually 
the matrix integral generates both connected and disconnected surfaces, so we have written 
on the left hand side of As familiar from field theory, the exponential of the 

connected diagrams generates all diagrams, so Z as defined above represents contributions 



^ Had we used real symmetric matrices rather than the hermitian matrices M, the two indices 
would be indistinguishable and there would be no arrows in the propagators and vertices of fig. 3. 
Such orientationless vertices and propagators generate an ensemble of both orientable and non- 



orientable surfaces |13]. 
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only from connected surfaces. We see that the free energy from the matrix model point of 
view is actually the partition function Z from the 2d gravity point of view. 

There is additional information contained in A^, the size of the matrix. If we change 
variables M — > M\/N in ( p..8| ), the matrix action becomes A^tr(— ^trM^ -\- gtxM^), 
with an overall factor of N% This normalization makes it easy to count the power of 
associated to any diagram. Each vertex contributes a factor of A^, each propagator 
(edge) contributes a factor of A^"^ (because the propagator is the inverse of the quadratic 
term), and each closed loop (face) contributes a factor of A^ due to the associated index 
summation. Thus each diagram has an overall factor 

^V-E+F ^ ^ N''-2h ^ (;L_9) 

where x is the Euler character of the surface associated to the diagram. We observe that 



the value N = makes contact with the coupling 7 in (IT). In conclusion, if we take 
g = e~^ and A^ = e''', we can formally identify the continuum limit of the partition function 
Z in ( p..8|) with the Z defined in The metric for the discretized formulation is not 

smooth, but one can imagine how an effective metric on larger scales could arise after 
averaging over local irregularities. In the next subsection, we shall see explicitly how this 
works. 

(Actually (|l]^) automatically calculates (|1T|) with the measure factor in (|1.2| ) cor- 
rected to J2s \G(S) \ ' ^^^^^ \ is the order of the (discrete) group of symmetries of the 
triangulation S. This is familiar from field theory where diagrams with symmetry result 
in an incomplete cancellation of l/n!'s such as in ( |1.4| ) and ( |1 . 7| ) . The symmetry group 



G{S) is the discrete analog of the isometry group of a continuum manifold.) 

The graphical expansion of ( |1.8| ) enumerates graphs as shown in fig. 1, where the 
triangular faces that constitute the random triangulation are dual to the 3-point vertices. 
Had we instead used 4-point vertices as in fig. 3b, then the dual surface would have square 
faces (a "random squarulation" of the surface), and higher point vertices {gk/N^/'^~^)ti 
in the matrix model would result in more general "random polygonulations" of surfaces. 
(The powers of A^ associated with the couplings are chosen so that the rescaling M 
M^/N results in an overall factor of A^ multiplying the action. The argument leading to 



^ Although we could as well rescale M M/g to pull out an overall factor of N/g"^, note that 
N remains distinguished from the coupling g in the model since it enters as well into the traces 
via the N x N size of the matrix. 
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( p. .91) thus remains valid, and the power of continues to measure the Euler character of 
a surface constructed from arbitrary polygons.) The different possibilities for generating 
vertices constitute additional degrees of freedom that can be realized as the coupling of 2d 
gravity to different varieties of matter in the continuum limit. 

1.3. The continuum limit 

From ( |1.9| ), it follows that we may expand Z in powers of A^, 

Z{g) = N^Zoig) + Z^{g) + N-^Z^ig) + ■ ■ ■ = N^-^'^Z^ig) , (1.10) 

where Zh gives the contribution from surfaces of genus h. In the conventional large 
limit, we take N —>■ oo and only Zq, the planar surface (genus zero) contribution, survives. 
Zq itself may be expanded in a perturbation series in the coupling g, and for large order n 
behaves as (see for a review) 

Zo{9) - Y.^'-'-'ig/Qcr - {9c - 9?-' ■ (1-11) 

n 

These series thus have the property that they diverge as g approaches some critical coupling 
g^. We can extract the continuum limit of these surfaces by tuning g ^ g^. This is because 
the expectation value of the area of a surface is given by 

(A) = (n) = ^lnZo(^7)~^— 
og 9 -9c 

(recall that the area is proportional to the number of vertices n, which appears as the 
power of the coupling in the factor g^ associated to each graph). As g ^ g^, we see that 
A ^ oo so that we may rescale the area of the individual triangles to zero, thus giving a 
continuum surface with finite area. Intuitively, by tuning the coupling to the point where 
the perturbation series diverges the integral becomes dominated by diagrams with infinite 
numbers of vertices, and this is precisely what we need to define continuum surfaces. In 
general, surfaces of large area are connected with the large order behavior of the Taylor 
series expansion in powers of g and therefore to the singularity of Zh{g) closest to the 
origin. 

There is no direct proof as yet that this procedure for defining continuum surfaces is 
"correct", i.e. that it coincides with the continuum definition ( |1 . 1| ) . We are able, however, 
to compare properties of the partition function and correlation functions calculated by 
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matrix model methods with those properties that can be calculated directly in the contin- 
uum (as in the early work of WM, and which we shall review in later sections here). This 



gives implicit confirmation that the matrix model approach is sensible and gives reason to 
believe other results derivable by matrix model techniques (e.g. for higher genus) that are 
not obtainable at all by continuum methods. 

One of the properties of these models derivable via the continuum Liouville approach 
is a "critical exponent" 7str, defined in terms of the area dependence of the partition 
function for surfaces of fixed large area A as 

Z{A) ~ a(^=*'-2)>^/2-i . (1.12) 

(Note that if we consider (|1.12|) restricted to genus zero, i.e. with x = 2, then we see that 
7str coincides with 7 of ( p. . 1 1| ) . ) To anticipate some relevant results, we recall that the 
unitary discrete series of conformal field theories is labelled by an integer m > 2 and has 
central charge D = 1 — 6/m(m + 1) (for a review, see e.g. [0), where the central charge is 
normalized such that D = 1 corresponds to a single free boson. If we couple conformal field 
theories with these fractional values of D to 2d gravity, the continuum Liouville theory 
prediction for the exponent 7str is (see subsec. ^5^ ) 

Tstr = ^{D-1- ^iD-l){D-25)) = . (1.13) 

The case m = 2, for example, corresponds to D = and hence 7str = —\ for pure gravity. 
The next case m = 3 corresponds to D = 1/2, i.e. to a 1/2-boson or fermion. This is the 
conformal field theory of the critical Ising model, and we learn from (|1.13| ) that the Ising 



model coupled to 2d gravity has 7str = ~\- Notice that ( |1.13| ) ceases to be sensible for 
D > \. This is the first indication of a "barrier" oX D = \. 

In sections we shall present the solution to the matrix model formulation of the 
problem, and the value of the exponent 7str provides a coarse means of determining which 
specific continuum model results from taking the continuum limit of a particular matrix 
model. Indeed the coincidence of 7str and other scaling exponents (see sec. |]) calculated 
from the two points of view were originally the only evidence that the continuum limit 
of matrix models was a suitable definition for the continuum problem of interest (note 
however a subtlety in the comparison for non-unitary models). The simplicity of matrix 
model results for correlation functions has spurred a rapid evolution of continuum Liouville 
technology so that as well many correlation functions can be computed in both approaches 
and are found to coincide, as we shall review in sec. |6|. 
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1.4- The double scaling limit 

Thus far we have discussed the naive N ^ oo Umit which retains only planar surfaces. 
It turns out that the successive coefficient functions Zh{g) in ( |1.1C1| ) as well diverge at the 
same critical value of the coupling g = (this should not be surprising since the divergence 
of the perturbation series is a local phenomenon and should not depend on global properties 
such as the effective genus of a diagram) . As we shall see in the next section, for the higher 



genus contributions (|1.11|) is generalized to 



n 

We see that the contributions from higher genus (x < 0) are enhanced as g g^. This 
suggests that if we take the limits N ^ oo and g ^ g^. not independently, but together 
in a correlated manner, we may compensate the large high genus suppression with a 
g ^ g^ enhancement. This would result in a coherent contribution from all genus surfaces 

To see how this works explicitly, we write the leading singular piece of the Zh{g) as 

Zh{g) - Mg - 9cY'-'^''^''/' . 

Then in terms of 

^-i^N(^g_g^)i2-y...)/2 ^ (1.15) 
the expansion ( p..lO| ) can be rewritten! 

Z = K-^h + fl + ^^f2 + --- = Yl ^'""^ fh ■ (1-16) 

h 

The desired result is thus obtained by taking the limits N oo, g ^ g^ while holding 



fixed the "renormalized" string coupling k of (|1.15|) . This is known as the "double scaling 
limit" . 



^ Strictly speaking the first two terms here have additional non-universal pieces that need to 
be subtracted off. 
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2. The one— matrix model 

In order to justify tlie claims made in the previous section, we introduce some for- 
malism to solve the matrix models. Since the integrand in (|1.8|) depends only on the 
eigenvalues of the matrix M, we can factorize the integration measure into the product of 
the Haar measure for unitary matrices and an integration measure for eigenvalues. The 
integration over unitary matrices is then trivial and we can rewrite the partition function 
in the form 

e^^/dMe-''^<^^'^/ndA.A^(A)e-2-''<^-' , (2.1) 

i=l 

where we now allow a general polynomial potential V{M) = M"^ + y^^^^ 
the Aj's are the eigenvalues of the hermitian matrix M, and 

A(A) = n(A, - A.) (2.2) 

i<j 

is the Vandermonde determinant.i (In appendix we give a more formal derivation for 
the appearance of the Vandermonde determinant based on the group metric.) Due to 
antisymmetry in interchange of any two eigenvalues, (|2.2| ) can be written A(A) = det A'^"^ 
(where the normalization is determined by comparing leading terms). In the case = 3 
for example we have 

1 A A^ 

(A3-A2)(A2-Ai)(A3-Ai) = det I 1 A2 x\ 

1 A3 Ai 



(|2.1| ) may be derived via the usual Fadeev-Popov method: Let Uo be the unitary matrix such 
that M = UqA'Uo, where A' is a diagonal matrix with eigenvalues A^. The right hand side of ( p. ID 
follows by substituting the definition 1 = J H.dXidU S{UMU^ - A) A^(A) (where J dU = 1). 
We first perform the integration over M, and then U decouples due to the cyclic invariance of the 
trace so the integration over U is trivial, leaving only the integral over the eigenvalues Ai of A. 
To determine A(A), we note that only the infinitesimal neighborhood [/ = (! + T)Uo contributes 
to the U integration, so that 

„ w 

1= l[dX,dU5''\uMU^ -A)A\X)= / dT 5^(^-1) ([T, A']) A2(A') . 

i=l 



Now [T,A']ij = Tij{X'j — X'i), so ( p. 2] ) follows (up to a sign) since the integration dT above is over 
real and imaginary parts of the off-diagonal Tij's. 
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2.1. The large N limit: steepest descent 

The large limit of the matrix models considered here was originally solved by saddle 
point methods in H. For this procedure, it is convenient to change the normalization in 



the integrand of (|2.1| ) and consider instead 



where the coupling constant g plays the role of the cosmological constant. To describe 
pure gravity, we recall that only one ctfc is needed, for instance one can use only triangles. 
More general models correspond to additional degrees of freedom on the surface. 

In the conventional large limit, in which according to ( p..lO| ) only surfaces with 
the topology of the sphere contribute, the integral (|2.3|) can be evaluated by steepest 
descent. The Vandermonde determinant leads to a repulsive force between eigenvalues 
which otherwise would accumulate at the minimum of the potential V . The saddle point 
equations that come from varying a single eigenvalue in (|2.3|) are 



This equation can be solved by the following method: We introduce the trace of the 
resolvent of the matrix M 

^U) = —ti—^ = —S" . (2.5) 

i 

Multiplying eq. ( p.4| ) by l/(Ai — z) and summing over i, we find 

^ 1 , ^ - J- x: ^^MM . ,,0, 

This equation is analogous to the Riccati form of the Schrodinger equation, the wave 
function -0 being related to u by Nu{z) + NV {z) / {2g) = i^' /i^. The eigenvalues are 
the zeros of the wave function. In the large limit, we can neglect N~'^uj'{z) (this is the 
WKB approximation, but note that the equation can more generally be used to study the 
convergence of the distribution of zeros toward its limit). 

In this limit the distribution of eigenvalues p(A) = Ylii '^(A — A^) becomes continuous, 

and 
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Note that the normahzation condition J p{X')dX' = 1 is the analogue of the Bohr- 
Sommerfeld quantization condition. Eq. (|2.4| ) can now be rewritten 



or equivalently 



u(z + iO)+u(z-iO) = --V'(z) . (2.i 

9 



Finally, eq. ( p.6|) becomes 

iv^z) + -V'{z) uj{z) + -^R{z) = , (2.9) 
9 4:9 

where 



R{z) =Agj A\p{\)^^^^^^^^^ (2.10) 

is a polynomial of degree / — 2 when V is of degree /. Note that the coefficient of highest 
degree of R is fixed by the normalization of p(A) while the remaining coefficients depend 
explicitly on the eigenvalue distribution. 

The eigenvalue density p(A) is extracted from u){z) via the relation 

p(A) = ^ {u{z + zO) - u{z - zO)) , (2.11) 

and from (|2.3|) , we can write the gravity partition function as 

Z = (^j dXdpp{X)p{p)ln\X- j dXp{X)V{X)^ . (2.12) 



The solution. The solution to eq. (|2.9|) i^ 



u,{z) = ^{-V'{z) + a{z)) , (2.13) 



2^7 



^ Recall that this solution could also be determined indirectly by first solving the homogeneous 
equation, i.e. with the r.h.s. of (|2.8| ) set to zero, by looking for a function that has a cut on the 
support of p and takes opposite values above and below the cut. Such a function has square root 
branch points and is therefore the square root of a polynomial. The inhomogeneous equation 
has particular solution —V'{z)/2g, and additional constraints arise from the condition that io{z) 
behaves like —1/z for \z\ large. This equivalently determines the solution given with the degree 
I — 2 polynomial R{z) having fixed highest degree coefficient. 
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where (j{z), up to the normahzation, is the singular part c<;sing(2;) of uj{z) 



a{z)=2guj,,^^{z) = ^{V'{z)) - R{z) . (2.14) 

Generically uj{z) has 2{l — 1) branch points corresponding to the roots of the polynomial 
— R. Therefore the support of p(A) is formed of / — 1 disconnected pieces. In the 
simplest case, when the potential has only one minimum, we expect a single connected 
support and thus only two branch points. It follows that the polynomial V'^ — R must 
have / — 2 double roots and this yields / — 2 conditions that fully determine R. 

For later purposes, it is convenient to give a more explicit representation of the 
one-cut solution. The simplest one-cut solution of the homogeneous equation (|2.8[ ) is 
a/(z — ai)(z — 02). Dividing uj{z) by this function, we can transform ( |2.8[ ) into a discon- 
tinuity equation, 

uj{z + iO) uj{z-iO) _ 1 V'{z) 

i\/{a2 - z){ai - z) -i^ [a-i - z){z - ai) ig^ [a-i - z){z - ai) 

It follows that 



27r£? Ja, A - z ^(a2 - A)(A - ai) 

The discontinuity equation defines a solution up to regular terms, and the large 1^1 
behavior implies their absence. The large \z\ behavior yields also two other conditions 
which determine ai and 02: if we expand the r.h.s. of ( p.l5|) for \z\ large we find first a 
constant piece which must vanish, so that 

"2 V'(\) 

dA ^ ^ = , (2.16a) 

ai v(a2 - A)(A - ai) 

and then we find a term proportional to 1/z whose residue is known, so that 

"2 \V'(\) 

dA ^ _ 27rg . (2.166) 



ai 



^/{a2 - A)(A - ai) 



These equations can be written in another convenient form by replacing the integrals 
over the cut with contour integrals, 

v/(A-ai)(A-a2) 

V'(\) ^ ^ ^ 

dAA ^ ^ = AiTvg . 

^/(A - ai)(A - 02) 
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After the change of variable 



A = 2 + - ai + a2 H — 

2 16z 



eqs. (|2.16a, ^ ) take the form 



f^V'{X{z))=g, (2.18a) 
— V'{Xiz))=0. (2.186) 



r/ie partition function in the large N limit. The partition function can be calculated 
from eq. ( |2.12| ) but it is more convenient to work with the derivative 

— = -(try(M)>~— ydAp(A)y(A), 

differentiated again after multiplying by g^: 

This leads us to consider the function 

dg dg J X — z dg 



whose real part vanishes according to eq. ( |2.8| ). Moreover since for z large, u behaves as 
— 1/2, Q{z) also behaves as —1/z. Finally since uj{z) behaves near ai,a2 as y/z — ai, its 
derivative behaves at most as 1/ z — a^. The unique solution is 

n^z) = -{iz-ai){z-a2)y^^^. (2.21) 

Transforming the integral in (|2.19|) into a contour integral, 



dg \ dg J 2nv J ^(A - ai)(A - aa) 

and differentiating a last time with respect to g gives 

d{gu{g)) _ iV2 fdai / y(A) 



4in\dgJ (A-ai)3/2(A-a2)i/2 



9^7 7 (A-ai)i/2(A-a2)3/2 
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Integrating by parts in eqs. ( p.l7|) , we can generate the two needed integrals so that finaUy 
we have 

= 2^2/ln|ai-a2| (2.23) 
dg dg 

which we shall use to determine the singular part of the partition function. 

Even potentials. For an even potential the cut end-points take equal and opposite 
values ±a. Eq. (|2.186|) is automatically satisfied and ( p.l8a|) becomes 

g = j^V'iz + a^/Az) . (2.24) 

For the quartic potential 

V{\) = \\^ + \\\ (2.25) 

for example, we find0 

uj{z) = ^(-z-z^ + {z^ + l + ia^) ^z'' - aA , (2.26) 
Ig \ / 

with 



2 

From (|2.11j ), we find that 



a' = I (-1 + ^1 + 12^7) 



p(A) = (A^ + 1 + W) Va2 - A2 . (2.27) 
ZTxg 



It is important to notice here that the partition function (|2.3| ) to leading order in large 
(the spherical limit) has an analytic continuation0 from g > Q to g < Q and the first 
singularity arises aX g = g^ = —1/12, at which point Z has a square root branch point. The 
existence of the (7 < region can be understood as follows: the number of planar diagrams 
increases only geometrically, while the barrier penetration effects for (7 < responsible 
for the divergence of perturbation theory behave as e~^^^^^ . The latter are therefore 
exponentially suppressed for large (for details see subsec. |7. /[ ) . 

Note that, as should have been expected, the singularity in g occurs at a point where 
a zero of the singular part a{z) of uj{z) coalesces with an end-point of the cut. At g^, a{z) 
becomes 



aiz) = {z'-alf' 



a" 



With respect to the conventions in eqs. (17a, 6) of we have g g/A, A^/gr, 
a? /Ag, and oj = —F. 

Due to our choice of normahzation the average matrix is anti-hermitian when g is negative. 
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with ttc = a{g^). 

The continuum limit. We want now to study the singular behavior of functions near 
(the continuum Umit). Let us blow up the neighborhood of the cut end-point. Due to 
the symmetry of the potential both end-points play a role, for generic potentials only one 
end-point would be relevant. We set x = 1 — q/Qc and z = a^l — s). Then for s and x 
small with s = 0{^/x), cr{z) has a scaling form: 



a{z) {4/3f/^ {s + ^^/x) ^^y/x- s and - Wsmg(s) = 2^6 a(2) . 

The singular part of the partition function can be calculated directly using equations ( |2.22| ) 
and ( p.23| ). Near at leading order, we can replace the explicit powers of g by g^, and 
integrate to find 

g'^Z" ~ 2N'^ In |ai — 02] + less singular terms . 



For the potential ( |2.25| ), this gives 

The partition function thus takes the form 

Z = -^iV^x^/^ less singular terms , (2.28) 

which implies 7str = —1/2. (Note that the Legendre transform from the fixed area partition 
function Z{A) of ([1121) gives / dA A(T-*'-2)>^/2-ie-^^ ~ ^{2-j,,,)x/2^ ^^^^^ ^ero 

Z{x) ~ a;2-T=''-.) 

2.2. Multicritical points 

We have seen that a critical point is the result of a confluence of a regular zero of the 
singular part a of uj{z) with a cut end-point. By taking potentials of higher degree, which 
thus depend on more parameters, we can adjust these parameters in such a way that m — 1 
zeros of a reach a cut end-point for the same value g = g^. At the critical point, (^{z) will 
then have the form0 

a{z) ^z'^-^/^z-by/^ . (2.29) 



^2 in a minimal realization because it could have additional irrelevant zeros 
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We have assumed here that at the critical point the end-point of the cut is at 2; = 0, 

and that the potential is generic (not even) such that the other end-point of the cut is 
z = b>0. 

Note that the form ( |2.2y| ) determines both the critical potential and since from 
( p.l3|) we have 

a{z) = V'{z) -2gz-^ + {z-"^) . (2.30) 

It follows that 

V'iz) = [z^il-b/zy/')^ , (2.31) 

where the subscript -|- means the sum of the terms with non-negative powers of z in the 
large z expansion. Equivalently, we can write 

= sdw) + <^ - I' - ■ 

where B is the usual ratio of F functions 

r(a)r(/3) 



B{a,l3) 



Note that the minimal potential goes as 2"^+^ for z large, and therefore the corresponding 
matrix integral can be defined only by analytic continuation for m even. From ( p.29| - |273T| ) , 
we also have 

Um+l 

g^= Bim + 1/2,3/2) . (2.32) 

For X = 1 — g / g^ small (but not critical), the zeros at z = and the end-point of the 
cut, now at 2; = a, will split. The other end-point (at z = b) also moves but it is easy to 
verify that this effect is negligible at leading order in x. For a and z small, with z = 0{a), 
a again assumes a scaling form. For z large, this scaling form must match the small z 
behavior of the critical form ( |2.29| ). Thus we find 

a,, = b'/'{a-zy/'a^-'P„,.^{z/a) , 

where Pm{z) is a polynomial such that ~ z"^ for z large. Moreover from ( |2.2U| ) and 
( ^ ), we infer 

^ = -26-V2(a_^)-i/2 . (2.33) 
dg 
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This imphes that aU terms with positive powers in the expansion of a^c for z large 
must vanish, except for the first term which is independent of a. It foUows that 



(where again the subscript means the sum of terms with non-negative powers in the large 
z expansion), and the term of order z~^l'^ in Usc is proportional to . Comparing with 
( p.33|) , we find that da^/dg is a constant. Eq. ( |2.33| ) then determines a{z) up to a 
multiplicative constant which is fixed by the large z behavior, and we obtain the integral 
representation 

aUz) = r ds {s - z)-'/h--' . (2.34) 



S(m,l/2) 

Expanding this expression in powers of z and comparing with the leading term of expression 
( ^) , we find 

Q m 

^ = -2mb-^B(m,l/2) , 
dg 

and thus the relation between a and x is 

2(m + l)(^-j =x{a) = x. (2.35) 



Combining this relation with eqs. (|2.13| , |2.32| , |2.34| ), we obtain a useful representation 
for the singular part o^sing of u in the scaling limit: 



" , dx 
ds 



Ceasing = 6-1/2 I ds — {s- z)-^/^ . (2.36) 



Prom ( p.22| ) and (|2.23 ), we then immediately obtain the singular part of the partition 
function: 

^s'ing(^) ~ -2iV2(a/6) , (2.37) 

and thus 

We conclude in particular that the susceptibility exponent 7str takes the value 

7str = -1/m 

(see comment following eq. ( |2.28| ) for the identification of 7str)- As anticipated at the 



end of subsec. 1.^, we see that more general polynomial matrix interactions provide the 
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necessary degrees of freedom to result in matter coupled to 2d gravity in the continuum 
limit. 

To compare the result ( ^.381 ) for m = 2 with ( |2.28| ), it is necessary to relate the 
normalizations of a; in both calculations. Note that J dzu){z) is normalization independent, 
because the related quantity J dAp(A) = 1, so that comparison between the scaling forms 
of the singular part of uj{z) in both calculations determines the relative normalizations of 
X and z. We find here for Zsing one half of the result (|2.28| ). This surprising result has a 
simple explanation [l^] : For even potentials both cut end-points contribute to the partition 
function and this yields an additional factor of 2. In what follows, we shall largely restrict 
for reasons of simplicity to even potentials, so we need to keep this peculiarity in mind. 
The relative normalization of Z^ing and a will also be useful when discussing large order 
behavior of the topological expansion. Finally we indicate in sec. ^ how the 0{n) gas loop 
on a random surface []T^ can be investigated by a generalization of the method presented 
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The loop average. Eventually we shall also consider correlation functions of quantities 
of the form L(s) = -^tr M*, which create in the surface a loop of length s. Here we examine 
the behavior of the loop average (L(s)), for large loop length s, in the case of a general 
critical point. Such an average represents a sum over surfaces with a loop of length s as 
boundary. 

In the large N limit, we find 



[L{s)) = I dXp{X)y 



Using the analogue of eqs. ( |2.iy| - pT^^ , we obtain 



dg 2m J ^(A - ai)(A - 03) 

For s large, the integral is dominated by the neighborhood of max{|ai|, |a2|}. Let us 
assume \a2\ > \ai\. Thus we have 

^{g{Lis))) ~ , ^ , 
dg ^ ' s^ca _ a2/ai 

We shall now assume that 02 corresponds to the cut end-point where all zeros of the 
resolvent coalesce in the critical limit (a condition satisfied by the minimum potentials 
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( p.31|) when A is shifted such that V'{X) ~ A when A — > 0), and we consider the case of an 
j^th Qj-(iQy- critical point. It foUows from eq. ( |2.35[ ) that 

a2 - [a2]c oc x^/"" , 

and therefore 



ln{|L(s)) ~ s lna2c + const ■ s 



X 



l/r 



The first term on the r.h.s. is a short distance effect and can be canceUed by a matrix 
renormahzation so that we have 

[ln{L{s))]^^^^sx'/^ . (2.39) 

This behavior shows that the distance should be rescaled by a factor x^/'" in the continuum 
limit. For m = 2 (pure gravity), this agrees with the area scaling as 1/x as we argued in 



subsec. For m > 3, however, the result is different, leading to a difficulty in identifying 
X with the cosmological constant (which by definition is coupled to the area). Let us define 
the length scale by ( |2.39| ) and call n ~ x^/"* the cosmological constant. If we then define 



the exponent 7 by the behavior (|1.12| ) of the fixed area partition function Z{A), we find 



(see comment following eq. ( ^.28[ )) 

and thus 7 = 3/2 — m. We see therefore that there are several ways to define a string 
exponent, depending on the reference parameter. A similar problem will arise in general 
for non-unitary models (of which the m > 3 multicritical one matrix models comprise a 
particular subclass), as will be discussed in sec. ^ 

Note finally that the normalized average ;^(tr e''^), which yields a weighted super- 
position of loops of different length, has the same behavior in the continuum as L{s). It is 
often used instead of L{s) because it has simpler algebraic properties. It is in particular 
related to the trace of the resolvent u!{z) by 

(Wis)) ^ ^(tr e^^O = • (2.40) 

2.3. The method of orthogonal polynomials 

The steepest descent method allows a general discussion of the large N limit. It is 
difficult however to calculate the sub leading orders in the expansion and therefore to 
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discuss perturbation theory to all orders. We now present another method that allows us 
to recover previous results and to extend them to all orders in 

This alternative method for solving (|2.1D makes use of an infinite set of polynomials 



P^(A), orthogonal with respect to the measure 

/oo 
dA e-^(^) P„(A) P^(A) = 5r,m ■ (2.41) 
-oo 

The Pn's are known as orthogonal polynomials and are functions of a single real variable A. 
Their normalization is given by having leading term Pn(A) = A" + . . ., hence the constant 
Sn on the r.h.s. of ( p.41| ). Due to the relation 



A(A) = det Ar^ = det P,_i(A,) (2.42) 



satisfied by the Vandermonde determinant ( p.2|) (recall that arbitrary polynomials may 
be built up by adding linear combinations of preceding columns, a procedure that leaves 
the determinant unchanged), the polynomials can be employed to solve (p7l]). We 
substitute the determinant det Pj_i(Ai) = 'Y1,{~^Y Y\k Pik-^i^k) for each of the A(A)'s in 
( p.l|) (where the sum is over permutations ik and (— l)'^ is the parity of the permutation). 
The integrals over individual A^'s factorize, and due to orthogonality the only contributions 
are from terms with all Pi(Aj)'s paired. There are A^! such terms so ( |2.1| ) reduces to 

Z 



£ 7i-,7r' k j 

N-1 N-1 



N-1 N-1 

eN-k 
'k 

i=0 k=l 



where we have defined fk = Sk/sk-i- The solution of the original matrix integral is thus 
reduced to the problem of determining the normalizations Sfc, or equivalently the ratios fk- 
In the naive large limit (the planar limit), the rescaled index k/N becomes a 
continuous variable ^ that runs from to 1, and fk/N becomes a continuous function 
/(^). In this limit, the partition function (up to an irrelevant additive constant) reduces 
to a simple one-dimensional integral: 

i^^ = ^E(l-^/^)l^^~ /"'de(l-01n/(0 ■ (2.44) 

k •'^ 

To derive the functional form for /(O, we assume for simplicity that the potential 



V{X) in (|2.41| ) is even. Since the P^'s from a complete set of basis vectors in the space 
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of polynomials, it is clear that AP^(A) must be expressible as a linear combination of 
lower Pi's, XPn{X) = Yl7=o '^i -^i-W (with = s^^ J e~^XPnPi). In fact, the orthogonal 
polynomials satisfy the simple recursion relation, 

XPn = Pn+l+rnPn-l , (2.45) 

with r„ a scalar coefficient independent of A. This is because any term proportional to P„ 
in the above vanishes due to the assumption that the potential is even, f e~^A Pn Pn = 0. 
Terms proportional to Pi for i < n — 1 also vanish since / e-^PnXPi = (recall AP^ IS a 
polynomial of order at most z + 1 so is orthogonal to Pn for i + 1 < n). 

By considering the quantity P^XPn-i with A paired alternately with the preceding or 
succeeding polynomial, we derive 

6 Pn X Pn—1 fn ■ 

This shows that the ratio fn = Sn/sn-i for this simple caseEl is identically the coefficient 
defined by ( p.45|) , fn = Vn- Similarly, since P!^ = nPn-i + 0(A"~^), 



nSn-l = j P'nPn-l = - j Pn^ie-"" Pn-l) = J e-"" V Pn Pn-1 ■ (2.46) 

This is the key relation that will allow us to determine r^. 

2.4- The genus zero partition function revisited 

Our intent now is to find an expression for fn = Tn and substitute into ( |2.44| ) to 
calculate a partition function. For definiteness, we take as example the potential 

V[X) = -(x^ + X^ + hX^) , 
^ ^ 2<7 V ; ' (2.47) 

with derivative gV'{X) =N {X + 2X^ + 36A^) . 

The right hand side of (|2.46| ) involves terms of the form J e~^ A^^~^ Pn Pn-i- According to 
( p.45|) , these may be visualized as "walks" of 2p — 1 steps (p — 1 steps up and p steps down) 
starting at n and ending at n — 1, where each step down from m to m — 1 receives a factor 
of rm and each step up receives a factor of unity. The total number of such walks is given 



In other models, e.g. multimatrix models, fn = Sn/sn-i has a more complicated dependence 
on recursion coefficients. 
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by (^^~^), and each results in a final factor of Sn-i (from the integral / e~^ Pn-i Pn-i) 
which cancels the Sn-i on the left hand side of ( p.46|) . For the potential ( |2.47| , p.46| ) this 
gives 

— = r„ + 2rn{rn+i + + r^-i) + 36(10 rrr terms) . (2.48) 

(The 10 rrr terms start with r^ir'^ + r'^+i + r'^^i + . . .) and may be found e.g. in [p!9| .) 

As mentioned before (|2.44| ), in the large limit the index n becomes a continuous 
variable ^, and we have r^ — r(^) and r„±i — r(^ ± e), where £ = 1/A^ (as in ( |2.44| ) we 



assume that for n,N large becomes a smooth function of n/N). To leading order in 
1/A^, reduces to 

g^ = r + 6r^ + 306r^ = W{r) 

(2 49) 

= g^+lW"\r=rMO-rcT + -- - . 

In the second line, we have expanded W{r) for r near a critical point r^ at which VF'|r=rc = 
(which always exists without any fine tuning of the parameter b), and = W{rc). We 
see from (|2.49|) that 

^-^c~ {9c -90^^^ ■ 

To make contact with the 2d gravity ideas of the preceding section, let us suppose 
more generally that the leading singular behavior of /(^) (= r(^)) for large is 

/(0//c-l~i^(^7.-^70"^ (2.50) 

for g near some g^ (and ^ near 1). (We shall see that 7 in the above coincides with the 
critical exponent 7str defined in (|1.12|) .) The behavior of ( |2.44| ) for g near g^ is then 



- 0i9c - 90-^+'\' + 2n ^^9 - 90'^+' (2.51) 

f(l-7) lo £?^(l-7)(2-7) 

~ 2n — -M - 9)-^^' - E --f^^^''~'i9/9cr ■ 

fc(l-7)(2-7) Vfcr(7) 



Comparison with (|1.12| ) shows that the large area (large n) behavior identifies the exponent 
7 in ( |2.50| ) with the critical exponent 7str defined earlier. We also note that the second 
derivative of Z with respect to x = g^ — g has leading singular behavior 

Z"ig) ~ Kg^-^g^ - ^7)-^- ~ g--" = 1)//, - l) . (2.52) 

27 



From (|2.50| ) and (|2.51 ) we see that the behavior in ( |2.49|) imphes a critical exponent 



7str = —1/2. From (|1.13|) , we see that this corresponds to the case D = 0, i.e. to pure 



gravity. It is natural that pure gravity should be present for a generic potential. With 



fine tuning of the parameter b in ( |2.47| ), we can achieve a higher order critical point, with 
W'\r=rc = W"\r=rc = 0, and heuce the r.h.s. of ( |2.49| ) would instead begin with an [r — rc)^ 
term. By the same argument starting from ( |2.5C1| ), this would result in a critical exponent 

7str = -1/3. 



General potential. For a general potential V{M) = '^pVpM'^P we find VF(r 

(2p-l)! 

p {{p-iy-r 



2 X]p (P-i)n2 'Vp^^: or equivalently in more compact form 



/dz 
V'{z + r/z) . (2.53) 
2z7r 



(In the "stairmaster" interpretation following ( p.47|) , we see that z corresponds to stepping 



up and to stepping down. The integral is non-vanishing only when there is an overall 
factor z~^, and correctly takes into account a factor of r for each step down.) Note that 
when r is identified with a^/4 (where a is the boundary of the eigenvalue distribution), 
eqs. (|2.49|) and (|2.24|) become identical (up to the change ^ g)- Eq. ( |2.53| ) also can be 



inverted to yield 

ds 



V{M) 



/ —W{sil-s)M^). (2.54) 
Jo ^ 



With a general potential V{M) in ( |2.1| ), we have enough parameters to achieve an 
^th order critical point ||2^ at which the first m — 1 derivatives of W{r) vanish at r = r^. 



The minimal potential has degree 2m and corresponds to a function W{r) of the form 

^W = %(C-(rc-r)'") , (2.55) 



m 
c 



and r{g) behaves as 



From ( |2.52| ) we then find 



rc-r{g) fgc-9 



(2.56) 



at2 2 / \ 2+1/r 

N'^m^ I g^ — g 



(2m + l)(m + l) \ g^ ' 
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identical, up to normalizations, to (|2.38|) . The form ( p.55|) for W{r) is thus associated with 



the critical exponent 7str = —l/m. Expression ( p.57|) suggests the existence of a scaling 
region in which becomes large with fixed 

Z={1- ^^jAr2m/(2m+l) _ ^2.58) 



The minimal critical model. For the minimal critical model, i.e. the model in which 
the polynomial V has the smallest degree possible at a fixed value of m, we find from the 
explicit form of W in ( 2.55 ) and the relation ( p.54|) that the critical potentials Vm{M) are 
given by0 



" Jo 



2\m-l 



(4r, 



^2.59) 



This also follows directly from the considerations of subsec. \2. 3^ , since in the even case 

a{z) = -62)m-l/2_ 

( p.59|) also shows that the term of highest degree of V{M)/g^ in the integral ( pT^ ) is 
positive for m odd, and negative for m even. In the latter case, the minimal potential 
integrals can thus be defined only by analytic continuation, as in the case of generic (not 
even) potentials. This is the source of important differences between the two parities. 



2. 5. The all genus partition function 



We now search for a solution to ( |2.48| ) that describes the contribution of all genus 
surfaces to the partition function ( p.44| ). 



Pure gravity. Let us for simplicity discuss first pure gravity. We shall retain higher 
order terms in in (|2.48|) so that e.g. (|2.49|) instead reads 



9^ = W{r) + 2r(0 (l + 156r(0) (r(^ + e) + r(^ - e) - 2r(0) + ■ ■ ■ 

^g^ + \W"\r=r^ (r(0 - rc)' + 2re(l + X^hr^) (r(^ + e) + r{i - e) - 2r(0) + 



(2.60) 



14 



This expression was established by another method in [21|. The + subscript is defined as in 



subsec. meaning the sum of terms of non-negative powers in the large M expansion. 
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where the dots mean terms of order (r — Tc)^, (r(^ + e) + ~ ^) ~ 2r(^)) (r — Tc) or terms 
which approach £^r""(^) in the smaU e hmit. We shaU eventually verify that the omitted 
terms are neghgible. 

As suggested at the end of subsec. we shaU simultaneously let ^ oo and g 



in a particular way. Since g — g^ has dimension [length]^, it is convenient to introduce a 
parameter a with dimension length and let 

9c~9 = gcO-'^z , (2.61) 

with a — > 0. Our ansatz for a coherent large N limit will be to take 

e = l/N = ^ 

so that the quantity = z^l'^ = {g ^ gc)^^^^ remains finite as g g^ and N —>■ oo. 

Moreover since the integral ( |2.44| ) is dominated by ^ near 1 in this limit, it is convenient 
to change variables from ^ to a;, defined by 

gc-9^ = 9ca^x. (2.62) 

Our scaling ansatz in this region is 

r(^) = r^{l-au{x)) . (2.63) 

If we substitute these definitions into (|2.49|) , the leading terms are of order and result 



in the relation ~ x. 

To include the higher derivative terms, we note that 

r(e + e)+ r{C - e) ~ 2r(0 = + O fe^ ^'"^ 



—rf,a^:r-^au{x) + O [a^) ~^ —rcO^u" 



where we have used 

e{d/di) = -{g/gja^/\d/dx) ~ -a^/\d/dx) , (2.64) 

which follows from the above change of variables from ^ to a: and g — gc = O{o?'). Substi- 
tuting into ( |2.6C1| ), the vanishing of the coefficient of a? implies the differential equation 

Kx = u^ - \u" with R-^ = {6rl + 906r^) /g^ . (2.65) 

30 



In the ordinary large limit, we have calculated the partition function in eq. (|2.44| ). It 
is however easy to verify that in the double scaling limit the corrections to (|2.44|) are of 
order (^N(g^ — g)) ^ oc a^, and thus remain negligible. After the changes of variables (|2.62 
2.61j ) and with the definition ( p.63|) , we then find 



Z'(z) ~ — / dx {z — x)u{x) =^ Z"{z) = —u{z) 



(2.66) 



generalizing (|2.52| ). In the double scaling limit, the second derivative of the partition 
function (the "specific heat") has leading singular behavior given by u{z) for z = {1 — 



The solution to (|2.65| ) characterizes the behavior of the partition function of pure 
gravity to all orders in the genus expansion. (Notice that the leading term is w ~ 2;^/^, so 
after two integrations the leading term in Z is z^/^ = ^"^j consistent with (|1.16| ).) Finally 
we note that the change of normalization 

z t— > pz, u{z) 1-^ p^u{pz) , (2.67) 

does not affect ( |2.66| ) but allows us to rewrite eq. ( p.65| ) as 

z = u^ - \u" . (2.68) 

The property that the r.h.s. of the equation is invariant is directly related to the property 
that r — Tc = 0{a) and e{d/d^) = 0(a^/^) (eqs. ( p.63| , p.64|) ), which was essential in 
selecting terms of the same order. This property will generalize to higher order critical 
points. 



The Painleve I equation. Eq. ( f2.68D is known in the mathematical literature as the 



Painleve I equation. One characteristic property is that its only moveable singularities (in 
the complex plane) are double poles |^^. In the normalization ( |2.68D they have residue 2, 
and, as eq. ( f2.66| ) shows, correspond to double zeros of exp Z. Since it is a second order 
differential equation, its solutions are determined by two boundary conditions. We are 
interested only in solutions that have an asymptotic expansion for z large (the topological 
expansion) that begins with the leading spherical result u{z) ~ ^/z. The perturbative 
solution in powers of = is then determined and takes the form 

u{z)^z'/\l-J2^kZ-"'/^) , (2.69) 
fc=i 
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where the Uk are all positive-lll This verifies for this model the claims made in eqs. ( p..l4| - 
1.16| ) of subsec. \1.4[ For large k, the Uk go asymptotically as (2/c)! (for details see sub- 



sec. |7.^ ), so the solution for u{z) is not Borel summable and thus does not define a unique 
function. Our arguments in sec. |l| show only that the matrix model results should agree 
with 2d gravity order by order in perturbation theory. How to ensure that we are studying 
nonperturbative gravity as opposed to nonperturbative matrix models is still an open ques- 
tion. Some of the constraints that the solution to ( |2.68| ) should satisfy are reviewed in ||2^ . 
In particular it is known that real solutions to ( |2.68| ) cannot satisfy the Schwinger-Dyson 
(loop) equations for the theory. 

Higher order critical points. In the case of the next higher multicritical point, with b 
in (|]4|) adjusted so that W' = W" = at r = r^, we have W{r) ^ + \ W"'\r=rAr - 
rcY + ■ ■ ■ and critical exponent 7str = —1/3. In general, we take g — Qc = ^^/'''^^"^"^^a^, 
and e = 1/N = a^"'^''"' so that the combination {g — g^y~'^''*-^^'^N = k^^ is fixed in the 
limit a 0. The value ^ = 1 now corresponds to z = ^^^/(Tstr-s)^ ^^le string coupling 
k"^ = 2;'''^*''"^. The general scaling ansatz is r(^) = rc(l — a~'^"^'""'u{x)), and the change of 
variables from ^ to x gives e{d/d^) ~ — a~'''^"^(9/9x). 

For the case 7str = —1/3, this means in particular that r(^) = rc(l — a^/'^w(x)), 
= and e{d/d^) ~ ~^^^^'§^- Substituting into the large limit of ( |2.48| ) gives 

(again after suitable rescaling of u and z) 

z^u"" -uu" -\{u'f + au"", (2.70) 

with a = j^. Note again that since |r—rc|^/^ = 0(a~"''="') = 0{e{d/d^)), the r.h.s. of ( |2.70| ) 
is invariant under the transformation ( p.67|) . The solution to ( |2.7C1| ) takes the form u = 
z^/^[l + Uk z~'^^^'^). It turns out that the coefficients uu in the perturbative expansion 
of the solution to ( 2.70|) are positive definite only for a < j^, so the 3**^ order multicritical 



point does not describe a unitary theory of matter coupled to gravity. Although from 
( p..l3| ) we see that the critical exponent 7str = — 1/3 coincides with that predicted for the 
(unitary) Ising model coupled to gravity, it turns out [|T9| , |2^ that ( |2.7C1| ) with a = 



The first term, i.e. the contribution from the sphere, is dominated by a regular part which 
has opposite sign. This is removed by taking an additional derivative of u, giving a series all of 
whose terms have the same sign — negative in the conventions of ( p.68| ). The other solution, with 

1 /2 

leading term — z , has an expansion with alternating sign which is presumably Borel summable, 
but not physically relevant. 
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instead describes the conformal field theory of the Yang-Lee edge singularity (a critical 
point obtained by coupling the Ising model to a particular value of imaginary magnetic 
field) coupled to gravity. The specific heat of the conventional critical Ising model coupled 
to gravity turns out (see subsec. \4-'% to be as well determined by the differential equation 



( p.70|) , but instead with a 



_2_ 
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For the general m^^ order critical point of the potential VF(r), we have seen that 
the associated model of matter coupled to gravity has critical exponent 7str = —l/m. 
With scaling ansatz r{^)/rc = 1 — a^/™'w(x), we find leading behavior u{z) ~ ^^Z™- (and 
Z ~ ^s+i/m = ^-2 expected). The differential equation that results from substituting 
the double scaling behaviors given before ( |2.7C1| ) into the generalized version of (|2.48|) turns 



out to be the m^^ member of the KdV hierarchy of differential equations (of which Painleve 
I results for m = 2). In the next section, we shall provide some marginal insight into why 
this structure emerges. 

In the nomenclature of so-called "minimal conformal field theories" (those with 
a finite number of primary fields) are specified by a pair of relatively prime integers (p, q) 
and have central charge D = Cp^q = 1 — 6{p — q)"^ /pq. (The unitary discrete series is the 
subset specified by (p, q) = (m-f l,m).) After coupling to gravity, these have critical 
exponent 7str = — 2/(p + — 1), as calculated in the matrix model. In general, the m}^ 
order multicritical point of the one-matrix model turns out to describe the (2m — 1,2) 
model (in general non-unitary) coupled to gravity, so its critical exponent 7str = — l/m 
happens to coincide with that of the w}^ member of the unitary discrete series coupled 
to gravity. The remaining (p, q) models coupled to gravity can be realized in terms of 
multi-matrix models (to be defined in sec. We shall see that the interpretation of 7str is 
slightly subtler in the non-unitary case since the identity is no longer the lowest dimension 
operator. 

2. 6. Recursion formulae more generally 



In subsec. |^. j[ , we introduced orthogonal polynomials and used them to calculate the 
matrix integral. We now introduce some additional formalism related to this polynomials 
which will prove useful in what follows. As an immediate application, it will allow us to 
discuss the problem of general (not necessarily even) potentials. 

Since the polynomials Pn{^) of ( p.41|) form a basis, we can write (in matrix notation) 

P' = AP , XP = BP . (2.71) 
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In explicit component form, these relations are equivalently written 

n-l 

P^=Y, ^r^mPm , (2.72a) 
m=0 
n+1 

XPn = BnmPm ■ (2.726) 

m=0 

The normalization of Pn (eq. (|2.41|) ) gives 

Note that the matrices A,B of ( |2.71| ) necessarily form a representation of the canonical 
commutation relations 

[B,A] = 1. 

They can be related to the polynomial V{X) by the relations 

J dn{X) XPmP-n = BnmSm = SnBmn (2.73) 

and 

= / dA ^ e-^^(^)/«)) = Ar^ms^ + s^Amn - ^Km^m , (2.74) 

where we employ the notation = [V'{B)]nm and d/j{X) = dXe~^^^^^^^ . Introducing 
the matrix S with matrix elements Snm = SmSnm, we can rewrite ( |2.73| ) and (|2.74| ) in 
matrix form, 

BS = SB^, 

(2.75) 

AS + SA^ = {N/g)V'{B)S . 



Eq. (|2.73|) , together with the defining relation (|2.726|) , shows that Bmn is different 



from only if |m — n| < 1. If we no longer assume that the function V{X) is even, the 
recursion relation ( p.72£|) between orthogonal polynomials becomes 

XPn = Pn+l + fnPn + r^^n-l . (2.76) 

where r^, are a short- hand notation for Bn,n-i, Bn,n respectively. 
Specializing ( |2.73| ) to m = n — 1, we recover 

dlJ,XPnPn-l = Sn = r^Sn-l , 
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and consequently 



Sn-l 



(2.77) 



Substituting into (|2.43|) , we find 

Z = In {s^N\) + ^ (iV - n) Inr, 

n=l 

SpeciaUzing also ( |2.74| ) to m = n — 1, we recover ( |2.46 ): 



7V-1 



^2.78) 



n 



Q— =V' 1 . 

N rijU— 1 



(2.79) 



As we have already explained, since in the case of even potentials can be expressed 
entirely in terms of r„ , ( p.79|) leads to a recursion relation for the coefficient r„ . For general 
potentials, to determine both r„ and r„, we need instead another equation, in addition to 
( p.79|) , obtained by specializing ( |2.74| ) to m = n. 



V = . 



(2.80) 



Example. Consider the potential 



V{X) = + 



(equivalent after translation to + A'^/3). The two equations ( |2.79| ) and (|2.80|) in this 

case are 

gn/N = rn{rn + rn-i) , 

1 = fl + rn + Tn+l , 



with critical values 



1/3 , r, = ±v/^ , = 2fl 



The sign of Tc is irrelevant (see however the remark at the end of subsec. [g. l\j . Setting as 
in ( p.63|) , r{g)/rc — ^ — N~'^/^u{N'^^^{1 — g/g^)), we obtain the equation 



3 2 1 // 

— = z . 



(2.81) 



We see that eqs. ( |2.68p and ( |2.81| ) are identical up to a rescaling of the function u and 
variable z. Moreover the normalization of the partition function depends only on the ratio 
of the coefficient of u" to the coefficient of u"^. We note here that this ratio is 1/6 instead 
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of 1/3 in ( |2.68| ). The partition function for a potential of the form —A + from eq. ( |2.66| ) 
is half of the previous partition function, and the double poles of the solution of eq. ( |2.81J ) 
have now residue equal to 1. In a direct saddle point technique, the interpretation of this 
result is that the partition function receives two identical contributions when the potential 
is even. 

General potentials in the spherical limit. To leading order in large A^, eqs. ( |2.79| ) and 
reduce to 

/dz 
V'(z + r/z + f) , (2.82a) 
2z7r 

/dz 
- — V'(z + r/z + r) . (2.826) 

We recognize eqs. ( |2.18a, 6|) , with r and f related to the end-points ai,a2 of the support 
of the matrix eigenvalue distribution by 

r = (ai + a2)/2 , r = (ai - a2)^/16 . 



The criticality condition for eqs. ( |2.82a, Ij^ reads 

dW dX dW dX 



Note that 



dr df dr dr 
dW dX 



(2.83) 



dr df 



and 

dx r dz ,,,, , , , r^^w 

ir = ^ (z + r/z + r) = (z^-) -—- 

or / Imz^ z r or 



so that the criticality condition (|2.83| ) can be rewritten 



dxV fdX^^ 



drr / \df 



C 



For g close to g^, we parametrize r and r as 

r = Tc ( 1 — w) , f = fc — \/rcU ■ 

Eq. (12.8261) then imphes 

dX ^ dX 

= rcU- h VtcU—z- , 

ore ore 

and thus 

u = ±u . (2.84) 

This relation is indeed satisfied in the example considered above (leading to ( p.81|) ). Rela- 
tion ( p.84| ) will be useful in next section for the general analysis of the one-matrix problem. 
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2. 7. Correlation functions 

We can also consider correlation functions of the form 



(trFi(M)trF2(M)---) 
For the simplest examples, we find 



N-l 



(trFi(M)>=5^[Fi(i?)]_, (2.85a) 



n=0 

N-l N-l 



(trFi(M)trF2(M)> = H [^i(^)]nn [^^(S 



pp 

n=0 p=0 



+ E [F,{B)] (2.856) 



I np L " ^ 'J pn 

0<n<N-l<p 



More generally, connected correlation functions may be interpreted as averages of one-body 
operators in a free A^-fermion state (for a review, see [|TT| ). 



In the large N limit, the sums are replaced by integrals. For example 

{trF^{M))r^N [ d^{^\Fi{B)\^), i = n/N. 
Jo 

In the scaling limit, the singularities come from the neighborhood of ^ = 1. In terms of 
xa^ = 1 — g/ Qc and ya^ = 1 — g^/gc {^/N = a^"*"^/™"), the above trace involves an integral 
of the form 

(trFi(M))~a-i/"^ / dy {y\Fi{B)\y) . (2.86) 

J X 

The upper bound of the integral goes to oo in the large (double) scaling limit. 

In the large limit, away from the scaling region, the average of a product of traces 
is always dominated by the product of their averages because, as can be seen from the 
example ( p.85fc|) , the additional terms are subleading by a factor due to the "locality" 

of B and the restrictions on the summation indices. 

2.8. Loop equations, Virasoro constraints 

The loop equations [^,^^ are obtained by performing the change of variables 
M M + eM'^ in the matrix integral (|2T^). We find the identity 

fc-i 

^(trM^trM'^-^-^) = (A^^) (trM'=l/'(M)) , (2.87) 

1=0 
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where (■ ■ ■) means average with respect to the matrix integral. 

It is convenient to rewrite these equations in terms of the operator G{z) = {M — 
Multiplying eq. ( p.87| ) by z~'^~'^ and summing on k gives 

N 



tvG{z)V'{M)) = -((tr G{z)) ) . (2.88) 

This expression can be further simplified by rewriting V'{M) = V'{z) + (y'{M) — V'{z)) 
and noting that tr (y'{M) — V {z))G{z) is a polynomial in z. Setting 

R{z) = ^{tr{r{M)~r{z))G{z)), 

we obtain an equivalent form of ( |2.88|) 

jV'{z){trG{z)) + ^R{z) = -{{tvG{z))') . (2.89) 

Note that the loop equations, once expressed in terms of the matrix B, can in principle be 
derived directly from eqs. ( p.75|) , but the derivation is not so straightforward. Introducing 
the loop average L{s) (eq. ( |2.4C1| )), we can also write eq. (|2.89|) as 



/ dt{W{t)W{s-t)) = —V'{d/ds){W{s)) 
Jo 9 



The l.h.s. can be interpreted as the operation of gluing two boundaries together. 

Virasoro constraints. Let us set g = 1 for convenience here and parametrize the 
potential V{M) as 

Calling the matrix integral 'd{t) = , we can rewrite the set of loop equations ( |2.87| ) in 
the form of a set of Virasoro constraints (see e.g. 



Ck^{t) = , k>-l 

with 



^ ^ dtmdtn ^ ^„ 



^ dtmdtn ^„ tm+fc 
m+n=k m=U 

With respect to the trivial variable to, we have the additional equation 

dto 
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Large N limit. In the large limit, ( 2.89 ) can be solved by noting from expression 
( p.856| ) that the average of the product of traces becomes asymptotically equal to product 
of averages. Setting 

1 



u,{z) = -{tvG{z)) 



we find 



1 



1 



iv'{z) + -V'{z) u{z) + -^R{z) = . 



9 ■ ■ V" 
We recognize the form (|2.9|) of the saddle point equations. 



3. A general method: the canonical commutation relations 

It is possible to study the recursion formulae ( p.79| ) for general polynomials V{X). In 
the continuum limit one finds a non-linear differential equation of more general type for 
a scaling function u. A simpler algebraic method has been found, however, which easily 
generalizes to the several matrix problem. 

It is convenient to introduce normalized orthogonal polynomials Un, 

satisfying 

y d;u(A) n^jll,, = ■ (3.1) 
We now redefine matrices A and B in terms of the 11^, 

n— 1 n+1 
Iln ~ ^ ^ Anm^rm -^Iln = ^ ^ -^nmllm • (3-2) 

With this new definition, the matrix B is symmetric. In the even potential case, in terms 
of the coefficients r„ introduced in eq. ( p.76| ) the recursion formula for the orthogonal 
polynomials becomes 

= y/rn+l Un+l + ^/r^Un-i ■ (3.3) 
Instead of as in ( p.75| ), the equation for A now reads 

A + A^ = —V'{B), (3.4) 
9 
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while the commutator relation remains [-B, A] = 1. It is convenient to shift A and introduce 
the matrix 

C^A-^^V'iB) = l{A-A^) (3.5) 

representing the operator d/dA acting on the orthogonal functions e"^^*-^-*/^^ H^. Then 
C is antisymmetric and satisfies the same commutation relation as A, 

[B,C] = 1. (3.6) 

A basic property. A remarkable property can be proven (for details see sec. ^): Let B 
be a Jacobi matrix B (a symmetric matrix with Bmn = for \m — n\ > 1) that satisfies 
the commutation relation ( p.6|) , with C an antisymmetric local matrix (i.e. with Cmn = 
for \m — n\ > I). Then there exists a lower triangular matrix A {A^ri = for n > m) and 
a polynomial V{X) of degree / + 1 such that 

jV'{B) = A + A^, C = ^{A-A^) ^ [A,B] = 1 . 

The diagonal matrix elements of the last equation yield an equation equivalent to the 
difference between (|2.79 ) taken for two consecutive values of n. 



Since the coefficients in V are implicitly determined by the criticality conditions, the 
original problem can thus be entirely reformulated in terms of the matrix B satisfying the 
commutation relation ( |3.6| ), from which the singular part of the partition function Z can 
be calculated. We now proceed to take the large and scaling limits directly in these 
expressions. 

3.1. The large N limit 

We now show that in the double scaling limit B, C become differential operators and 
that the commutation relation ( |3.6|) determines the string equation. 

In the limit of large (and thus n), ^ = n/N can be treated as a continuous variable 
and then the matrix B can be expanded for r„ near Tc, 

= 2n^ H + e + O (r^ - n^, e 



with e = = a^~^=*^ As we have discussed in subsec. \2.d^ for a general critical point 
we take 9 = Qc^cl"^ so that the combination (1 — g/g^)N'^/^'^~~'''^''^ = 2; is fixed in the 
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limit a — s> 0. The change of variables from ^ to a; = a~^(l — gi/g^ gives e{d/d^) ~ 
— a~'^=*'^(9/9x). 

The general scaling ansatz is then 

r(^) = rc(l - a"^'^=''M(a;)) . 

The leading term in the expansion of B, proportional to the identity, does not contribute 
to the commutation relation ( p.6| ). The leading corrections are of order a"^'''^"^. Only two 
terms contribute which together form a second order differential operator Q: 

[{r-'/^B - 2) n]^ ~ Qn„ , (3.7) 

with 

Q = d^ - u{x) (3.8) 
(where d is a notation for d/dx). Note that the rescaling ( p.67|) , 

X 1-^ px, u{x) ^— s> p^u{px) , 

transforms Q into p^Q. The formal hermiticity of the operator Q follows directly from the 
symmetry of B. 

General potentials. The recursion relations between orthogonal polynomials, general- 



izing (|3.3| ), are 



In the large N scaling limit, the matrix B now becomes the differential operator 

(sn). = (2V?; + + I {^-^^ + n„ + ^n;:} + ■ ■ ■ . 

Recalling now the relation u = ±u of eq. (|2.84| ), we see that one sign yields a trivial 
Q operator and the other replaces u by 2u in Q. The sign which yields a trivial Q has 
the following interpretation: In this case it is not the polynomials Un which are smooth 
functions of n but instead (— l)"!!^^. Taking into account this property, one finds that in 
all cases 2u will now be the solution of a differential equation described below, and the 
only overall effect is to multiply Zsing by a factor of 1/2. 
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The matrix C . In the double scaling limit, the matrix C also becomes proportional to 
a differential operator P, which moreover is formally anti-hermitian because C is antisym- 
metric. The degree of the operator P is seen to be at most 21 — 1 when the multicriticality 
conditions are met. Let us prove this property for the minimal critical potential. At the 
critical point we expect C oc d^'~^ = P, terms of lower degree in d arising from devia- 
tions from criticality. At the critical point the dependence on n of -Bn,n and Bn,n+i can 
be neglected which means B^^n+k ~ bk- Therefore powers of the matrix B become very 
simple after Fourier transformation. Let us introduce the variable z — e*^ and consider 
the Fourier transform of B 

J2 Bn,n+kZ^ ~ hkZ^ = biz) . 
k k 

Then we have 

J2F{B)n,n+kz''-^F{b{z)) . 
k 

In the same limit, the Fourier transform of the matrix C has to be proportional to 
021-1 £qj. g-^oii ^}iug (^2; — For the minimal potential, the matrix C is then 

determined up to a multiplicative constant: 

/ 

J2 Cn,n+kZ^ ~ C{Z) OC Z-\z - lf'-\z + 1) . 
k=-l 

The function V'{h{z)) differs from c{z) only by the terms with negative powers 

[c{z)]^ = \[V\b{z))]^ , 
i.e. V is determined by the expansion of c{z) for z large. Let us then set 

t^\{z + l/z) ^ b{z) 2Bn,n+lt + Bn,n • 

Conversely we choose z = t + — 1 such that z large corresponds to t large. In the 
variable t, V'{b) is a polynomial which can be determined by expanding c{z) for t large. 
Noting \/2z = y/t-1 + Vt + 1, we find 

c{z) oc (V^ - (v^+ oc 2'(t + l)V2(t _ 1)^-1/2 . 

It follows that 

v'{bit))o.it+iy/\t-if-'/\ 
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which, up a shift of variables, is identical to ( p.31|) . 

In the double scaling limit, it will be convenient to normalize P by 



P = d^^-^ + O (d^^-^) . 

The operator P, obtained by collecting all terms of the same order in a, must transform 
multiplicatively under (|2.67| ). This implies P ^ p^^~^P. Moreover writing explicitly the 



powers of a in the relation between C and P, we find 
The first factor a'^+T-tr = jsf 

comes from the relation (^]^) between C and P, and the 
second factor a"*-^'"^'''^^"^ from the change of variables ed/d^ t-^ d/dx. To check the 
consistency of the scaling ansatz, note that the powers a cancel in the commutator [P, C] 
as they should, — 27str — (2 — 7str) — {21 — l)7str — which indeed yields 7str — — !//• 
The commutator [P, Q] has the form 

[P.Q] = K , 

where iiT is a constant which can be set to 1 after a rescaling of the form ( |2.67| ), i.e. of the 
variable proportional to the deviation from the critical coupling constant. 

The problem of formulating the double scaling limit of matter systems coupled to 2d 
gravity is thus reduced to finding solutions of the canonical commutation relation 

[P,Q] = 1, (3.9) 

a problem we shall discuss in the next section. 



Finally starting from eq. (|2.44|) and introducing the quantities x, z, u, we verify that 



the singular part of the partition function is always given by eq. ( 2.66 ) 

d^Z 



^smg 



dx^ 



—u{x) 



the result for the pure gravity case. Note that this equation is left unchanged by a rescaling 
of the form (|2.67|) , a property which will be used systematically later on here. 



Remark. In this formulation, the argument of the orthogonal polynomials, i.e. the 
variable A, becomes the eigenvalue of a Schrodinger-like differential operator. The spectrum 
of Q yields the asymptotic distribution of the eigenvalues of the hermitian matrices. 
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Application: the resolvent in the large N limit. We can use eq. (|2.86| ), applied to 
[M — z]~^, to compare the expression for the singular part of the resolvent obtained from 
a steepest descent calculation (eq. ( |2.34|) ) with the resolvent of the operator Q at leading 



order. Considering eq. (|3.7| ), it is natural to set a'^^^z — 2 — rc ^^'^z. We then find 



with 



uj{z) = -itr [M - z]-^ ~ -r;^/^a''/^iJo{z) 



u{z) = / dy{y\Q + z\y) . 

J X 



At leading order, the function u{x) can be replaced by the leading term u{x) ~ x^/^ and 
the resolvent can be calculated in the semi classical limit, i.e. the non-commutation between 
d and x can be neglected. It follows that 



dp 1 1 y'{s)ds 



2% p'^ + u{y) - z 2 - 5 ' 

with u{y) = s, and c = w(a~^). We recognize that the singular part of u), for c ^ oo, is 
identical, up to normalizations, to expression ( |2.3ti| ) with z \—>- z. 

The scaling limit beyond the spherical approximation. We now consider the scaling 
limit of (m^) and find 



y dt y dt' {t\{Q - z)'^\t') {t'\{Q - zy^\t) = regular function of z . (3.10) 



It follows quite simply that the spectrum of Q cannot be discrete . Indeed at an isolated 
eigenvalue of Q, corresponding to a normalizable eigenvector (pri{x), the resolvent (Q — 
z)~^ has a pole located aX z = —Qn- The l.h.s. of ( p.lO| ) thus has a double pole with residue 



4(t)) + I dt I dt' ^lit) ^lit') = I dt^l( 



where we have used the normalizability of fn- This residue does not vanish except in the 
large x limit, where it decreases as exp(— 2 dt ^^/u(t)). The r.h.s. of ( |3.1C1|) , on the other 
hand, is a regular function of z: a discrete spectrum for Q thus results in a contradiction. 
Note, however, that since the residue vanishes faster than any power for x large, this effect 
is invisible in perturbation theory. 
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3.2. Construction of the differential equations 

The differential equations following from (|3.6| ) , [P, Q] = 1 , may be determined directly 
in the continuum by relating this problem to the corresponding KdV flows (see appendix 
A). The operator P that can satisfy this commutator is constructed as a fractional power 
of the operator Q = — u of ( |3.8| ) (or more generally we can take Q to be a q^^ order 
differential operator in which case with P a p^^ order operator, we construct the (p, q) 
minimal models mentioned at the end of subsec. \2.3i coupled to gravity) . The differential 



equations describing the {21 — 1, 2) minimal model are given by P = ^/2.|i6| 

[q'-'/\ Q] = 1 , (3.11) 

where Q^,. indicates the part of Q'"^/^ with only non- negative powers of d. In terms of 
the Gelfand-Dikii polynomials Ri, we find 

[Q'-'/',Q]=4R[ . (3.12) 

After integration and rescaling of x, the equation [Q^j^ = 1 thus takes the simple 

form 

{I + ^)Ri[u] = X . (3.13) 
The quantities Ri can be for instance calculated from the recursion relation 

R'i+i = lR!{' -uR'i-^u'Ri . (3.14) 

An action principle. The Ri^s satisfy as well a functional relation that allows us to 
write eq. (|3.13| ) as the variation of an action. Indeed it is shown in app. A. 2 that 



du 



j dxRi+M = -{l + \)Ri[u] . (3.15) 
The differential equation (|3.13|) therefore results as the variational derivative with respect 



to u of the action 



S = j dx {Ri+i + xu) . (3.16) 

(We treat the above integral formally here and ignore throughout that physically relevant 
boundary conditions on u typically preclude existence of such integrals.) 
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In subsec. -^.j, we give a generalized discussion of this "string equation". 
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3.3. Relevant perturbations and interpolation between multi critical models 

Instead of fixing to critical values all but one parameter in the potential (i.e. the 
one corresponding to the most singular perturbation), it is possible to let them approach 
criticality at a rate related to in the large limit. The partition function then becomes a 
function of a set of new relevant scaling parameters and this defines a general "massive" 
model interpolating between multicritical points. The operator P then takes the form 



I 

k-l/2 

+ 

k=l 



k=i 

and the corresponding string equation is 

I 

= -x + J2 {k + l)t(^k) Rk[u] 
k=i 

I 



(3.17) 



fc=0 



with t(o) = — 4x. Note that the constant i(i-i) can be eliminated by an irrelevant trans- 
lation of u{x) (adding a constant to u means adding a regular contribution oc to 
the partition function). 

Equivalently, using (|3.15|) , eq. (|3.17|) can be seen to follow from the action 



S = 2 I dx( Vt 



I 

Rk+i[u\^ 



'^^^ (3.18) 



k=0 



with t(o) oc X, the trace of a pseudo-differential operator being defined in app. A. 3. We 
shall see that the form of the action ( |3.18| ) generalizes to [p, q) models. 

The scaling parameters can be thought of as dimensionful couplings to certain 
operators 0(^k) added to the pure gravity action, 

^ = ^grav + 5^t(fc) [ d^CVgO^k) (3.19) 
fc>0 



(where the pure gravity action ^Sgrav is as in (lO), suitably defined by matrix model tech- 



niques as in sections or by topological gravity or Liouville methods as in sections |],^ 
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One of these operators is the identity operator coupled to the (renormahzed) cosmologi- 
cal constant. Other operators are so-caUed gravitational descendants, arising from the 
scaling limit of traces of certain polynomials of the matrix in the one hermitian matrix for- 
mulation. The most general one-matrix model solution corresponds to some perturbation 
of the pure gravity model by these operators: eq. (|3.17|) describes the space of all these 



perturbations. (See for the transformation from the basis of the t(fe)'s to the scaling 
operators of the theory.) 

KdV flows. As discussed in app. the dependence of the specific heat u on the 
parameters t(^k) is given by the higher KdV flows 

^ -Q= -J-u = ^Rk+i[u] = Q] . (3.20) 



Using the commutativity of the higher KdV flows, it is straightforward to verify consistency 
of (pl^ with ( CT ). 



As usual in fleld theory, differentiation of the partition function with respect to a 
parameter of the potential generates a correlation function with the insertion of the 
corresponding operator In particular, the string susceptibility —u = log Z is a two- 
point correlator of the operator C(o) • —u= (C'(o)C'(o))- We also have (9/9t(fc)) ((9(o)C^(o)) = 
{0(0)0 (0)0 (k)) = {d/dx){0(o)C^(k))j where 0(^k) is as above the appropriate scaling oper- 
ator that couples to ( p.20|) therefore identifles Rk+i[u] = {0(o)0(k)) as the 2-point 
function of the operator O(o) with 0(k) , and we can rewrite the string equation ( |3.17D and 
the action ( |3.18| ) in terms of these 2-point functions. 

Scaling dimensions. Previous relations can be used to calculate correlation functions 
on the sphere. The string equation ( |3.17| ) becomes an algebraic equation for the string 
susceptibility u. The deviation x = to from the critical coupling constant provides a scale 
to the theory and we can assign it the dimension 1. Then u = ((9(o)0(o)) oc x~^^^'' has 
dimension — 7str, i-e. l/l for the l-th critical point. Each coupling constant carries a 
dimension 1—l/k. Correlation functions of the operator O^^) are obtained by differentiating 
u with respect to and integrating twice with respect to x. It follows for example that 
the one- and two-point functions exhibit the scaling behavior 

{0(k)) oc x^-^-^+^^Z' {0(k)0(k')) oc 0;-^-+'=/'+^'/^ .... 
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This terminology is borrowed from topological field theory, see sec. ^ for details. 
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This behavior is related to the KPZ scahng of conformal operators "dressed" by gravity 
(see subsec. 1^ 



The identity or puncture operator. We see from the previous scahng relations that 
the operator (9(o) has the smallest dimension and is thus the most singular. As we have 
defined the coupling constant g, i.e. multiplying the whole potential V{M), it couples to all 
operators. It is therefore not surprising that the most singular survives in the double scaling 
limit. From the analysis of sec. ^we would be tempted also to identify (9(o) as the puncture 
or identity operator coupled to the cosmological constant, as in the pure gravity case. As 
we have already seen in subsec. however, this identification is inconsistent with the 



behavior of macroscopic loop averages. The behavior in (|2.39|) suggests instead that t;_2 is 
coupled to the identity, C^(z-2) = -^5 because t(z-2) = At ~ x^/'. In terms of /x, the partition 
function scales as Z ~ /j'-^^^'^ (correspondingly 7str(A*) = 3/2 — /), and (Oj-^)) ~ 
so we assign to 0(fc) the dimension ^(k) = {I + k + l)/2 — {I + 1/2) + 1 = {k + 2 — l)/2 in 
terms of /x. This conclusion is supported by a direct calculation in the continuum Liouville 
formulation (see sec. 

In what follows, when we refer to coupling constant, without other qualification, we 
shall mean the parameter coupled to the operator of lowest dimension, here 0(o): which 
may or may not be the renormalized cosmological constant. 



4. Multi-matrix models 

4.1. Solving the multimatrix models 

By a method of orthogonal polynomials, it is also possible to solve models involving 
integration over several matrices M^°'\ The basic identity is 

(2) 



f -trVifM^i)) +ctrM(i)M( 



oc 



A (A(2)) 



IS 



where the M*^") are N x N hermitian matrices, the A^" their eigenvalues, and A (A) 
the Vandermonde determinant det(A'^~^) (the proof of this identity is much more involved 
than in the one-matrix case, see appendix ^ for a sketch of the proof). When the action 
5" has the form 

a=l ct=l 
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it follows that 



, -tr>5(M(«)) 



a=l 



n dA;")A(AW)e ^A(A('^-i)), 

a=l,!^— 1 
1=1, AT 

generalizing ( |2.3D . 

The result ( [4.3| ) depends on having Cq,'s that couple matrices along a line (with no 
closed loops so that the integrations over the relative angular variables in the M*^")'s can be 
performed). Via a diagrammatic expansion, the matrix integrals in (|4.3| ) can be interpreted 
to generate a sum over discretized surfaces, where the different matrices M*-"-* represent 
u — 1 different matter states that can exist at the vertices. The quantity Z in (|4.3| ) thereby 
admits an interpretation as the partition function of 2D gravity coupled to matter. Note 
that matter has only a finite number of states (by taking u ^ oo, one can represent a 
D = 1 model, i.e. a single free boson, coupled to gravity). Furthermore the only possible 
symmetry which can be implemented corresponds to reflecting the line about its center. 
Therefore these matrix models can only implement an Ising-like Z2 symmetry on a random 
lattice. 

Generalized orthogonal polynomials. To solve the matrix models we define generalized 
orthogonal polynomials Iln{X), satisfying 

I d^(AW,...,A(^-i))n^(AW)n4A(^-i)) =5^, , (4.4) 

(generalizing (|3.1j)) where the measure dfx is defined by 

d/x(AW,...,A(^-i))=e J+2.aC«A -Qd;,(a)_ ^45) 

a=l 

To derive recursion formulae, we insert A^"^ and d/dA*-"-* respectively in the integral ( [4.4|) . 

Let us denote by dfXa-i the measure ( |4.5|) in which the integration is restricted to the 
first a — 1 variables. It is then convenient to introduce a matrix associated with A*^"-* 
and defined by 

A(") j d^^-iUm{X^^^) = [Bo.]^^j d;Ua_in,(A(i)) . (4.6) 
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We can also define B^, using a similar expression but in which the roles of A*^^) and A*^'^ 
are exchanged, i.e. by integrating over A("+i), . . . , A(^-i). Then, multiplying by UniX^''~^^) 
and n^(A'^^)) respectively, and integrating over all A's, we obtain 

Finally we define matrices Ai and A^-i, 

ntn = [^i]_n„ , K = [^-i]_ff. , (4.7) 

and an additional matrix = A^_^. 

Inserting d/dA*-"-* in ( |4.4| ) gives 

Ai + ciS2 = y/(Si) , 

C„_lSa_i +Cc,Sa+i = 1/^ (Sa) for 1< CK < - 1 , (4.8) 

with 



(4.9) 



[Si,Ai] = l, 
= 1 . 

Eqs. ( liTSD imply that 

Ca S„] = Ca-i = 1 , (4.10) 



and thus that [A^-i, B,j-i\ = 1, consistent with (|4.9| ). 

Inspired by the one-matrix case (|3.5| ), we can also introduce matrices Cq, defined by 



Ci = - iy/ (Si) , 

C« = Ca_iS«_i - (S«) for a > 1 . 

It follows from these definitions and eqs. ( [4.10| , f4.9| ) that 

[-Bq, Cq] = 1 , 

C« = (B^) - c^B^+i for a < - 1, (4.11) 



Let us call 1^ the degree of the polynomial V^. It is then easy to verify from ( ^^q ) 
that [Ba]mn is non- vanishing only for 

1^ — 1 a — l 

n - 1) < ^ - ^ < n (^/^ - 1) • (4.12) 

/3=a+l /3=1 
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This uniform bound on the "range" of the operators Sq,, henceforth on the operators Cq, 
through (|4.11| ) , is the heart of the result obtained by Douglas . In the large "double- 
scaling" limit (with all but one of the couplings in ( [4 .21 ) tuned to critical values), all these 
operators become differential operators. The uniform bound on the range above induces a 
bound on the order of the differential operators. 

Z2 symmetry. When Vq, = and = Ci,_i_q, the matrix problem has a Z2 

symmetry corresponding to the mapping of matrices M^") ^^ M^^""). We can then 
choose n„ = n„. It follows that A^^i = Ai, B^^i = Bi. In addition this symmetry yields 
the relations 



B„ 



and hence Cj = — Ci^-q, 



The free energy. We now assume for simplicity ^2 symmetry. We normalize, as in the 

the relations 



one-matrix case, n„(A) = X^/ y/s^ -f- . . . . We then find from eqs. (^?6 

[Bl]n,n+1 = = {Sn+i/ Sn)^^'^ , 

As in (|2.78|) , the matrix model free energy is given by 

F = In Z = In (A^! sqSi ...sn-i), 

N-l 

= In (A^! s^) + ^ (A^ - n) Inr^ . 



(4.13) 



n=l 



The two-matrix model. Since it has been shown that it is sufficient to consider the 



two-matrix model to generate the most general critical point [Q, let us just write the 
previous equations in this special case 



Ai + cB2 = VliBi) , A2 + cBi = V^{B2), c[Bi, B2] = 1 . 



In the ^2 symmetric case the equations become particularly simple 



(4.14) 



A + cB^ = V'{B), [A,B] = 1. 



Note that in the various relations (|4.14|) the matrices Ai,A2 can be completely eliminated 
by writing them 



c[B: 



2 I mn 



[v;{B 



1 )\mm 



[V^{B^)]mn forn> 



m 



(4.15) 
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and c[Bi,B2] = 1. 

Let us now briefly indicate, generalizing the arguments of subsec. 1\ , how in the large 
limit the potentials Vi and V2 can be chosen (multicriticality conditions) such that the 
matrices Bi , B2 become linear combinations of two differential operators P, Q of degree 
p, q respectively. We recall that at criticality the matrix elements of the various matrices 
depend only on the differences between their two indices, in such a way that a Fourier 
transformation renders the previous relations particularly simple. Let us introduce two 
functions 61(2), 62 (^) 



^l(^) ~ X] [Bi\n+k,nZ'' , 62(2) ~ [B2]n,n+k 

fc=-l k=-l 



Both functions behave as for z small. Then the relations (^4.15| ) become 



c[b2{z)U = [v;mi/z))U , c[h{z)U = [V^Ml/z))U , 

where as elsewhere the subscript + means the sum of terms of non-negative powers. 

If we impose that in the large limit at criticality Bi and B2 become proportional 
respectively to d^, d*^ (d means for example differentiation with respect to the index), then 
for 2; — > 1 we have bi{z) oc {z — 1)^, b2{z) oc {z — 1)'^. The simplest example is 

bi{z) = {z-ir/z, b2{z) = {z-iy/z . 

We introduce two new variables ti, ^2 via 

h=bi{l/z) t2=b2{l/z) . 

Inverting these relations, we choose the roots such that z large corresponds to ti, t2 large. 
Then the two critical potentials Vi, V2 are determined by 

Vi{t^)=cb2{z{h))^, Vi{t2)=cbi{z{t2))^ . 

4-2. The continuum limit 

In this subsection, it is convenient to return to the normalization of previous sec- 
tions, i.e. multiply the whole potential by a factor N/g as in (|2.3|). To study the double 
scaling limit, we introduce the renormalized deviation from the critical coupling constant 
X = a~^{l — q/qc) where e = 1/N = a'^~'^=*^ The arguments which follow are then a 
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simple generalization of those of subsec. 3.1. The matrices Bi and Ci generate differential 



operators Q, P of finite order say q > p, which still satisfy the "string equation": 

[P,Q]=1, (4.16) 



(as in ( |3.9| ), the commutator has been normalized for convenience to 1). The operators P 
and Q are obtained by collecting the terms of order a"^'''^"^ and a~^"'"^="^a~'^"''="^ respectively. 
Since the commutator is independent of a, we obtain the value of the string susceptibility 
exponent: 

2 

-P7str + (-2 + 7str) - ?7str = ^ 7str(P, <?) = -7 ; TT • (4-17) 

The operator Q can always be written 

Q = d'? + Vq-2{x)d'i-^ + ■■■ + vo{x) . (4.18) 

(By a change of basis of the form Q —>■ f~^{x)Qf{x), the coefficient of d*^"^ may always 
be set to zero.) The continuum scaling limit of the multi- matrix models is thus abstracted 
to the mathematics problem of finding solutions of ( [4.16| ) . 

The various coefficient functions involved in eq. ([4.18| ) are scaling functions of the 
coupling constant x. In units of x, the "grade" of d = is —1, so that the gradelil of 
Vq-a is —a for an operator Q of overall grade —q. 

The function Vq-2 can be identified (up to normalization) in the continuum scaling 
limit with the second derivative of the free energy with respect to the coupling constant 
(here proportional to x). Equivalently we can write Vq-2 oc d^logZ = (CminC'min) in 
terms of the 2-point function of the operator coupled to x, which, according to the analysis 
of subsec. \3.3^ , is the operator Omin of smallest dimension (i.e. the most singular in the 
continuum limit). 
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This notion of grade is related to the conventional scahng dimensions of operators, see for 



instance the end of subsec. 3.S. It can also be used to determine the terms that may appear in 



many equations, since these wih only relate terms of overall equal grade. 
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4-3. String equation solution and generalized KdV flows 

The differential equations ( [4.16| ) may be constructed as follows. Let L = Q^/"^ denote 
the 1*^* order pseudo-differential operator (i.e. a formal series with inverse powers of d) 
whose q^^ power is equal to Q. Then, by a theorem of the most general form P can 
have is 

p 

P = Y,^,{U)+ , (4.19) 

i=i 

where the subscript +, as before, indicates that we truncate the pseudo-differential op- 
erators to their differential part (non-negative powers of d). We will soon see that the 
(p, q) minimal model {p and q relatively prime) is described by the critical equation with 

[(gp/^)+, g] = 1 . (4.20) 



Let us concentrate on this case first. The equation ( |4.20| ) is a coupled differential 
system for the q — 1 coefficients fg_2,...,fo of Q. Using the fact that Q = L'^ and 
= {L^)+ + {L^)- commute, we write 

[P,Q] = -[(LP)_,L^] , (4.21) 

and expand the commutator in formal powers of d. If (L^)- = ci{x)d~^ + ■ ■ ■, then the 
leading term qcid'^~^ -(-■■■ in ( [4.21|) must vanish, giving rise to an integration constant 
ci = i^q-i. We can now proceed by writing (L^)- — i'q-iL~^ = C2d~^ + ■ ■ ■ . By the same 
argument, we get a second integration constant: C2 = and so on, until we reach 

the d*^ term in the commutator ( [4.21j ). Taking (L^)- = Yli=i ^q-iL~^ + c^-id^""^ + ■ • • , 
qc'g-i = 1 integrates to Cq-i = x/q, where we absorb the last integration constant in a 
shift of X. 

The occurrence of these integration constants is crucial for the following reason: the 
dependence of L (therefore of all the coefficients of Q) on these constants Ui {i = 1, . . . , q—1) 
can be shown to be given by the first q — 1 generalized KdV fiows. Namely, in terms of 
the "times" ti = qfi/i (z = 1, 2, . . . , g — 1), L satisfies the evolution equations: 

duL=[{U) + ,L\ . (4.22) 

We refer the reader to appendices A,B and PH] for a complete proof of this property. It 



relies on the Jacobi identity for pseudo-differential operators. 
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Going back to the most general form ( |4.19|) for P, this result has a beautiful general- 
ization. Extending the set of KdV times to U = —qfXi-q/i for i = q + 1, . . the evolution of 
the solution to the string equation w.r.t. these variables is given by the higher KdV flows 
( [4.22|) , for i = q + 1, . . . . This remarkable fact turns out to be very useful for studying the 
physical content of the string equation ( [4.1tj| ), which now reads (see app. B.l) 

[P,Q] = 1, with P= Yl -(l + A;/g)tfc+gQt^^ . (4.23) 

fc>l; 
fc^O (mod q) 

Note that this general string equation can be derived as in the one-matrix case from an 
action principle ||3^, generalizing the result (|3.18|) , with the action functional 



S, = J2iktrQ''/'^ , (4.24) 



k=i 



(and the trace operation deflned in app. A. 3). The direct proof is tedious but the result 
also follows from the corresponding action principle for the discrete matrix model which 



we present in subsections \8. 1\ , pT; 

Remark. It is obvious that in ( [4.24|) the terms proportional to t^, m = (mod q), 
give no contribution and can thus be omitted. This seems to break the symmetry between p 
and q. If we assume p > q, however, only one term tp tr Q^^'^ is such that m = (mod p). 
It is clear that this term can be eliminated by shifting Q by an irrelevant constant. More 
generally, in app. B.2 an equivalence is established between the two actions expressed in 
terms of Q and P respectively. (See also for more on p-q duality in these models.) 



To summarize, the solutions of the string equation are not unique, but can be deformed 
in an inflnite number of directions to reach any {p, q) with p, q relatively prime. These 
directions are identifled with KdV time evolutions (see appendix |BD, but from a field 
theoretical point of view they correspond to fiows along RG trajectories between various 
critical theories, identified with (p, q) minimal conformal models coupled to gravity. Such 
trajectories are explored by adding an infinite number of relevant matter operators dressed 
by gravity to the original critical action. Therefore, the KdV time variables provide us 
with a natural definition of the operators in the theory: let the insertions of operators (pj^ 
into correlators be defined dually by differentiation w.r.t. the KdV fiows tk, 
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(4.25) 



The first q — 1 operators have been singled out in our approach, as dual to the integration 
constants introduced to rewrite the string equation. They play a special role in the (p, q) 
model picture: they correspond to the gravitational dressing of the order parameters of 
the conformal model. 

It is very useful to make an analogy with the Landau-Ginzburg (LG) picture [BH] of 



the minimal models in the particular case where \p—q\ = 1 (unitary series). From this point 
of view, the order parameters are normal ordered powers of the basic LG field, specifically 
the first q — 2 such powers in the {q + l,q) case (we exclude the identity). On the other 
hand, the LG action reads = ( Vj. )^ + + . . . + $2g-2 _^ $2g-4^ $2g-4 therefore 
the thermal operator $(1,3), which drives the {q + 1, q) model to the (g, q — I) model. We 
can see that a formal correspondence 1-^ L"^+^ exists between the LG picture and 
its gravitational counterpart in the framework of KdV fiows (with m G [0, — 2] ||30|| ). 
We can also expect subtleties for higher operators, since even the LG definition becomes 
ambiguous due to the interplay with the equation of motion 99$ oc (for which the 

corresponding fiow, generated by L^'' = Q^, is trivial). 

Scaling dimensions. It is rather easy to extract directly some information from the 
string equation ( |4.16| ) in the spherical limit. Then the string equations become a set of 



coupled algebraic equations for the coefficients Vi. The solutions are scaling functions of 
the variables tk- Giving dimension 1 to the variable x oc ti coupled to the most singular 
operator, we find by a simple counting argument: 

In particular for the string susceptibility u oc ^^-2, we find from which we 

recover the string susceptibility exponent ( |4.17D . 

Technically, once the string equation has been solved and an operator L is constructed, 
we have access to all the correlators ( [4. 25] ) by repeated use of the fiow equations ( [4.22| ). 
Consider the evolution equation 

dt^L=[(Ln + ,L]=-[{Ln-,L] . (4.27) 



Equating the coefficients of d ^ on both sides of eq. ( [4.27| ) gives 



^dt^v,.2 = d,ResL^ , (4.28) 
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where the symbol Res stands for the residue of a pseudo-differential operator (see 
app. A. 3), i.e. the coefficient of d~^ in its formal power series. Integrating once w.r.t. 
X gives an expression for the 2-point functions, 



) oc ResQ"^/'? . 



(4.29) 



This gives in the spherical approximation the scaling behavior 



{p+q+m}/{p+q-l) 



(a result consistent with the dimension of the corresponding parameter tm in ( [4.26| )). As 
we have already discussed in the one-matrix case, to compare with the KPZ result 0] for 
the gravitational dressing of conformal dimensions we must first discuss the problem of the 
identification of the identity operator and of the cosmological constant. 
Example 1: the critical Ising model (4,3) 

The Ising model has a natural realization as a two-matrix model ||39|| in which the 
two matrices represent the +/— states of an Ising spin. The two-matrix model has been 



first solved directly using recursion relations [|T9| , |2^ . A simpler derivation follows from 



considering the commutation relation ( 4.16 ) [32,36] with 



Q 



-{u,d} + -w 



1^1 \3/2 



-w 



-P = QT + \hQT = - + {«;, d} + ^u" 



6 6 



(4.30) 



(we have set t^ = 3/7). The corresponding string action (|4.24|) leads to two equations: 



1 // 3 5 

-W --UW + -t5W + t2 



0, 



1 

12 



-u 



(4) 



3 2 1 
4 8^ ^ 2 12 



(4.31) 



u 



+ 2^ +^1 



0, 



with ti oc X. We can identify the operators dual to the parameters ^1,^2, and we recover 
the identity {cpi = V), spin {(f)2 = cr, the order parameter of the model) and energy 
((^5 = e) operators from the ^1,2,5 fiows respectively: t2 is the exterior magnetic field and 
the temperature shift from the critical value. We find the scaling properties ^5 ~ [x]^/"^. 

Note that, as explained above, the flows t^ and t4 have no significance for the model. 
Note also that in the limit — > 00, the system (|4.31| ) reduces to the pure gravity equation. 
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This is a perturbative manifestation of the weU known RG flow of the Ising model to pure 
gravity in the high temperature hmit. 
Example 2: The tricritical Ising model (5,4) 

We parametrize now the two operators as (tg = 4/9) 

Q = (d"^ -u)^ + {w,d} +v , (4.32a) 
-P = gf = {d'-uff + l{wA'} + l{v,d} - , (4.326) 

where w is again a ^2 breaking field that results in coupling to a magnetic field. Again 



varying the corresponding action (|4.24| ) we find (with tg = 4/9) 



4i?3 + ^v" - ^uv + ^w^ + lts = 0, (4.33a) 
o 4 4 z 

lyjW _ -(uw)" - -uw" - -vw + -u^w - 4t2 = , (4.336) 
2 4^ ^ 4 2 4 ^ ' \ J 

16 8 8 8^ ' 

5 5 3 

- -ww" + -w'^u tsu + ti=0, (4.33c) 

4 4 2 

where a; oc ti. 

We can also perturb the tricritical Ising model in the direction of the next lower model, 
i.e. the Ising model. Instead of ( [4.326| ), we use 



where Q remains as in ( 4.32a ). The equations that follow from [P, Q] = 1 are the tricritical 



Ising equations plus tj times respectively the critical Ising equations ( [4.31| ) (with critical 
temperature = 0). (For "large" in some sense, the equations cross over to the Ising 
equations.) The most general perturbation of the tricritical Ising model also includes the 
"pure gravity" piece, and one simply has to consider P = P — |t6(d^ — u) as in ( [4.19| ). 

The genus zero equations for the tricritical Ising model are given by ignoring the 
derivatives in the tricritical equations. From these equations, we read off the scaling 
properties u ~ [a;]^/"^, v ~ [x]^^^, w ~ [x]^/^, so that t2 ~ [x]'^/^, ts ~ [xf/'^, tj ~ [x]^/"^. 
7/8 is the gravitationally dressed weight of the spin field (the (2,2) operator with undressed 
conformal weight 3/80) in the tricritical Ising model, 3/4 is the dressed weight of the 
energy operator (the (3,3) operator with undressed conformal weight 1/10), and 1/4 is 
the dressed weight of the vacancy operator (the (3,2) operator with undressed conformal 
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weight 3/5). Derivatives of the free energy with respect to the parameters ti generate 
correlation functions of the associated operators. 

More generaUy, the unitary minimal models {q,p) = (n, n + 1) flow onto each other 
when one perturbs by the $(1,3) "thermal" operator generalizing the Ising energy. In our 
language, this is transparent in the string equation 



rn+l 2n 1 j-n-1 rr 
l2n-l-t^+ 5 



1 , (4.35) 



which interpolates between the [n + l,n) and (n, n — 1) models, and identifies the field 
4'2n-i dual to t2n-i with the gravitationally dressed $(1,3) operator.lli 

4-4- Solution of the unitary (n + l,n) models on the sphere 

In this section, we carry out completely the program sketched above in the case of the 
unitary series p = n + 1, q = n. By restricting to spherical correlation functions, we ignore 
the possibly non-perturbative information contained in the string equation (which anyway 
is subtle, as we will see in sec. P), and concentrate on their leading perturbative behavior. 

This amounts to retaining in the solution L = Q^l"^ only the leading power-like be- 
havior as a function of x, or equivalently rewriting the correlators as successive powers of 
the string susceptibility u oc Vn-i- Moreover the expected IL1 parity of the solution further 
restricts the general form of L. In general, however, the spherical solution is not completely 
determined by ( [4.16|) . Instead we have an algebraic equation (polynomial of fixed degree) 



which admits several solutions for the coefficients fo, ■■, ^'^-2- To resolve this discrete am- 
biguity, it is necessary to introduce some physical constraint, as suggested by the results 
of the preceding subsection: We impose the vanishing of all the order parameter one-point 
functions on the sphere. We now describe the particular solution this determines. 
Let us introduce for m > 1 the pseudo-differential operator 

where 

* - ^ ' ^ — , (4.37) 



p J p\ 



19 



Note, however, that this is a formal relation that can fail non-perturbatively [26|. 
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and u denotes the string susceptibility u = —d"^ log Z. Then, up to terms involving two or 
more derivatives of u, we have 

[Qm,Qp]=0, (4.38) 
for all m,p. In addition the operators Qm satisfy by definition the crucial identity: 

{Qm)- = -(w/2)"^d-"^ + 0(d— 2) . (4.39) 

In it was proved that the unique solution to the string equation with vanishing 
order parameter one-point functions on the sphere takes the form P = {Qn+i)+ and 
Q = (Qn) + - The properties ( [4.38| , [4.39| ) suffice to prove that this is indeed a solution. 



Uniqueness is a little more subtle and we refer the reader to for details. 
The string equation at genus zero is 

(n + l)(w/2)" = a; , (4.40) 

giving the critical exponent 7str(?^ + l, 'n) = —1/n. Correlators are given as explained above 
by the KdV flows: 

{<p^(j)^) = 2ResL^ (4.41a) 
((/>;(/>^(/>i) = 2Res[LV,L!!^] (4.416) 
{c^ict>^ct>,ct>,) = 2 Res [[[LL^L^ILL] - [L^,[L'+,LL]]) . (4.41c) 

In the computations, we approximate the true L = Q^/" by M = qI/"" = Qi, which 
results in the systematic error 

-L^ = --iu/2Yd''-^^ + ■ ■ ■ , (4.42) 

n 

so that for k < n we can say that = —{u/2)^d~'^ + . . . . 

This suffices to calculate all the correlation functions of the order parameters, which, 
up to a field redefinition (p^ (n + l)'-"^~^-*/^"^(/)^ and a partition function renormalization 
logZ — {n/2){n + l)~^/"'logZ, can be written 



) = 5i,i(n+ l)x^+^/" , l</<2n-2 



(</'2n-i) = ]{n + l){2n-l)x^ 

4' (4.43) 

{<Pi(pm) = Si,mnlx^/'' , l<Z,m<n-l 

(j)^ (j)^) = Imr Ni^r a;-i+('+"^+'^-i)/2" , 1 < /, m, r < n - 1 . 
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We recover the KPZ scaling exponents A; = (/ — l)/2n |^ (see sec. ^) attached to each 
operator (f)^ . Note the particularly simple form of the three-point function as compared with 
the rather cumbersome flat space expressions coupling to gravity has resulted in a 
dramatic simplification of the structure. The coefficients Nimr G {0, 1} are the restrictions 
to the order parameter space of the ordinary CFT fusion rules. They arise naturally in 
the calculation as follows. 

Using the symmetry of the three-point function under the interchange of I, m, r, let 
us take r > l,m. Then the first term in ( |4.41(^) vanishes (the leading power of d is too low 



to produce any residue), and we are left with 

{<P, <P^ <PM = -2 Res ( [LV, LL]] ) • (4.44) 

For the residue to be non-zero, we must have 

I + m>r + 1 , (4.45) 

and by symmetry, the cyclical permutations of ( |4.45D hold as well. 

Comparing this with the restriction of the fusion rules of the CFT to the diagonal of 
the Kac table, an additional condition, l + m + r<2n — l, must be satisfied for the three- 
point correlator not to vanish. To derive this, let us take again some general I, m, r with the 
restrictions ( [4.45|) and permutations. We still have to evaluate (|4.44|) , where we can drop 



the -|- subscripts from L"^, L\ since their negative part does not contribute to the residue. 
The expression is further simplified by noticing that, on the sphere, the derivative d acts 
only once, so that [P,LL] = W-^[L,U_l and [L'^, [L^ L!!]] = ml L'^+^-'^[L,[L, LL]]. 
This gives the general form for the three-point correlator by symmetry, 

{(pi (f)^ (f)^) = mlr F{1 + m,r) = mlr G{1 + m + r) . 

It is now easy to compute the function G for some convenient values of I, m, r. In the 
case r = n for example, the flow is trivial so the correlator vanishes. This means that 
G vanishes whenever its argument I + m + r > 2n (recall I + m > r). This is exactly the 
relevant fusion rule. To calculate the value of G when it is non-zero, let us take l+m = r+1 
(the boundary of the fusion rules). The leading contribution to (|4.44D reads 



2[d"^, [d^ (M/2)M-n]' = 2lm {{u/2yy"d-^ , (4.46) 
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which, after integration over x, yields G and (|4.43 ) follows. 



A slight adaptation of this line of arguments yields the general A^-point function of 
order parameters on the boundary of the fusion rules, i.e. {(prm ■ ■ ■ ^mjv)' with mi + . . . + 
rriN-i = niN + N — 2. We evaluate 

{^.■■■<Pn.J = -'2 J dx Res [L^\[L^M..,L^--\L^-]...] 

= -2mi...mN J dx Res (L^^"^ [L, [L . . . [L, L!?^/miv] . . .] ) (^•4'^) 

= mi . . . TTlTV G'Ar(mAr) . 

Let us now take mi = ■ ■ ■ = niN-s = 1- Then we use the three-point function result ( |4.43| ) 
to identify Gn = so that 



mi ■ ■ ■ Tmpf 



) = mi . . . mNd^-^x'+^-^ , (4.48) 

where we recognize the KPZ scaling exponent s = Xli^i ^mi— 7str— A+2 = 2—N+{mj\[ /n) 
for the A-point function (see sec. |^). 

4-5. An alternative method for solving the {p,q) models on the sphere 

The tree level solution for the Q differential operator of the (n, n+1) models of previous 
subsection can be easily expressed in terms of a variant of the Chebychev polynomials of 
the first kind (see also |3^), as 



Q = WX(d/v^) , (4.49) 

where 

T^(2cos^) = 2cos(n^) , (4.50) 

and we have set v{x) = —u{x)/2. This can be proven simply and compactly, leading to 
interesting possible generalizations. 

At tree level, the differential operators P = P(f(x),d) and Q = Q[v{x), d) , with 
respective degrees p and q, can be rewritten 



P = yP/^ P(l, d/v^) , Q = v'l/^ Q(l, d/v^) , 
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since we can neglect the non-commutation of d and v{x) at this order. In the spherical 
(i.e. semiclassical) approximation, the commutator in the string equation ( [4.16|) can be 
replaced by Poisson brackets and thus becomes 



dP dQ _ dQ dP 
dd dv dd dv 



v' = l. 



Let P{z) = P{l,z) and Q{z) = Q{l,z) denote the two corresponding polynomials, 
with respective degrees p and q. The string equation now reads 

qP'iz) Q{z) - pQ'iz) P{z) = 2pq . (4.51) 

Suppose p > q, then the polynomial P is given by P{z) = {Qiz)P/i)^, where again the 
subscript + means we retain only the polynomial part of the formal series expansion in z. 
On the other hand, due to the differential equation ( [4.51j ), we can also write 



P = 2pQP/i{z) [ dtQ-^-P/'?(t) . (4.52) 
Jo 



To find solutions to (|4.51| ), let us try Q = Tq{z) where the Chebychev polynomial 



Tq{z) is defined in (^4.501 ). The formula (|4.52|) yields a polynomial expression for P if and 



only if p is of the form p = (2m + l)q ± 1, for some integer m, in which case we have 

Piz) = f2l^^{'^yp-2i,{z) . (4.53) 



We recover the unitary solution of previous subsection when m = 0, p = q±l. In addition, 
we see a very particular pattern emerge for the theories (g, (2m-|-l)q'±l) [jT^,5I]. A general 
study of the equation ( [4.51| ) remains to be performed. 



The resolvent in the spherical limit. A surprising and useful identity can be proven in 
the spherical limit. Let us calculate the trace of the resolvent uj{z) of one of the differential 
operators, e.g. Q: 



UJ[Z,X] 



tr(Q-;2)-i~- j\y j 



dk 1 



27r Q{ik, y) - z 

The integral over k is given by a residue. We thus find the relation 

duj{z, x) dQ(ik, y) 



dx dd 
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1, with Q{ik,y)=z. (4.54) 



The derivative of uj with respect to x is taken at z fixed. We can transform it into a 
derivative at d = i/c fixed, 

duj{z,x) dd 
dx 



duj{z, x) 



dx 



du{z, x) 



dx 
du){z, x) 



+ 



dx 



dd 

duj{z, x) dQ 
dd dx 



dQ 
dd 



Using this relation, eq. ( [4.54|) can be rewritten 

duj dQ duj dQ 



1 



dx dd dd dx 

We recognize the Poisson brackets of uj and Q. We conclude that the Poisson brackets of 
UJ + P and Q vanish. Therefore a; + P is a function of Q, i.e. of z. Comparing the analytic 
properties in z of t^; and P we conclude that they can differ only by analytic terms. We 
thus obtain the curious relation: 



5(2;, x) = — P(d, x) with Q{d,x) 



(4.55) 



The loop average and the problem of the cosmological constant. In the multimatrix 
models, a surface with a boundary of length s can be created by considering any non- 
vanishing average of a trace of the product of s matrices. In the double scaling limit, 
all matrices become a linear combination of the identity and the two matrices P and 
Q. If we assume p > q, then Q is the leading operator and in the generic situation, 
the renormalized loop average for a loop of macroscopic length s is thus proportional 
to tre~'*''°"^ with Sj-en = sa"'^'''^*'^ (see the arguments of the end of subsec. \2.3i} . In the 
spherical approximation, the dimension of Q is Q ~ 2;9/(p+9-i) Therefore the coupling 
constant coupled to the area, i.e. the cosmological constant, should scale as x'^'^/^'p^'^~^\ 
From the scaling relations ( [4.26| ), we thus conclude that the parameter = fJ. plays 

the role of the renormalized cosmological constant. In terms of /x, the partition function 
in the spherical approximation behaves as Z ~ h^^p^'^ and thus the associated string 
susceptibility exponent is 

7str(/u) = 1 - p/q for p> q . (4.56) 
The parameters tm have a dimension tm ~ ^(p+'J-"i)/2i7 fpj^g dimension of the corre- 
sponding operators (jjm is then 

= 1 - {p + q - m) /2q = {q - p + m)/2q . (4.57) 

These conclusions are consistent with analysis of the continuum Liouville equation of sec. ^ 
Note however that it is as well possible to create boundaries corresponding to the 
second operator P, and their interpretation is less clear. 
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5. Topological gravity 



We saw in subsec. how operators could be defined in the framework of the one- 
matrix model description of some two dimensional quantum gravity theories. They should 
correspond to polynomials in traces of powers of the hermitian matrix of the model, with 
leading behavior oc tr (M^"^). In a pioneering paper [^2|, Witten compared the structure 
of correlators of such operators, obtained from the KdV flows and string equation of 
one-matrix models, to his reinterpretation of Donaldson's cohomological theory in terms 
of topological field theory, referred to as topological gravity. The latter has a precise 
mathematical definition as the theory of intersection of classes (topological invariants) of 
the moduli space of (punctured) connected Riemann surfaces with arbitrary genus which 
we describe below. The Virasoro constraints in the context of these topological models 



were studied in P3| , ^ and the addition of minimally coupled matter was first provided 
in 



^ (where the topological system that corresponds to the multi-matrix models was 



presented) . In this section we shall not directly follow this historical approach, but rather 



we shall focus on the analysis of these models allowed by the later work of Kontsevich |46 



5.1. Intersection theory of the moduli space of punctured Riemann surfaces 

More precisely, let A^g,n denote the moduli space of Riemann surfaces E of genus g 
with n marked points (or punctures) xi, . . . ,Xn. The basic operators of the theory are 
built from exterior powers of the first Chern class ci(£(i)) of the holomorphic line bundle 

whose fiber is the cotangent space to E at the point Xi, considered as a complex one 
dimensional vector space. The intersection numbers attached to these operators are just 
integrals over some compactification Mg,n of M.g,n of their exterior product, and can be 
thought of as correlators, 

{cTd^aa, ■ ■ ■ aaj = [_ ci(/:(i))^^^ ■ ■ ■ ci(/:(„))^'^" . (5.1) 

These numbers are positive rationals, and vanish unless the total degree of the integrated 
form equals the dimension of the moduli space i.e. 6(7 — 6 + 271 = 2^2 di. The punctiire 
operator is P = ctq. The effect of its insertion in a correlator is not trivial, since it 
corresponds to fixing a point Xi (i.e. in generally covariant terms, one divides out by 
the reparametrization group which leaves Xi fixed). The topological definition ( |5.1| ) leads 
to a number of recursion relations between correlators at various genera, which make 
their computation a finite exercise. Actually all correlators can be expressed in terms of 
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correlators involving only the puncture operator ao- In this sense, the other operators cr^, 
n > 0, are (gravitational) "descendents" of the puncture operator. 

Due to striking resemblances to the one-matrix model correlators as defined through 
KdV fiows and string equation, Witten was lead to conjecture the equivalence between 
topological gravity and two dimensional quantum gravity in the one-matrix model form 



5.2. The Kontsevich matrix model 

Although the abovementioned recursion relations enabled computation of all the in- 
tersection numbers in principle, it seemed very difficult to prove the relation to the KdV 
framework. Dramatic progress emerged from M. Kontsevich's combinatorial treatment of 
the problem, in which, upon introducing a "cell decomposition" of the moduli space ^Ag,n, 
he was able to compute the intersection numbers by explicitly performing the integrals of 
exterior products of chern classes. Skipping details of the mathematical construction ||46|| , 
we shall write just the final answer for the "free energy" generating function F(to, ti, . . .) 
of the intersection numbers ( |5.1j ), with arbitrary genus and punctures, as follows. For a 
given correlator as in ( |5.1|) , define Ui = i^{dk, \ dk = i}, and rewrite (H'^i"'') = (11 ^d,)- 
The generating function F reads 



for arbitrary real parameters to, ^i, • • ■ ■ The result of Kontsevich allows F to be rewritten 
in the form 



where the sum is over connected fatgraphs, or ribbon graphs (i.e. a set of vertices joined 
by non-intersecting double lines, or ribbons, each line carrying a "color" index i = 1, 2, • • •, 
also summed over, as described after (|1.6D ), v{r) and |Aut(r)| denote respectively the 
number of vertices and the order of the automorphism group of the connected fatgraph F, 
and the product extends over all the double- lines in F, carrying colors The collection 
of parameters is related to the ti by 




(5.2) 




(5.3) 




(5.4) 



k>i 
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where (2z — 1)!! = (2z — 1) ■ (2z — 3) ■ ■ ■ 3 ■ 1, and (—1)!! = 1. The expression ( p73|) is of course 
very suggestive of the Feynman graphical expansion of the one hermitian matrix integral 
discussed in sec. || 

Suppose the number of colors is fixed, and as earlier take to be the size of the her- 
mitian matrix. Introduce the (real) diagonal matrix A with diagonal entries (Ai, . . . , Aat). 
Kontsevich proved that the N x N hermitian matrix integral 

/dMe*=^(^^'/'^-^^^'/2) 

has a well-defined N ^ oo limit S(A) as an asymptotic series of the variables ti = —{2i — 
1)!! tr A~^*~^, coinciding with e^^^-\ Moreover it is a so-called r-function for the KdV 
flows, i.e. its evolution in terms of the parameters ^2i+i = tr (A~^*~^)/ (2z + 1) = —ti/{2i + 
1)!! is given by the KdV flows 

Q = d^ -u u = -2dl logS(^) 

(5.6) 

and it satisfies in addition the one-matrix model string equation 



[P,Q] = 1, with p = ^-(A; + i)^2fc+i(Q"-'/'^ 

k>l 



These various statements will be proved in next subsections. This enables the identi- 
fication = ^2fc+i between the scaling constants of the one-dimensional matrix model 
and the 6''s, which are proportional to the couplings tk of the topological gravity theory. 

Note that the fat graphs involved in the summation ( [5.3| ) are connected, hence the 
"free energy" F is related to the matrix model partition function ( ^.5|) through the usual 
logarithm, which extracts only the connected fatgraph contribution to the Feynman expan- 
sion. Of course analogous disconnected correlators can be defined, as counting intersection 
numbers of the moduli space of possibly disconnected Riemann surfaces, an interesting ob- 
ject of study for mathematicians. The disconnected free energy is then just the Kontsevich 
integral S. 

Note also the important distinction between the matrix integrals considered earlier 
(e.g. (|1.8|) , ( |2.1|) or (|4.3|) ), in which a double scaling limit had to be taken (tuning the 
overall coefficient of the polynomial potential and sending oo), and the integral 

( p.5| ), in which the potential is fixed and only cubic, and only oo. Despite these 

apparent differences, Witten's conjecture is equivalent to identifying the operator content 
of the double-scaled one matrix hermitian model with that of Kontsevich's integral ( |5.5| ). 
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5.3. Computing the Kontsevich integral 



The numerator of the Kontsevich integral ( |5.5|) is best reexpressed after a change of 
variable M M - zA asi°l 



An{A) = J dMe^t-^^^'/^+^'^Z^) 



(5.7) 



up to a multiplicative factor e*'^ independent of M. The now customary integration 

over angular variables leaves us with an integral over the eigenvalues mi of M. A direct 
computation of the denominator of (^]^) , which is just a trivial Gaussian integral, enables 
us to rewrite 



N 

Siv(A) = JJ / d;UA,(m,) Yl 

i=l^ l<i<j<N 



mi — rrij 
i{X^ - Xj) 



(5.8) 



with the measure 



djjxi'm) = dm ( ^ I e 



i(m^/6+A^m/2+(iA)^/3) 



Let us use the notation 



|„0 „1 ^N-li 

tAy ^ ^A-J ^ ■ ■ ■ ^ ^Xj 



1 Xi 
1 X2 



X2 



X 



,N-1 
N 



A{x) 



1 Xn ■ 

for the Vandermonde determinant (p.2[) . By rebuilding a Vandermonde determinant out of 
integrals over the m^'s, the multiple integral ( ^.8|) can be recast in a single determinantal 
expression 



S^(A) 



(5.9) 



where the function z{X) = E]i(A = A) is the one-dimensional Kontsevich integral, and D 
denotes the differential operator A + 2^ ~ J^>^- particular, the function z{X) satisfies 
Airy's equation 

(D2-A2)z(A) =0 . (5.10) 



The function z has the well known asymptotic expansion 

z(A)=f;c,A-3\ with c, = (-l/36)^^^^"^^" 



A;=0 



(2/c)! 



(5.11) 



We summarize here various techniques appHed to the computation of the Kontsevich integral. 
For more details, see for instance especially |48|, and [p9|-^. 
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Let us define also the function ^(A) = X~^D z = ^(ifcA"'^'^ witli dk = Cfc(6/c + 1)/(1 — 6/c). 
In tlie determinant appearing in tlie numerator of (p.9[ ), we can replace D ^^z and D 2^+1^; 
respectively by X'^'^z and X'^^^^z, up to linear combinations of columns of lower rank. When 
reexpressed in terms of the variables Xi = A~^, this results in 



S,v(A) 



I JV— 1 N — 2- 

\rf,N-l rf,N-2 r 
•X" ^ »Aj ^ • • • 



Expanding both functions, we get the final result in terms of the coefficients a^n^^^^ = 



CLn ^ — dn '■ 



|™3ni+A'" — 1 ^3n2+N—2 ™3njv | 

tAy ■ fAj • • • ■ • tAy 



N 

hn(a)= j: in-'g) '" "bN-'iV-^ ",0; (5-12) 



ni,...,njv^O \i=l 



Recall now the definition of the characters of irreducible representations of GL(A^), 



indexed by Young tableaux (/i, /i, . . . , /jv), with the i row of length li {li > I2 > . . . > I 
When evaluated on a matrix X G GL(A^), the characters read 

l^h+N-l ^l2+N-2 ^InI 
Xi„...m(^) = l^jV-l' jV-2 '""^01 = '^''^[Ph+J-^(^)]l<^<J<N ' (^-IS) 

•Ay ^ tAy ^ • • • ^ tJy 

where x = diag(a;i, . . . , xat) is the diagonalized version of X. The result has been reex- 
pressed in terms of the Schur polynomials PmiX), themselves characters corresponding to 
the totally symmetric representations of GL(A^), pm{X) = Xmi-^) i^^^ boxes). 
These can also be expressed in terms of traces of powers of the matrix X , 

ni>0 , i=l,2,... i>l 



E 



and Pm = if m < 0. The definition of these characters can be extended to any ordering 
of the /i's by using the determinantal expression in terms of Schur polynomials ( |5.13| ). By 
a slight abuse of notation, we still call these functions Xi-^ «jv("^)- readily see that 
S^r of eq. (|5.12|) is expressed in terms of these generalized characters 

Siv(A)= Yl (n«S)x3n„3n„...,3n.(A-') • (5-15) 

ni>0, i=l,...,N \i=l / 

In this form, the result is obviously a function of the traces tr A"'*, i = 1, 2, . . . (see the 
definition of x ( CT , |5T4D , with X = A'^). 
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Large N limit. The first question one may ask concerns the large limit: is it well 
defined? The answer is surprisingly simple. Note first that there is a natural gradation of 
the expansion of Sat (A) in terms of the traces tr A~*, to which we assign degree i. The 
sum is over the characters X3m,...,3nN^ which have total degree S^n^. For instance the 
first few terms in ( |5.15| ) up to degree 3 read (recall that cq = do = 1 and ci = —5/24, 
di = 7/24) 



Po Pi 

P2 P3 



Sat (A) = Co + ciIpsI + codi 

_^ tr(A-3) tr(A-i)3 





Po 


Pi 


P2 




+ codoci 





Po 


Pi 


+ 




Pi 


P2 


P3 





(5.16) 



+ ■ 



24 6 

for > 3 (we omitted the remainder of the N x N determinants because they are just 
formed by some upper diagonal blocks with diagonal elements po = 1 and do not alter the 
result). Note also that no determinant of size 4 or more with as last diagonal element 
contributes to the degree 3 piece of S^r, because one would have two identical lines in 
the corresponding matrices. Note finally that the terms of the form tr (A~^) tr (A~^) have 
been cancelled out automatically, leaving us with only negative odd powers of A in the 
traces. 

More generally, it is clear that the degree 3k contribution to Sjv (A) is independent of 
the value of A^, provided > 3k: the matrix determinants which actually contribute must 
have Po as i^^ diagonal term for i > 3/c, otherwise they would have two identical lines. 
This means that as A^ grows, the contributions of given degree 3k to S7v(A) stabilize as 
soon as N > 3k. This enables to define order by order the N ^ oo limit S of Sat (A), now 
considered as a function of the traces tr A~*, as 



ni>0, i=l,2,... \i>l 



(5.17) 



where only finitely many n's are non zero for each term. 

Dependence on odd traces only. As already evident in the above example ( |5.16| ), the 
traces of even powers of A~^ cancel out of S. This is a general fact. To see why, introduce 
the variables 9i = tr(A~*)/i, and compute the derivative do^-'E{6) of S, considered as a 



function of the ^'s. From the definition ( |5.14|) , it is easy to see how the Schur polynomials 
behave under differentiation: 

d62^Pm = Pm-2i ■ 
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Let us differentiate eq. ( |5.17|) w.r.t. 92i, term by term. Each determinant yields a sum of 
terms of the form 



P3ni+1 



P3ni+2 



P3ns-s+l-2i P3ns-s+2-2i PSus-s+S-i 
P3nM-N+l P3nN-N+2 P3npf-N+3 



P3ni+N-1 
P3ns-s+N-2i 



for any N > 3k = 3^ni (the degree of the term considered). Now take N > s + 2i (and 
> 3k). In that case, the (s + 2iY^ line of this determinant is obtained from the s^^ hne 
by exchanging the roles of Us and ns+2i', but the prefactor ■ ■ ■ aii] ■ ■ ■ cinX^l'^ ■ ■ ■ is the 
same for both terms, and symmetric in the interchange s ^ s + 2i (the corresponding a's 
are both c's or d's), whereas the determinant is antisymmetric under the interchange of 
lines s and s + 2i. Therefore the total contribution vanishes, and S does not depend on 
the traces of even powers of A~^. That the step be 2i is crucial for the prefactor to be 
symmetric and not mix c's and (i's: this is why only odd powers survive. 

The compact expression ( ^.17 ) certainly leads to a straightforward computation of any 
intersection number ( ^.1|) (see for even more powerful variants, and many examples). 
Still the link with one-matrix model quantum gravity is far from transparent at this point. 



5.4- Equivalence between topological gravity and one-matrix model 

The KdV flows again. We now turn to the derivation of the KdV evolution equations 
for the large limit of the Kontsevich integral S(t.), as a function of the ti = —{2i — 
1)!! tr A~^*~^. First notice that eq. (|5.17| ) allows us to rewrite SAr(A) as a single Wronskian 



determinant 



SAr(A) 



/l 
/2 



N 



.(7V- 


-1) 


1 






-1) 


2 




.(7V- 


-1) 


N 





(5.18) 



where /' denotes de^^f = dtrh-^f and 



(5.19) 



n>0 
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This form suggests that we introduce a differential operator Ajv, whose action on a 
function / is just the normahzed (A^ + 1) x (A^ + 1) Wronskian 



AA./ = X:«;.d^-^(/) = ^4AT 



where d = do^, and we used the fact that 



/l 


r ■ 
f'l ■ 


j-iN) 

■ fr 


(5.20) 


fN 


In ■ 


AN) 
■ ■ J N 





dOiPm = Pr 



dl.Pr 



(5.21) 



Write Atv = Wn(^^ , Wn = X]r>o^^'^ '^'^^ functions Wr are just minors of the 
determinant ( |5.20| ). In particular, wq = Ihj the choice of normalization, and we find 



wi = -dg^ log Sat (A) 



(5.22) 



as a direct result of the action of de^ on the columns of the determinant ( |5.18D . Each of 
the coefficients Wr enjoys for obvious reasons the same property as Sat to be stabilized 
order by order as grows. As the limit oo of VFat, the pseudo-differential operator 

W = Wrd~'^ of degree is therefore well-defined. 

Anticipating the final result, we plan to show that the degree 1 pseudo-differential 
operator L = W(iW~^ satisfies the KdV evolution equations ( ^^ ) with Q = L^, in the 
variables 621+1 = trA~^*~^/(2i -|- 1). Let us show first that 



dOi^N — Qi,N^N — A^^d^ 



(5.23) 



with Qi^N some differential operator of order i. For the right-hand side of this equation 
to be a differential operator of degree at most A^ — 1, it is necessary that 

Q,,Ar = (AArd*A^l)+ . 



Moreover for the equation (|5.23|) to hold, it is necessary and sufficient that it be satis- 
fied when acting on A^ linearly independent functions. Choosing the basis of the /j's of 
eq. ( p.l9|) , which are by definition in the kernel of A^r, gives 



^^^(Aa.) 4- AnSI^ - Q^,N^N)fj = deA^Nfj) = 
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where we used again the property ( ^.211) to trade for dg. . This completes the proof of 
the statement ( p.23|) . 

Translating eq. ( ^.231) in terms of W and L = WdW~^, and taking N ^ oo, results 



m 



de,L=[{U) + ,L] . (5.24) 

For the same reasons as above, the coefficients Wr of W do not depend on the even traces 
^2i, which amounts to the KdV restriction that Q = is itself differential. Writing 
Q = — and using the values of t^o = 1 and wi = —dg^ logS of eq. ( |5.22|) , with 
^ oo, we finally identify 

u = -2 dl log S . 



This completes the proof of eq. ( |5.6| ). The evolution of the generating function F as a 
function of the ^'s is therefore dictated by the KdV flows. 

The string equation. Let us start from the expression ( |5.7| ) of the matrix Airy function 
An (A), numerator of the Kontsevich integral, after the change of variables M ^ M — iA. 
The string equation will appear as a particular case of equations of motion for this integral. 
We first have 

= 2 / dM ^ ^rtr(M3/6+MAV2) 

dMkk 



dM{{M^)kk + Xl) e^f(MV6+MAV2) 

On the other hand, expressing the invariance of the integral under the infinitesimal change 
of variables M ^ M + ie[X, M], with X^j = Si^k M^i, gives to first order in e. 



= {^Ma - Mkk + 2 - Xi)MkiMik 

where the first term is the contribution of the Jacobian J = 1 + ie{Mu — M^jS) + O(e^) 
and the second comes from the exponential. Collecting both equations, we deduce 

= (Afc + (Mfcfc) -2z2_^— 2 rr- 



To see why, one can restrict the expressions of fj to depend only on the odd 9's, by setting 
02s = 0. This does not alter the result ( ^.18| ), but allows writing de2s{^N fj) = de2s{^N)fj = 0, 
which implies that c^eg^Ajv = and consequently that ^eg^W = 0. 



73 



An insertion of M^k in An{A) is generated by differentiation w.r.t. A^. More precisely, 
we have 

(a^. - iK'dxJ' - 2 E JT^iK'dx, - X7'dx,))A^{A) = , (5.25) 



referred to as matrix Airy equations, generalizing the ordinary Airy equation (|5.10| ). This 
results straightforwardly in a differential equation 



I l^k ^ ^ 



(5.26) 



for the quantity 

= lim 



iv^oo / dM e-t-^ (^^V2) • 

Using Afc as expansion parameter gives an infinite set of constraints, order by order in 
Afc, of the form 

2 E A-^("^+^)l^S = 0, 

m=— 1 

with some differential operators of the A's, expressible in terms of the ^2i+i = 
trA~^*~^/(2z + 1). These operators were found to form a piece of a representation of 
the Virasoro algebra, and the above constraints are often called Virasoro constraints. Let 
us stress again that they express nothing but the equations of motion (reparametrization 
invariance) of the matrix model, and as we have seen in subsec. their appearance is 
quite general in the framework of matrix models 

Concentrate on the first constraint L-i, obtained as the leading order term in ( |5.26| ), 

L_iS = [lei - ia,, - EZiik + i)^2.+ia,,,_ J s = o . (5.27) 

Dividing by S and differentiating twice w.r.t. 9i gives 

{-ids, - EZiik + ^)02k+ido,,^,)i-2dl logs) = 1 , 

which after the trivial translation 63 ^3 — 1/3 can be recast by application of the evolution 
equations ( |5.24| ) into the familiar one- matrix model string equation ( |3.9| ): [P, Q] = 1, with 



k=l 
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This completes the proof of the equivalence between topological gravity as intersection 
theory of the moduli space of connected punctured Riemann surfaces and one hermitian 
matrix models in the double scaling limit. We stress that the equivalence at this point 
holds only at a perturbative level since we demonstrated only an identity between two 
perturbative expansions, namely that of the free energy interpolating between multicritical 
points of the one-matrix model and the generating function F for intersection numbers 



5.5. Polynomial averages and observables 

The Kontsevich integral (|5.5| ) possesses a number of other remarkable properties. 
Among those, two concern interpretations of polynomial averages taken with the Gaus- 
sian part of the Kontsevich weight. Such quantities are the natural counterparts of the 
operators usually defined in the context of matrix models, as insertions of sources in the 
Feynman expansion of the model. Typically, inserting trM'^ amounts to restricting to 
"triangulations" surrounding a given polygon with k edges. The net effect of this is in- 
terpreted in the double scaling limit as an operator creating a microscopic hole in the 
Riemann surfaces. It is therefore interesting to identify such operators in the framework 
of the Kontsevich integral. 

Gaussian polynomial averages. We consider the average 

/dMP(M)e-^-(^^V^) 
where P{M) is some polynomial of the traces of odd powersH of M. 



In ||4g], Kontsevich proved the following statement (see also |52] for an alternative 
algebraic proof). For any polynomial P of odd traces of M, the quantity Q{A~^) = (P) 
defined in eq. ( |5.28| ) is itself a polynomial of the odd traces of A~^. This defines a map 

K : n[X] rr[X] 
P^g(A-i) = (P(M)) 

from the set 7t[X] of polynomials of odd traces of X to itself. In [^, this mapping 



was constructed explicitly using a convenient basis of 7r[X], formed by generalized Schur 



The restriction to odd powers is crucial in the following. It is indirectly related to the 
underlying KdV structure induced by the Gaussian weight, which only retains traces of odd 
powers of A^^. 
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polynomials (see also [53] for another presentation). Attaching the degree i to the variables 
9i{X) = trX^/i, the polynomial Q has a total degree deg{Q) = deg(P)/2. 
Now we see that the Kontsevich integral (|5.5| ) may be expressed as 

(-1/36)' 

n=0 

Order by order in M, we therefore find an asymptotic expansion with degree 3n polynomials 
in odd traces of (the third order term is computed in eq. ( ^.161) ). This gives an 
alternative proof of the fact that (|5.5[ ) only depends on odd traces of A~^. 

Observables versus polynomial averages in the Kontsevich model. From our knowledge 
of the usual matrix models, we might expect that the actual correlators of topological 
gravity be realized as polynomial averages within the framework of the Kontsevich integral. 
The final picture is a little more involved, however, since this turns out to be true only in 
the disconnected case (i.e. where the free energy of the theory is directly the Kontsevich 
integral ( |5.5|) and not its logarithm). We have the following result, conjectured in and 
proved constructively in [^^, using the same basis as above: For any polynomial R{d0 ) of 
derivatives w.r.t. ^2i+i = tr A~^*~^/(2i + 1), there exists a polynomial P{M) of odd traces 
of M such that 

R{d0)^{e) = (P(M)e^*''(^^'/^)) 

/dMP(M)e*''(^^'/^-^^'/2) (5.29) 
= JdMe-t"- (AMV2) • 

This shows that the intersection numbers of the moduli space of possibly disconnected 
(punctured) Riemann surfaces have a representation in the Feynman diagrammatic expan- 
sion of their generating function S as insertion of polynomial sources. This renders even 
deeper the connections between the double-scaled one-matrix model and the Kontsevich 



matrix model. In particular, it was pointed out in that these polynomial representa- 
tions of observables lead to an interesting generalization of short distance operator products 
in the topological framework. 

5.6. Generalization: multi-matrix models and topological field theory 

The Kontsevich integral may be generalized [H6[ to include higher degree potentials. 



Actually one gets a r function for the q^^ reduction of the KP hierarchy by taking 

^^tr (M-fz'?+iA)'J+i| 



l>2 



^'W = . (6.30) 

/dM( 
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|^tr(M-fz'?+iA)'J+i| 



where the subscript > 2 in the numerator means we omit the constant and hnear terms in 
M in the expansion of the polynomial, and the subscript = 2 in the denominator means we 
only keep the terms quadratic in M in the expansion of the polynomial. The logarithm of 
this integral, considered as a function of the traces of negative powers of A, is the generating 
function of a set of correlators generalizing the intersection numbers (|5.1|). These form a 



topological field theory coupled to topological gravity as axiomatized by Witten ||4^,p4 
Instead of just one "primary" operator cxo and its gravitational descendents cr^, n > 0, the 
theory possesses a set of — 1 "primaries" . . . , ^^-i, the first of which is the puncture 
operator = V, and their "gravitational descendents", whose correlators are expressed 
through recursion relations in terms of correlators involving only primaries. An explicit 



realization of these theories was given in the framework of = 2 superconformal 



theories, made topological (i.e. correlators do not depend any longer on the points of 
insertion) by a twist of the stress tensor . 

Remarkably, these gravitational primaries can be identified with the q — 1 order pa- 
rameter fields of the minimal (p, q) models coupled to gravity in the framework of the 
double scaling limits oi q — 1 matrix models discussed in sec. ^. Actually the whole set 
of correlators defined through F = limj^^ao logS^"* coincide with that of the gravitation- 
ally dressed operators of the (p, q) minimal theories coupled to gravity. The proof of this 
equivalence goes essentially as in the previous q = 2 case. 

By computing the integral (|5.3C1| ) in terms of generalized GL(A^) characters (|5.13 



or 



alternatively as a Wronskian determinant, one can establish the following results. 

First of all, the N ^ oo limit of ( |5.3C1| ) is well defined since the same phenomenon 
occurs as in the q = 2 case: terms with given degree in the variables 6i = tr A~*/z stabilize 
as A^ grows (they become independent of A^ for sufficiently large A^). 

q-reduced KP flows. The integral (|5.30| ) is a r-function for the g*'^ reduction of the 



KP hierarchy: i.e. suitably defined, the differential operator Q = dfl — |ud'^ ^ + . . . with 
u = — 29|^ logS*^"^) is found to evolve with the 9i = tr A~Yz as 

This immediately implies that the generalized Kontsevich integral ( |5.3C1| ) does not depend 
on the 9qi^ z = 1, 2, . . . 
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q-string equation. Moreover, the equations of motion for the matrix model ( p.30| ) can 
be recast into Virasoro constraints, the first of which imply the q — 1 matrix model string 
equation [P, Q] = 1, as in eqs. (|iJ6| , p^ , ^^2^, with 

P= Yl -{l + k/q)e,+,{Q'/'^)+ . 
k>i- 

fc/O (mod q) 

There is therefore an exact perturbative equivalence between (p, q) models coupled 
to gravity within the framework of double-scaled q — 1 matrix models and the q^^ theory 
of topological field theory coupled to topological matter realized through the generalized 
Kontsevich integral ( |5.30|) . 

There are still many open questions regarding these generalizations. In particular, no 
generalization of the disconnected case observables as polynomial averages is known yet, 
nor any equivalent of Kontsevich's observation about the Gaussian averages in the case 
q = 2. More topological field theories are known, providing possible candidates for the 
coupling of minimal models with W-symmetry (larger symmetry algebras, including the 
Virasoro symmetry, for a review see e.g. |^^) to W-gravity (a higher tensorial form of 
ordinary gravity, which is a metric theory) , and no matrix models have been found yet to 
represent the corresponding free energies. 



6. The continuum approach: Liouville gravity 

Thus far we have solved some toy examples of 2D quantum gravity formulated as 
matrix models. The original 2D string theory with a target space comprised of a D < 1 
conformal field theory (coupled to a Liouville mode) is not guaranteed to admit such a 
formulation. In the pure gravity {D = 0), Hard Dimer [D = —22/5) |5^], Ising {D = 
1/2) Ijl^,^, and tricritical Ising cases, there is a complete line of arguments leading 
directly from the discrete string formulation of the statistical models coupled to gravity 
to a matrix model solvable by the orthogonal polynomial techniques of sections In 
general, however, the link between the matrix model solution and the continuum string is 
not necessarily proven. 

In support of this identification, however, a number of features already known from 
continuum calculations, such as scaling exponents, initially helped to identify which matrix 
model was a good candidate to describe which string theory. But experience showed that a 
limited set of exponents was not sufficient to distinguish even between solutions to different 
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matrix models (recall that the Ising model and the Hard Dimer model share the same string 
susceptibility exponent). 

In this section, we will derive general spherical correlators for the operators of the 
continuum theory (a string in D < 1 dimension), using the framework of Liouville theory. 
Comparison with the results of subsec. U-A will show a complete agreement for the unitary 
minimal series coupled to gravity, further establishing the link between 2D strings and 
double-scaled matrix models. 



6.1. Liouville gravity and conformal matter 



Starting from the original Polyakov string action (OJ.) in flat D-dimensional Euclidean 
space, 

SM{X,g) = ^g'^'daX^d.X' (6.1) 

(i = !,...,£)), the most convenient prescription to quantize this generally covariant two 
dimensional system is to fix a conformal gauge gab = ^^9ab- The system is then described 
by the Liouville mode ip and space coordinates X*, living in the background metric g. The 
gauge fixing also introduces reparametrization ghosts 6, c with spins 2,-1 respectively, 
which we henceforth omit since we are interested only in zero ghost number operators. 
The effective action is a sum of three pieces: the ghost action, the matter action ( |6.1| ) in 
the background metric {g g), and the Liouville action: 



5l = ^ 1 (^r'da^db^ - ^R<P + 2^e-+^^ . (6.2) 

In the following we will concentrate on peculiar matter theories, made of one free 
boson, with action: 



^M = ^ V9 irdaXdbX + -^RX ) . (6.3) 



This is nothing but the Feigin-Fuchs ||5^ representation of a conformal field theory with 



central charge c = 1 — 12ao^, therefore meriting the description as D = c < 1 conformal 



matter. For comparison with the results of matrix models of subsec. |^.^| , ao has to take 



the discrete values 1 / \/2n{n + 1). In general, minimal matter is obtained by taking 
rational. (Note that the above representation also requires an explicit truncation of the 
spectrum to treat only the self-contained set of states of interest.) 
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The various parameters in eqs. (|6.2| , 6.3 ) are fixed by requiring that the total action 



S = Sl + Sm + be independent of the choice of g. This can be recast into a BRST 
invariance condition on S, with the BRST charge §cT, where T is the total energy- 
momentum tensor of the system, and c the dimension —1 ghost. The result is:ii 



The center of mass of the string is described by the "tachyon" field 

Tk = e*'=^+^('=)^ , (6.5) 

where the Liouville momentum /3 is related to the matter momentum k by BRST invariance 
(i.e., that the 2D integral Tk = J Tk not depend on the choice of background metric g). 
This translates to the condition that 7^ is a dimension (1,1) operator: 

h{k-2ao)-^f3{(3 + Q) = l , (6.6) 

resulting in the mass shell condition: 

f3ik) = -^ + \k-ao\ . (6.7) 

Note that the quadratic equation (|6.6[ ) admits two branches of solutions. In the semi- 
classical approach of [^l|], the choice with the plus sign in front of the absolute value is 
chosen as follows. The Liouville coordinate can be considered as a time variable, and 
the (Wheeler-DeWitt) wavefunction reads \E' oc exp(Q(/?/2) Tk. In this language, we have 
retained only the states of positive energy E = P + The insertion of such operators 

in a correlator results in local disturbances of the surface, due to the infinite peak of the 
wave function at small geometries, f oo. The negative energy states do not correspond 
to local disturbances of the surface, and it was argued in that they cannot exist. (For 
a recent survey of the situation, see e.g. [jll|. The truncation is as alluded to after (|6.3|) .) 



As an example, the cosmological constant or "identity" operator present in the Liouville 
action, is Tk=o, with /3(0) = a+. 
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For recent review, see |11,30|. 



Analogously, the matter part of the wave function has momentum p = k — ao, so the mass- 
sheh condition ( |6.7D imphes (^P{k) + Q/2^ — (k — ao)^ = —p^ = 0, and hence describes massless 
propagation as appropriate for a "tachyon" Tk at D = 1. For D < 1 the gravitational dressing as 
well admits an interpretation as a target space on-shell condition, with E'^ —p^ = m? = {1 — D)/2A. 
By abuse of terminology, we refer to the massless mode at = 1 as the tachyon since it becomes 
tachyonic for D > 1. 
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6.2. Dressed weights 

It is now straightforward to derive some of the continuum results that compare with 
the matrix model results of sections Note that in ( |6.2| ), a_|_ is determined by the 

requirement that the physical metric he g = ge°'+'^. Geometrically, this means that the 
area of the surface is represented by / d^^v^e"+'^. is thereby determined by the 
requirement that e°'+'^ behave as a (1,1) conformal field (so that the combination d^^e"+'^ 
is conformally invariant). For the energy-momentum tensor T = —^dipdip — ^d^(p derived 
from (|6.2| ), the conformal weightH of e"+'^ is 



A(e"+^) = A(e"+^) = -ia+(a+ + Q) (6.8) 
(as in ( |6.6| ) for (3). Requiring that A(e"+'^) = A(e"+'^) = 1 determines that 

Q = -2/a+ - a+ . (6.9) 
Substituting Q = (25 — D)/3 from ( |6.4|) and solving for Q!_|_ gives 

\=.{^fW^ - ^fY^) = -% , (6.10) 



verifying ( |6.4|) . 

A useful critical exponent that can be calculated in this formalism is the string sus- 
ceptibility 7str of ( |1.12| ). We write the partition function for fixed area A as 

Z{A) = JvipVX S(^jd^^ v^e"+^ - , (6.11) 

where for convenience we group the ghost determinant and integration over moduli into 
VX. We define a string susceptibility 7str by 

Z{A) ~ A(7str-2)x/2-l ^ A ^ oo , (6.12) 

and determine 7str by a simple scaling argument. (Note that for genus zero, we have 
Z{A) ~ A'^^*-'~^.) Under the shift {p + p/a+ for p constant, the measure in ( |6.11|) does 
not change. The change in the action ( |6.2| ) comes from the term 



Recall that A is given by the leading term in the operator product expansion T{z) e°+'^'-'"^ 
Ae°'+'^ /{z — w)'^ + ... . Recall also that for a conventional energy-momentum tensor T 
— \dpdp, the conformal weight of e'^'^ is A = A = jl. 
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Substituting in (|6.11| ) and using the Gauss-Bonnet formula ^ / d^^ y^R = X together 
with the identity 6{Xx) = S{x)/\X\ gives Z{A) = e+^P^/'^'^+-P Z{e-PA). We may now 
choose e'' = A, which results in 

Z{A) = A+<5x/2a+-l Z(l) = ^(7=t.-2)x/2-l ^^^^ ^ 

and we confirm from (|6.4| ) and ( |6.10|) that 

7str ^2 + -^ = ^{D-l- ^{D-25){D-1)) . (6.13) 



Recall from the comments at the end of subsec. that minimal conformal field 
theories are specified by a pair of relatively prime integers (p, q) and have central charge 
D = Cp^q = 1 — 6{p — q)^/pq- The unitary discrete series, for example, is the subset 
specified by {p,q) = (m + l,m). After coupling to gravity, the general {p,q) model has 
critical exponent 7str = —2\p — q\/{p + q — \p — q\) (as calculated in Liouville theory, i.e. 
with respect to the area dependence of the partition function). Notice that 7str = — 1/m 
for the values D = 1 — 6/m(m + 1) ) of central charge in the unitary discrete series. Notice 
also that ( |6.13D ceases to be sensible for D > 1, an indication of a "barrier" at D = l.il 



Dressed operators / dimensions of fields 

Now we wish to determine the effective dimension of fields after coupling to gravity. 
Suppose that $o is some spinless primary field in a conformal field theory with conformal 
weight Aq = A($o) = ^($o) before coupling to gravity. The gravitational "dressing" 
can be viewed as a form of wave function renormalization that allows $o to couple to 
gravity. The dressed operator $ = e^'''$o is required to have dimension (1,1) so that it 
can be integrated over the surface E without breaking conformal invariance. (This is the 
same argument used prior to (|6.10| ) to determine «+). Recalling the formula (|6.8[ ) for the 



conformal weight of e^*^, we find that (3 is determined by the condition 

Ao-^f3{(3 + Q) = l (6.14) 



The "barrier" occurs when coupUng gravity to -D = 1 matter in the language of non-critical 
string theory, or equivalently in the case of d = 2 target space dimensions in the language of 
critical string theory. So-called non-critical strings (i.e. whose conformal anomaly is compensated 
by a Liouville mode) in D dimensions can always be reinterpreted as critical strings in d = 
D + 1 dimensions, where the Liouville mode provides the additional (interacting) dimension. 
(The converse, however, is not true since it is not always possible to gauge-fix a critical string and 
artificially disentangle the Liouville mode.) 
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(which becomes the mass sheU condition eqns. (6^, 6]7) in the string target space inter- 
pretation) . 

We may now associate a critical exponent A to the behavior of the one-point function 
of $ at fixed area A, 

F^{A) ^ ^ I V^VX ^(/d^e v^e"+^ - Jd'^^e^^ $o ~ A'-^ . (6.15) 

This definition conforms to the standard convention that A < 1 corresponds to a relevant 
operator, A = 1 to a marginal operator, and A > 1 to an irrelevant operator (and in 
particular that relevant operators tend to dominate in the infrared, i.e. large area, limit). 

To determine A, we employ the same scaling argument that led to ( |6.13|) . We shift 
(p ^ (p + p/a+ with e'^ = A on the right hand side of ( |6.15| ), to find 

4Qx/2a+-l+/3/a+ 

where the additional factor of e^^/"+ = comes from the e^*^ gravitational dressing 

of $0- The gravitational scaling dimension A defined in ( |6.15| ) thus satisfies 

A = 1 - (3/a+ . (6.16) 

Solving ( |6.14| ) for (3 with the same branch used in ( |6.10| ), 

/3 = -ig + -^ig2 _ 2 + 2Ao = (V25 -D-^l-D + 24Ao) (6.17) 



(for which —(3 < Q/2, and /3 ^ as D ^ — oo). Finally we substitute the above result for 
(3 and the value ( |6.10|) for into (|6.16|) , and findil 



V25-D- 

Assurance that the procedure of this section for identifying operators in these the- 
ories is consistent comes from calculations of the toroidal partition functions which 
essentially just count the states. 



We can also substitute /3 = a+(l — A) from (|6.16D into (|6.14D and use — |a+(a+ + Q) = 1 
(from before ( |6.10| )) to rederive the result A — Aq = A(l — A)q+^/2 for the difference between 
the "dressed weight" A and the bare weight Aq [@. 
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6.3. Tachyon amplitudes 

We shall now evaluate tachyonic amplitudes of the form 



A{ki,k2,...,kN) = {Tk,Tk,...Tk^) , (6.19) 

where the expectation value is taken with respect to the total action S. In the rational 
ao^ case, the spectrum of /c's in the CFT can be restricted to a finite set of degenerate 
representations of the Virasoro algebra, closed under operator product. Following the 



flat space correlators for these minimal CFT are computed by inserting a suitable number 



of screening operators of dimension 1, exp{ia±X), where a± = ao ± y/oo^+l = Tct±- 
We see that these are pure matter operators, since their Liouville momentum vanishes. In 
principle the computation of the most general amplitude ( |6.19| ) involves the insertion of 
say n T^^'s and m Ta_ 's. 

The flrst step in computing ( |6.19D consists of the integration over the zero modes 
of the flelds. Splitting {p = {p^ + and X = Xq + X, with f \fg(p = f \fgX = 0, and 
integrating over ip^ and along the real line, gives 



(6.20) 

/i=0 



where the scaling factor s is deflned by 

TV 



a+s + ^^{ki) = -Q{l-h) . (6.21) 



i=l 



h denotes the genus of the world sheet surface, and the above average is performed over 
the non-zero modes of the fields with the free action 

^=^j V^^"'(9aW + a„lc>,l) (6.22) 

(from now on, we drop the hat symbol on the free fields). The integration over the zero 
mode of X yields the electric neutrality condition, including the screening operators: 



N 

na^ + ma- + ki = 2ao(l — h). (6.23) 

i=l 



What has been gained in this process? First of all, we have reobtained the KPZ 
scaling exponent s for the dependence on the cosmological constant fi of the correlators as 
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in the preceding subsection. Namely to each tachyon operator T^, we associate a scahng 
dimension Aj. = 1 — (3{k)/a^, and on genus h surfaces the bare partition function scales 
like iJ.'^~^^'^'- , with the string susceptibility exponent 7sfj.'* = 2 + Q{1 — h)/a+. Notice that 
the screening operators of the matter sector do not contribute. For the (p, q) minimal CFT 
coupled to gravity (cto^ = {p ~ /'^PQ)^ we find 

litKp^ Q) = 7str(p, q) + h{2 - 7,tr) , (6.24) 
where the genus zero string susceptibility exponent reads 

7str(p, g) = 2 + = - , , (6.25) 

q;+ P + q-\p- q\ 



confirming eq. (|6.13|) . Recall that the linear behavior of the string susceptibility exponent 
( |6.24|) is the basis of the double scaling limit procedure of subsec. |j.4 

Note also that this result is in general different from that of the matrix models 
eq. ( [4.171) (in the non-unitary case \p — q\ 7^ 1), because the Liouville calculation al- 
ways selects a definition of the "cosmological constant" n as the coupling to the "dressed 
identity" operator, different from that of the matrix models (in which the scaling variable 
couples to the dressed lowest dimension operator, see eqs. ( [4.561 , ^4.57] )). These two defini- 



tions coincide only in the unitary case |p — q| = 1, where we will be able to compare the 
Liouville results directly to those of the matrix model. 

We can moreover find the dressed scaling dimensions of the tachyon operators T^^ ^ , 
where each momentum is quantized in units of a± according to the number of matter 
screening operators of each kind, 

kr,s = ^(1 - r)a+ + ^(1 - s)a- 

{{r + s - 2)\p - q\ + {s - r){p + q)) 



2.J2pq 

^ (^(r + g) - l)\p-q\ + \{s - r){p + q) 
p-\- q-\p - q\ 

Taking into account the restrictions over r, s, the last expression can be rewritten 



-\p-q 


+ 


ps — 


qr\ 


p + q 


-\p-q 





^r,s = '1 „i ' , (6.26) 
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with l<s<q — l,l<r<p — l, and with the symmetry = Aq-s,p-r- These coincide 
with the matrix model result for unitary theories, because the identity r = s = 1 is then 
the operator of minimal dimension, while for non-unitary theories the minimal operator 
which corresponds to r, s such that \ps — qr\ = 1 differs from the identity operator. 

Beyond this, we would like to compute the factor of in the general amplitude. Our 
strategy is as follows: 



(i) We perform the computation for non negative integer s, where (|6.20| ) becomes a simple 
free field correlator, with s insertions of the tachyon Tq, integrated over all positions of 
the operators. Although tractable in principle, the free correlators have an ugly form 
for genus h > 1, and it is not evident how to perform explicitly the integrations over 
the positions of the tachyons, and the moduli of the surface. We will only present the 
spherical results {h = 0) here. This last integral is particularly well-behaved only in 
the D < 1 case: the kinematics allow for a finite region of convergence in momentum 



space. We will calculate ( |6.20|) in this region, where the result is a simple polynomial 



m s. 

(ii) analytically continue this result to arbitrary (real) Here we have to resort to 
a more physical argument: the polynomial amplitudes above will be interpreted as 
locally described by a two dimensional effective field theory, which for large momenta 
gives an algebraic growth of the amplitudes. Requiring that all amplitudes be poly- 
nomial will fix them uniquely. 

(iii) knowing the expressions for A{ki, k2, . . . , kjy) for arbitrary s and ki in the conver- 
gence domain of the integral, we still have to analytically continue them to arbitrary 
momenta. This last step is not necessary in the D < 1 case, since the momenta even- 
tually take discrete values inside the convergence domain. We will comment later on 
the D = 1 case, where this last step becomes crucial. 

An alternative (and sometimes more powerful) approach to these calculations is to 
use the ground ring structure , as implemented in ||6^ . 



28 



As originally suggested in |6S]. In [|64| , it is argued that this analytic continuation agrees 



with the asymptotic behavior for complex s. 
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6.4- Three-point functions 

We first want to evaluate the three-point function without screening: 

^(^1,^2,^3) = i-nf (^)'r(-.)(T,,(0)T,,(oo)T,3(l)To^> , (6.27) 

where we have used the SL(2,(C) invariance to fix the positions of three tachyons, and 
rescaled the partition function by a constant factor for technical convenience. The momenta 
are subject to eqs. ( |6.21 , |6.23 ) with m = n = h = 0. Performing first the free field Wick 
contractions, with the propagators {(j){z) (j){0)) = {X{x) X{0)) = — log {z]"^, we are left with 
the integral 



{T,,T,,T,,T^) = 11 1 d^wjlw^ni - wj\'^ n 



\Wi — Wj 



|4p 



j=l l<i<j<s 

X(a,/3;p) , 



(6.28) 



where 

a = -a+(3{ki) , (3 = -a+(3{k^) , p = ±- . (6.29) 

The integral ( 6.28| ) was first computed by Selberg by analytic continuation from 
integer p. Introducing the function A(x) = r(x)/r(l — x), the result reads 

J(a,/3;p) = s! (7rA(-p))^ 

W (6.30) 
■ JJ A((z + l)/9) A(l + (3 + ip) A(-l -a-(3-{s + i- l)p) . 

With no loss of generality, we can take /si > cio, ^2 > cto and ks < ao- This enables us to 
solve for ks using (|6.21J ) and ( |6.23| ), or equivalently solve for /3, 

U(l-.) forao>0 
^ \ -1 - for ao < . ^ ^ 

This expression for P implies many cancellations in ( |6.30|) , and we find 

ao > A{ki,k2,ks) = 

ao<0 A{ki,k2, ks) = -7rA(-s) (pA(-p))' JJ -7rA( 



2 

—n/\(mi 

i=l 



(6.32) 



where 

B'^ - k'^ 

m, = ^^-^ . (6.33) 
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In the latter case, ao < 0, we have ma = —s, so the result (|6.32| ) can be put in a more 
symmetric form, 

3 

A{h, k2, ks) = {^A{-p)Yl[ -7rA(m,) . (6.34) 

A few remarks are in order: 

(i) The apparent divergence of the ao < result due to poles of the F-function at negative 
integers should be understood as a finite contribution for fixed area; the area A and 
the cosmological constant n are related by a Laplace transformation 

1 

n'A{-s) = - dAe-^'^^A-'-^ (6.35) 
■s' Jo 

(recall that the KPZ scaling is derived for fixed area amplitudes). 

(ii) This explains why we seem to get two qualitatively different results for the two signs 
of ao. Taking a closer look at (|6.34| ), we see that all the factors are finite for ao > 0, 
whereas ma = —s yields a divergence for ao < 0. Therefore ( |6.34| ) gives a vanishing 
finite area amplitude when ao > and a finite one for ao < 0. Both answers are 
correct, but the ao > vanishing is just an artifact of the resonance condition imposed 
on the momenta. ( |6.34|) is thus the general result, independent of the sign of ao. 



The f actor ized form of the result suggests defining "renormalized" tachyon operators: 



-7rA((/32-A;2)/2) 



(6.36) 



(This includes the cosmological constant operator, and we redefine /i as the coefficient of 
To in the Liouville action (|6.2|).) In terms of the operators (|6.36| ), the amplitude is simply 

{n,nM = ii' . (6.37) 

The situation here seems to be much better than in ordinary string theory, where tachyon 
amplitudes contain poles corresponding to the massive modes of the string and integrating 
them out leads to a highly non-local effective action. Here the tachyon interacts with an 
infinite set of massive modes, existing only at discrete values of the momentum (the poles 
of the renormalization factor A(m)), and integrating out these modes has the mild effect 
of renormalizing the tachyon field, therefore described by a 2D effective field theory. This 
suggestion is compatible with higher amplitudes growing algebraically with the momenta, 
i.e. polynomial in momenta. Under this assumption, (|6.37|) holds trivially for arbitrary s 
(1 is the analytic continuation of the polynomial P{s) = 1). 
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If we include arbitrary numbers n and m of screening operators, ( |6.3C1| ) is replaced 
by generalizations of the Selberg integral computed by Dotsenko and Fateev For 
the special kinematics chosen above, correlators simplify owing to many cancellations that 
occur. After some algebra, the result is 

3 

A^n{ki,h,h) = (^A(-p))'(-7rA(-p+))"(-7rA(-p_))"^ JJ-ttAK) , (6.38) 

i=l 



where p± = aj_/2. Performing the field redefinitions (|6.36|) also for the screening operators 



Ta^, we find the same result (|6.37|) , trivially extended to arbitrary s. 



This answer is not quite final, however, since a subtlety arises in the discussion of the 
selection rules for the three-point amplitudes. In the present case, we have ki = k^^^Si = 
i(l - ri)a+ + i(l - Si)a- for i = 1, 2, and ks = kj^^ = ^{1 + rs)a+ + i(l + Ss)a-, 
with a?_ = 2p/q, p < q. We have introduced the Kac indices 1 < ri < 1 < Si < q, 
such that Viq > Sip. The only constraint on the integers r^, Si arises from ( |6.23| ), and 
amounts to ri + r2 > r^, si + S2 > S3, and ^ = 1 mod 2, ^ = 1 mod 2. This cannot 
be the only constraint on the Kac indices for the three-point coupling to be non-zero, 
because it would violate the symmetry under permutation of the tachyons. Moreover we 
seem to have lost the truncation of the CFT fusion rules: 'Y^Vi < 2p, "^Si < 2q. In fact 
the correct result must be a consistent projection of (|6.37|) . In the fiat space CFT, this 
projection is given by the factorization properties of the four-point functions onto three- 
point correlators. This result does not arise directly from the Feigin-Fuchs integrals for 
the three-point correlators, which in addition have to be truncated by allowing any fiip 
of two vertex operators 14,s(= expikr^s^)—^ ^-r,g-s- We expect a similar phenomenon 
here, therefore leading to the same fusion rules as in the ordinary CFT case (recall that 
those arose naturally in the matrix model solution of sec. ^). The final result should then 
read: 

with the CFT fusion numbers N(r. G {0, 1}, and reproduces the KdV results of sub- 
sec. \(.4[ apart from different field normalization factors. 
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6.5. N -point functions 

We now turn to the computation of arbitrary A^-point functions of the tachyon field. 
We beheve the screening to be just a decoration of the amphtude, responsible only for a 
global renormalization factor and for implementing the CFT fusion rules. We will therefore 
concentrate on the A^-point functions without screening. The conservation laws ( |6.21 , 6.23| ) 
now read: 

N 



^^ki = 2ao 

=1 

N 

El V 
\ki - ao\ = {-^ - 1)Q 



(6-40) 

N 



1=1 

Fixing three tachyon positions by SL{2,(C) invariance, the A^-point amplitude reads 
for integer s: 



s ' r ^ r 

A{h,k2,...,kM)=-n' {^yn-s)ll / d^Wall / d^- 

s N 

{o)%, (oo)r,3 (1) n Mwa) n 



AT 

n 

a=l i=4 

and the free field integral reads: 

TV s „ N 



(6.41) 



^ -A 1*^ ^ A 



1=4 a=l 1=4 

a<b 4<i<j i,a 



3 



(6.42) 



25i 

I "^a I 5 



with ^ 

5,^-2a+l3{h) , 9,,j=hkj~P{h)P{kj) , p=-^-. (6.43) 

We now restrict to the (A^ — 1,1) kinematics (it is possible to show that the amplitude 
vanishes for any other choice (n, m) just as in the three-point case): 

/ci > ao , /c2 > tto 5 ■ ■ • 1 k^-i > tto , and k^ < aQ ^ (6.44) 
which enables us to solve for k^ using ( |6.4C1| ), 

^ + ^-3 1 

kN = a+ + -a- ■ (6.45) 
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It is natural to trade the momenta ki for the variables = {^l3(kiY — A;f)/2, in terms of 
which we have 

Si = p — rrii , i < N 
6n = ~1- {N + s-3)p 

(6.46) 

9ij = -rui -ruj , i < j < N 
9,,N = -I + {N + s - 3)m, . 

The complexity of the integral (|6.42[) resides in its pole structure.ii As in general 
Veneziano string amplitudes, poles arise when intermediary channels have momenta in 
certain discrete sets. From the point of view of the integral, those correspond to integration 
regions where some number of Zi approach each other. We can make this more precise in 
an example: suppose we want to investigate the behavior of ( p.42| ) when 24, Z5, . . ., Zp 
approach simultaneously. Then we can redefine 

Z4 = e, ^5 = e 2/5 , . . . , Zp = eyp . (6.47) 

Performing the integral over e in a small disk, a pole arises at the intermediate state of 
energy E = ^ + J2^=i(^^ momentum k = J2^=4^i^ such that E"^ — [k — ao)^ = 2/, 
where / a non-negative integer labelling the excitations (Z = is a tachyon intermediate 
state, / = 1 a graviton, etc.). What are the residues at these poles? It is easy to see 
that the answer is given by the factorization of the amplitudes. For instance the tachyon 
residue is: 

(T.,T,, . . .T,J oc E^-{k-aoY ' ^ ^ ^ 

and a straightforward generalization holds for non-zero /. 

The miracle which occurs in our case is the vanishing of most of these residues at in- 
termediary states, leaving us with almost completely factorized n-point amplitudes. Con- 
centrating first on the tachyonic poles (/ = 0), we are going to show that at least one 
of the two pieces of the residue ( |6.48| ) vanishes identically, due to the fact that the in- 
termediate state is a special state in the wrong branch. Namely, one has in this case 
m{k) = {l3{kY — /c^)/2 = p — 3, and P = —Q/2 — \k — ao\. Note that in this case the 
momentum takes the value k = (a_ — pa_|_)/2, and therefore corresponds, from the matter 



29 



In cases where the pole structure is insufficient to resolve ambiguities, one can refer as well 



to the ground ring structure possessed by the operators (see [35,& 
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CFT point of view, to a representation with Kac indices (0,p — 1), outside of the minimal 
table. 

We proceed by induction. Let us suppose that all the M-point functions with M < 
N — 1 have the form: 

M 

{Tk.Tk, . . .TkJ = P{ki,k2, . . . ,kM)l[Mm,) , (6.49) 

1=1 

with some polynomial P in the momenta. This is certainly satisfied for all M < 3. Plugging 
this form into the residue of ( |6.48D , we find that the first bracket is finite, whereas the second 
one vanishes due to the factor A(m(/c)) = A(p — 3) = 0, and therefore the residue vanishes. 
More generally, one can consider the situation in which any subset of 2:4, ... , zn-i is taken 
simultaneously to 0. The corresponding residue can be shown to vanish as well using the 
recursion hypothesis. More interesting are the poles obtained from integration regions 
where a subset of ^4, . . . , zn-i, together with zn, approach 0. In that case, one finds that 
the second bracket is finite, but the first one vanishes, due to the factor A(^m(2ao — k)) = 0, 
enforced by the kinematics. We are therefore left only with poles where zn approaches 
other vertices. But as a consequence of the relation ( |6.45|) , these turn out to be factorized 
poles in the individual momenta, involving s only. The residues of the latter poles do not 
vanish, due to the simultaneous presence of a special state of the wrong branch (which 
would in principle induce the vanishing of the residue), and of a special state of the right 
branch (which would give a divergent contribution), whose competing effects cancel each 
other. 

The poles at excited intermediary states (/ > 0) have the same fate. The latter can 
be seen as gravitational descendents (Liouville dressing of the Virasoro descendent fields) 
of the (wrong branch) tachyons. Typically such an operator is built by multiplying an 
operator ^ with (3 in the wrong branch, by a differential polynomial of the matter field 
X . Starting from a momentum k = kr^s of the minimal conformal grid, one obtains the 
first excited state at level rs, etc. The point is again that any insertion of such an operator 
in a generic correlator yields a vanishing result. The proof proceeds by induction, and we 
leave it as an exercise to the reader. This shows, just as in the / = case, that all residues 
vanish except for a set of factorized poles in the individual momenta. 

Where are these last poles located? The poles in mi come from the region where zn 
approaches 0. The poles in 1712 come from — 00, but can be recast into poles in mi, 
using the kinematic relation: 

N 

= l+ps. (6.50) 

i=l 
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In principle poles in mi occur whenever = — I = 1,2,..., but it is easy to see 

that among those only the ones corresponding to mi = — n, n = 1,2, . . ., have a non- 
vanishing residue: we know already that the locations of the poles only depend on s, so 
that we can take /cs, /C4, . . . , kN-i — > and be left with the three-point function ( |6.34|) , 
differentiated — 3 times w.r.t. /i (corresponding to — 3 insertions of the cosmological 
constant operator To), from which the poles in mi are obvious. Similarly, we find poles at 
7712 = —P: p = 1, 2, . . ., or equivalently at mi = 1 + ps — XliLs^ P = I5 2, . . . , . 
We now form the ratio 



(Tfc^ ■ ■ ■ TfcAf ) 

nil 



/s(mi;m3, . . . , mAf-i) = ' " W . (6.51) 



As a function of mi, fs can have poles only at the zeros of the denominator, i.e. at 
mi = 1712 = n, /, n = 1, 2, . . . , . To show that the numerator also vanishes at these 
points, we apply the recursion hypothesis ( |6.49|) to the residue of the individual pole at 
rriN-i = —n, which vanishes for mi = / or m2 = n. Therefore fs is an entire function of 
mi. Moreover we find that fs is bounded when |mi| 00. For large |mi|, the integral 
( |6.42|) is dominated by the region where all the points are close to 1. To blow up this region, 
we perform the change of variables Zi = exp(xi/mi), Wa = exp(tt;a/mi), and estimate the 
integral. The result is simply lim^^^oo fs = const. Repeating the argument for the other 
rrii, we find that fs is only a function of and s. Again we can send /C3, . . . , /cat-i 0, 
and read the result from the three-point function (|6.34|) , differentiated — 3 times w.r.t. 
11. This gives 

N 

{n,n, ...n^) = {-nA{-p)y [n -nA{m,)] (a^)^-3^^+^-3 , (6.52) 

i=l 

or after the redefinitions ( |6.36| ) of the tachyon and cosmological constant operators, 

(T,,T,,...T,„) = (a^)^-3^^+^-3 . (6.53) 



Assuming, as explained above, that the result is polynomial in the momenta, (|6.53|) still 
holds for arbitrary s, not necessarily integer. 

This concludes our calculation of the A^-point amplitudes on the sphere without screen- 



ing, showing a perfect agreement with the matrix model (KdV) results of subsec. 
(eqn. (CT))- 
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The general amplitudes of the rational case involve insertions of screening operators, 
enforcing the fusion rules of the CFT, but we believe the form of ( |6.53[ ) is unaffected by 
these decorations. 

As indicated earlier, we have to be careful with the problems of convergence of the 
defining integral ( |6.42| ). It is easy to see that the convergence domain is simply: rrii > 0, 
i = l,2,...,A^ — 1, with the kinematical constraint ( |6.50| ). Therefore, taking the momenta 
at discrete values corresponding to the conformal grid ki = k(ri,si)j i = 1,2, ...,A^ — 
1 and /cat = /c(p_rj^^g_s^) does not violate convergence, and we do not have to worry 
about analytic continuation of the results in momentum space. However, in the interesting 
limiting case ckq ^ (c ^ 1), all real values of the momenta are permitted, and we have 
to worry about this problem. The structure of the amplitudes turns out to be much richer 
in that case, although identical in the (A^ — 1, 1) kinematical regime (see [BS for details). 



7. Large order behavior and Borel summability 

In this section we show how one can determine the large order behavior of the topolog- 
ical expansion of the d <1 models by a straightforward analysis of the differential equation 
satisfied by their partition functions. After recalling some standard facts about divergent 
series, we discuss in detail the simplest case of pure gravity, and sketch the generalization 
to other cases. 



7.1. Divergent series and Borel transforms 

We start by recalling some standard features of divergent series, Borel summability, 
and summation methods (see, for example, pp. 840-842 of for a recent treatment 
with physical applications). Consider a function f{w), analytic in some sector S (say 
|Argi(;| < a/2, \w\ < \wo\) in which it has an asymptotic expansion 

oo 


This means that the series diverges for all no n- vanishing w, but in 5" there is a bound of 
the form 

N 
k=0 



< Cn+1 for aU N , (7.2) 
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and for definiteness we assume that Cat = M {(3N)l. Though the series diverges, it 
can be used to estimate f{w) for \w\ smaU by taking the truncation of the series ( |7.1[ ) at a 
value of that minimizes the bound (|7.2|) . This gives the best possible estimate of /(w), 
with a finite error e{w) = minj^^v} Cn \w^\ ~ exp —{A/\w\)^^/^. An asymptotic series does 
not in general define a unique function since we can always add to it any function analytic 
and smaller than e{w) in the sector 5". 

When the angle a defining the sector 5" above satisfies a > rcf], however, a classical 
theorem of analytic functions applies to show that a function analytic in S and bounded 
there by s{w) must vanish identically. This is the only case in which the asymptotic 
series defines a unique function /(tf), and for which there exist methods to reconstruct 
the function from the series. One such method is based on the Borel transform Bf{w) of 
f{w), defined from the expansion ( [7.1| ) by 

oo oo „ 







According to the assumptions following ( |7.1[ ), we have \fk/{(3k)\\ < M A~'^ , so Bf{w) is 
analytic at least in a circle of radius A and uniquely defined by the series. Then the integral 
representation 

/"OO 

f{w) = dt e-^Bf{wt^) (7.4) 
Jo 

converges in the sector |Argi(;| < a/2 for \w\ small enough, and yields the unique function 
which has the asymptotic expansion ( [7.1|) in the domain S. 

In general, the function Bf{w) may have poles and cuts running from its singularities 
off to infinity on the complex plane. If there is a singularity on the positive real axis, we 
say that the original series for f{w) is not Borel summable since we cannot run the contour 
in (|7.4|) along the real axis. Cuts in Bf{w) are indicative of possible "non-perturbative" 
effects (i.e. exponential in an inverse string coupling, the notation of 

( p.l3| )), and generally the choice of contour in ( [7.4| ) from the origin to oo in the Ret > 
half-plane refiects possible "non-perturbative" ambiguities. 

We now recall how the large order behavior of f{w) may be extracted from such 
non-perturbative behavior. The coefficients of the series defining Bf{w) satisfy 

1 f ds 1 

^'=2mtc7^^f^'^ .^oo 2m I ^'T^^'^'^^f^'^ ' (^-^^ 
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where the contour C encloses the origin. For k large, the behavior of is related to values 
of Bf{s) near the point on the contour where |s| is minimal, so by deforming the contour 
to run along the cuts of Bf{s) to infinity, we see that the integral above is dominated by 
the discontinuity of Bf{s) along the cut corresponding to the singularity r closest to the 
origin. After Borel transformation, it follows that the large order behavior of the original 
series is given by 

oc / die"*/ -^discS/(st^) oc / / die"* disc (7.6) 

Jo Jr/t ^ Jo S Jr 



k — >oo 



where f±{s) are the Borel transforms corresponding to integrations in (|7.4|) on opposite 
sides of the cut. 

In the cases of interest to follow here, both f± will satisfy the same differential equa- 
tion, and the (exponentially small) difference 

/•OO 

e{w) = f+{w) - f-{w) = dte-^ disc Bf{wt^) (7.7) 

J r 

will be determined by the corresponding linearized equation. Knowledge of the leading 
behavior of e{w) can be used to infer a great deal about f{w) and Bf{s). For example, 
when e{w) has leading behavior 

.H-.-V^e-^-^/*")'", (7.8) 

we find from (|7]^) that fk ~ r(/3/c + b) for k large. Moreover we see from ( [7.7| ) that the 
above behavior for €{w) results from the singular behavior Bf{s) ~ {l — s/A)~^ near s = A, 
where A is the singularity nearest the origin in the Borel plane. For large k this means 
that bk ~ A~'^k^~^ and hence we have the refined estimate fk ~ {(3k)l A~''k^~^ . This 
is typical of large order behavior of perturbation theory in quantum mechanics and field 
theory, where A is a classical instanton action.ll§ In what follows, we can now bypass the 
intermediate steps and use eqs. (|7.6| - [7.8|) directly to determine the asymptotics of fk and 
the locations of singularities of Bf{s). 



30 



As stressed by Shenker |70], however, the value of (3 provides an important distinction between 



string theory and field theory. In field theory, /3 = 1 and the non-perturbative effects in the 
exponential in ( |7.8| ) go as the inverse loop coupling, = In string theory, on the other 

hand, we shall see that (3 = 2, leading to much larger non-perturbative effects in the exponential 
that go as the inverse square root of the loop coupling, w^^^'^ = 1/k. 
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1.2. Pure gravity 

For pure gravity, the differential equation satisfied by tlie second derivative of the 
partition function is (( |3.13| ) with 1 = 2, after suitable rescaling)il 

u'{z)-^u"iz)=z . (7.9) 

If u{z) has an asymptotic expansion for z large, it satisfies u{z) = ±y/z + O {z~'^). The 
solution that corresponds to pure gravity has a z large expansion of the form 

u{z) = z-s'^/^) , (7.10) 

fc=i 
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where the Uk are all positive.! 

To determine the large order behavior of the expansion we argue as in (|7.3| - [7.6| ). We 
consider the Borel transform of the expansion, defined by 

oo 

in which P is chosen so that the series ( |7.11[ ) is convergent in a circle of finite radius. Then 
a solution (in general complex) to eq. ( [7.9| ) is obtained from the integral 



uiz] 



jdte-'B{tPz-'>/^)) , (7.12) 



provided there exists a suitable contour of integration from the origin to infinity in the 
Ret > half-plane on which the integral converges. 

The functions u±{z) defined respectively by integration in ( [7.12| ) above and below the 
cut both satisfy eq. (|7.9|). For z large, their difference e = uj^ — u- is exponentially small 
compared to their average uq = {u+ + u-)/2. e{z) is therefore a solution of the equation 
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We exchange x for z in what follows since we continue to the complex plane. Our normal- 



ization in (7^) corresponds to a matrix model with an even potential; for an odd potential the 



second term is instead —hu". 

6 

"^^ The first term, i.e. the contribution from the sphere, is dominated by a regular part which 
has opposite sign. This is removed by taking an additional derivative of u, giving a series all of 



whose terms have the same sign — negative in the conventions of ( [7.9| ). The other solution, with 
leading term — z^^'^, has an expansion with alternating sign which is presumably Borel summable, 
but not physically relevant. 
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obtained by linearizing (|7.9| ) . Taking the differences of the equations (|7.9| ) satisfied by 
and w_, we find that e satisfies 



e"{z) - 6uo{z) e{z) = , (7.13a) 

where uq is determined by 

uliz)-lu'^{z) + e'{z)=z (7.136) 



(and e^{z) in ( [7.13fc|) can be ignored to leading order in large z). To leading order, the 



function e is also proportional to the difference between any Borel sum of the series and the 
exact non-perturbative solution of the differential equation (up to even smaller exponential 
corrections corresponding to multi-instanton like effects). 

Eq. ( |7.13a| ) can easily be solved by the WKB method for z large. Substituting the 



ansatz e'/e = ru^^ + 6wq/wo, we find = 6 and b = —1/4. Dividing by the leading term 
Uq ~ 2^/^ to remove the overall factor z^^"^ in ( [7.10|) gives 



e{z)_ 



(x z ' e 



5/8^ 5 



.±v/6^5/4 



(1 + ■•■)■ (7.14) 



In terms of the expansion parameter (string loop coupling) — z the ratio (|7] 
reads 

The above solution is valid for z large, i.e. n small, so we may apply ( [TBI ) to find that the 
large order behavior in (|7.1C1| ) is given by 

- (j7i)"r(2*-^). (7.16) 

(The constant of proportionality in the above cannot be determined by this method.) The 
asymptotic T(2k — i) behavior is a slight refinement of the {2k)l behavior determined in 

mi- 

From the discussion following (|7.8|), we can see directly from ( [7.15| ) that (3 = 2, and the 
reality of has implied a singularity on the real axis in the Borel plane. This obstruction 



to Borel summability is consistent with the large order behavior in (|7.16|) , in which all 
terms have the same sign. 
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Remarks. Note that other singularities of the Borel transform are related to higher 
order corrections in the expansion of e. Because the coupled equations ( [7.13a, 6| ) for uq and 
e have a well-defined parity in e, the exponential behavior of successive non-perturbative 
corrections to e will be of the form exp(— n C^^/"^), where C = (4/5)-\/6 and n is an odd 
integer. This follows from iterating ( [7.13a, ex.p{—nCz^^^) terms in e in ( [7.13fc| ), with 
n odd (including the leading exp(— C^^/^) piece), result only in exp{—mCz^/^) terms in 
Uq with m even, and vice versa, in ( [7.13a[) . 

In [28 1, it is confirmed that the exponential in ([7.14[ ) coincides with the action for a 



single eigenvalue climbing to the top of the barrier in the matrix model potential, allowing 
us to interpret the exponential piece of the solution to ( [7.9[) as an instanton effect. In 



subsec. [7. ?[ here, we shall reproduce and generalize this observation to arbitrary one- matrix 
models. 

7.3. Ising / Yang-Lee 



We now consider the fourth order differential equations [p!9[ , p4[ , [34 for the Yang-Lee 
edge singularity and the critical Ising model. After suitable rescaling (different for the two 
cases, ( [3.13[ ) with / = 3 and ( [1.31[ ) with t2 = ^5 = 0), the equations are written 

- uu" ~ ^u'"^ + au^^^ = z , (7.17) 

where = jq^ ^ respectively for Yang-Lee and critical Ising. 
The asymptotic expansion takes the form 



u{z)=z'/'{l + J2^,z-''/') 
k=l 



Substituting u = uq + e (where now uq ~ z^^^), we find that the linearized equation for 
the discontinuity e{z), at leading order for z large, reads 

(4) I o 2/3 1/3 // -'- -2/3 / n 

ae^ ' + 6z ' e — z ' e — -z ' e = . 

O 

1 /2 

The ansatz e'/e = tuq + Buq/uq now gives b = —3/4 and yields a solution with the 
asymptotic form 

oc .-/-e~^''''\ (7.18) 

where r satisfies 

ar^ - -f 3 = . (7.19) 
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The solutions are = 9/2,9 in the case of the Ising model and = 5 ± in the 
Yang-Lee edge case. 

In terms of the expansion parameter k = we have ct(k) oc exp(— |(r/K)), 

and the large order behavior according to (|7.6|) is 

Remarkably enough, the large order behavior takes the same form, T{2k — as in ( [7.16| ). 

From ( |7.19D , we see that is complex when a > 1/12, and from ( |7.2U| ) we see that 
this is related to the non-unitarity of the Yang-Lee edge case, since it leads to asymptotic 
coefficients Uk that are not positive definite. This is also related to the Borel summability 
in the Yang-Lee edge case, since from the discussion following (|7.8| ) we see that the poles 
nearest the origin are a finite distance off the real axis in the Borel plane. There thus 



exists a real and physically acceptable Borel sum, presumably equal to the solution of |7l 
In the (unitary) Ising case, on the other hand (with real), the terms of the series are 
all positive, there is a singularity on the real axis in the Borel plane, and the series is not 
Borel summable. 

We shall shortly generalize the conclusion concerning Borel summability in the [l = 3) 
Yang-Lee case to all the / odd one- matrix models. We can also show that it is not affected 
by higher order exponential corrections. If we look for such corrections by expanding in e, 
we find that the coefficient r in the exponential is replaced by n+r + n_r*, but because 
the equations have a well-defined parity in e, n-\- + ri- must necessarily be odd.^il AH 
such allowed terms will not occur (some of them can correspond to singularities of the 
Borel transform B in other sheets of the complex plane), but we see in any event that 
no singularity will appear on the positive real axis and an integral like ( [7.12|) may be 
performed to define a unique real function. This function solves the differential equation 
of interest and is therefore the natural candidate for the partition function of the original 
matrix problem. 



"^^ The solution ei of the Unearized equation is a Hnear combination of the two exponentials 
involving r and r*, each multiplied by power series. The leading order corrections to uq involve e\ 
and thus correspond to exponentials involving 2r, 2r*, and r + r*. The argument then proceeds 



by iteration, just as remarked at the end of subsec. 7^ concerning the analogous properties of 
solutions to eqs. (|7.13a, 6| ). 
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7.^. The tricritical Ising model 

In the case of the tricritical Ising model, we perform the analogous computation by 
substituting u ^ u + tm v ^ v + ue^ into eqs. ( [4.33a-c| ) with T = h = 0. We introduce the 
ansatz e^/e^ = rw^/^ and e!^/e„ = rw^/-^, and remember that at leading order v ~ —v?l1. 
We then find a system of two linear equations for and e„, 

(-10 + - \r^) tu + (5r2 - lO) e„ = , (7.21a) 
(40 - 30r2 + - \r^) tu + (20 - lOr^ + r^) e„ = Q . (7.216) 

Writing that the determinant of the 2x2 matrix vanishes yields an equation for r^. 
Actually two roots are obvious because in eq. (|7.216|) the polynomial 20 — lOr^ + can 



be factorized. The equation for the other roots is then simply — 5r^ + 5 = 0, and the 
four solutions for are 

= 5± VS, i(5± VS) . 

Since these give real values of r, the theory is not Borel summable. This is again as 
expected for a unitary theory, in which the coefficients in the asymptotic expansion have 
fixed sign. The large order behavior is, up to the value of r^, the same as found in the 
earlier cases, eqns. ( [7.16| , |7.20|) , 



(|:)''r(2fc-i) . (7.22) 



2k 

Uk oc 



7.5. The general 'problem: Preliminary remarks 

In this subsection we shall explain the structure of the large order behavior in the 
general case. In the next section we shall then discuss explicitly the Borel summability of 
a large class of (p, q) models. 



The one-matrix models. We consider the string equation ( 3.13| ), Ri[u\ ~ z. To examine 



the leading large order behavior of perturbation theory, it is only necessary to know the 
terms in Ri [u\ of the form 



1 

Ri[u] = -Auu' + Y. n^~^n^^'~^^^ + ■ • • (7.23) 



(i.e. that contain at most one derivative of u factor. The next leading contribution is 
given by terms such as i.e. with a single factor of u' as well). From the 
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recursion relation (|3.14|) , we have for example that the coefficient of the highest derivative 
term u^^'"^) in ( [7.23|) is given by An = — and the coefficient of u'- is Au/l = (— 1)'(2/ — 
l)!!/(2^+i/!). 

Denoting the discontinuity of u{z) by e{z) and substituting u = uq + e as before (now 
uq ~ x^/'), to leading order for z large we find that e satisfies 

I 

= J2^ij<~'^^^'~^'^ ■ (7-24) 



Substituting the WKB ansatz 



^ = raj/' (7.25) 



then gives 

I 

= J2Aijr-''~^' , (7-26) 

an (Z — 1)^* order equation for with real coefficients. We shall determine the coefficients 
Aij and discuss the solutions of the equation in the next subsection. 

The subleading terms in Ri[u] mentioned after (|7.23|) are immediately deduced from 
the leading terms by noting that since Ri[u] is derived from an action (see eq. ( |3.15|) ), the 
operator acting on e is hermitian. Therefore the operator u^~^cP''~'^^ should be replaced 
by the symmetrized form ijM-'"^, d^^~^-^}, correcting ( 7.24 ) to 

= J2Aij{u'-'e'^''-'^^ + 1(2/ - 2j){j - 1)^^-2 w'e(2'-2i-i)) . (7.27) 
i=i 

To characterize more precisely the large order behavior, to next order we set 

= ru^l-^ + h- , (7.28) 
e u 

from which it follows, to the same order, that 

i_ = rV/2+r'=-iw('=-3)/2^,'A;(6+^(A;-l)j . (7.29) 
Substituting into (|7.27| ), we find 

= E + "'(^ - b' - 1 + 2 (6 + i(2/ - 2j - 1))] ) . (7.30) 
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We see that r remains a solution to ( [7 .261 ) and 6 = (3 — 2/)/4, independent of r. Dividing 
e by the leading term u oc z^/^ results in 

21 ^J2l+l)/2l 



OC z-(2^+i)/4^e , (7.31) 

generalizing ( |7.14| ) and (|7.18| ). In terms of the expansion parameter k = ^-(2^+i)/2i^ -^g 
find 



/ \ 1/2 21 



(7.32) 



We see that the r{2k - i) factor in ( fTT^ , [T?^ , [T]^) is general, owing in the case of one- 
matrix models to the special form of the subleading term in the equation ( [7.27|) satisfied 
by e (due to the fact that the original equations descended from an action principle). 

General {p,q) model. In the case of the general {p,q) model (eqs. ( ^4.16| ) and ( [4.24] )) 
there results a system of coupled linear differential equations for the variations e^, u^^evi (x) 
associated with the functions w(x), Vi{x) (the power Si of u is determined by the grading). 
As in the case of the tricritical Ising model considered in subsec. \7.4i at leading order we 
set e'^/eu = ru^/"^ = e^./e„. . We obtain, taking into account the leading relations between 
u and the Vi , a linear system for and e^. . Expressing again that the determinant of the 
linear system vanishes, provides an equation for the coefficient r. 

To determine more precisely the behavior of we have to consider subleading terms. 
As in the one-matrix case they can be determined by a hermiticity argument. Since the 
equations for u,Vi derive from an action ( |4.24|) , the linear equations for define a 

hermitian operator. Eliminating for example all the e„. yields an equation to next leading 
order for which can be expressed as a hermitian operator acting on (as was the case 
leading to ( [7.27| )). The coefficient of the subleading term that led to the T{2k — ^) behavior 
found in (|7.32| ) for the one-matrix models, since it only depended on the hermiticity of the 
operator acting on e, is thus universal for all the [p, q) models.!^ 



The (2fc)! large order behavior is also the generic behavior for D = 1 models coupled to 
gravity. More precisely 1 72], we found the contribution from genus k surfaces to go as f2k ^ r{2k — 
1). This was based on an instanton analysis that allowed an understanding of this behavior as a 
result of barrier penetration effects, typically of the form kl/A'', where the A is an instanton action 
given by the integral J dx ^yV(x) — E between the turning points. The perturbative expansion by 
itself for D = 1 does not fully determine the partition function and instead misses some essential 
non-perturbative feature of the problem (see also [pT|). 
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7.6. Borel summability 

Having shown that perturbation theory at large order has for aU (p, q) models the 
generic A''r{2k - i) behavior, we shall discuss for a class of models the existence of real 
positive values A, relevant for the Borel summability of the theory. For this purpose we 
need all equations only at leading order. We will show here that the relevant equations can 
be directly derived from the canonical commutation relations [P, Q] = 1 in the semiclassical 
limit The important remark which simplifies the analysis is that no derivative of u 
contributes to the equation for e at leading order (e.g. ( |7.24[ )). Therefore u appears only 
as a scale parameter and can be eliminated from the equations. 

We have seen that in the semiclassical limit, the operators P, Q take the formil 

P(d, u) = uP/^P (du-^/^^ , Q(d, u) = u'^/^Q (dw-^/2 j . 

Prom now on we call P{z), Q{z) the two polynomials P{z = dw~^/^, 1), Q{z = dw~^/^, 1). 
We recall that in subsec. \4.3i we have determined P, Q for all models such that p = (2m + 
1)9 ±1- 

One-matrix models. Before discussing the general p = (2m + 1)^ ± 1 models, let 
us return to the one-matrix models. From the analysis of the corresponding non-linear 
differential equations, we have learned that the variation e of the specific heat u{x) has for 
X large the asymptotic form 

e'/e ~ rv^ , (7.33) 

where r is a constant determined by an algebraic equation. Since the function u can be 
treated at leading order as a constant, we can rescale d, i.e. set u to 1. Equation (|7.33| ) 
can then be written as a commutation relation 



de = e(d-fr) ^ /(d)e = e/(d + r) . (7.34) 
Then the operators P, Q are simply 

Q = d2-2, P = P2z+i(d) = (d2-2)'+'/' . 



Note that in this subsection the normalization of u{x) corresponds to generic potentials, i.e. 
the minimum residues of double poles is 1. 
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The equation for e is obtained by expanding at first order in e the commutation relation 
[P,Q] = 1. Setting 

k=0 

({, } means anticommutator and this form takes into account the antihermiticity of P), we 
obtain 

[{e,i?(d)},d2-2] + [P,-2e] = 0. 
Using the commutation relation (|7.34|) to commute e to the left, we find the equation 



-(2rd + r2) (i?(d) + i?(d + r)) -2(P2Z+i(d + r)-P2;+i(d)) = . 

The first term vanishes for d = — r/2, so this must give as well a zero of the second 
term. Taking into account the parity of P21+1, we obtain 

P2Z+i(r/2) = ^ {r^-8f+'/^ = 0. (7.35) 

The polynomial -R(d) is then determined by division. The function (z'^ — if^^^"^ is also 
proportional to C'^^j_i(^) where C2;_|_i is a Gegenbauer polynomial defined by analytic 
continuation in 0. Note that the number of zeros is exactly the same as the number 
of operators in a (p = 2/ — 1, 2) minimal conformal model [Q. This is a property we shall 
meet again in the general case. In the one-matrix case there is a natural explanation: 
the steepest descent analysis shows that the number of different instantons (see the next 
subsection) is related to the degree of the minimal potential corresponding to a critical 
point. This degree in turn is also related to the number of relevant perturbations.!! 
The l.h.s. of the equation has a useful integral representation: 



For I even, eq. (|7.35| ) is an odd-order equation that will have at least one real solution 
for r^, positive as is obvious from the integral representation (|7.36D. The series therefore 



cannot be Borel summable. 

For I odd, on the other hand, the equation ( [7.35| ) for has no real solutions and 



therefore we expect the solution of the differential equation to be determined by the per- 
turbative expansion. 
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We thank F. David for this remark. 
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General (p, q) problem. In the general (p, q) case, in classical limit as above and after 
the same rescaling, we have: 



where P, Q are polynomials of degrees p, q, respectively, and -R, S are polynomials with 
the same parity as P, Q but with degrees p — 2, q — 2. 
The equation for e then leads to 



The polynomial P(d + r) — P{d) has degree p — 1 in d, while R only has degree p — 2. 
An equivalent property is true for Q, S. Thus the polynomials P(d + r) — P(d) and 
Q{d + r) — Q{d) must have at least one common root. Note that the first polynomial has 
p — 1 roots and the second q — 1- Moreover these roots are symmetric in the exchange 
d I— > — r — d. The existence of a common root thus leads to (p — l){q — 1) values of r, up 
to the symmetry. Note that the number of zeros is again exactly the same as the number 
of relevant operators in a (p, q) minimal conformal model Q of gravitationally dressed 
weights Am,n (see eq. ( |S.26D ) with l<?7,<q' — l,l<m<p — 1, and possessing the 
symmetry Am,n = Aq_n,p-m- The explanation of this relation is probably again that 
the number of different instanton actions is related to the degree of the minimal potential 
needed to generate a critical point in the multi-matrix model, and thus to the number 
of different relevant operators. Also we note that we are studying a general deformation 
of a critical solution and therefore the appearance in some form of the relevant operators 
should be expected. 

The above condition determines the possible values of r when the polynomials P and 
Q, i.e. the differential operators, are known in the classical limit. Examples are provided 
by the models p = (2m + l)q±l where these polynomials have been determined explicitly. 
The simplest examples are provided by the {q + 1, q) models, i.e. the unitary models which 
we examine below. 

Finally we verify that we can indeed find the polynomials R, S. We call a the common 
root and assume first that a ^ —r 12. Then the parity properties imply that — r — a is also 
a common root. Setting 



g = Q(d), P = P(d) = g^/'^(d), 5Q = {^(d),e}, 5P = {P(d),e}, 




(P(d + r) - P(d)) [S{d) + Sid + r)) 

- (Q(d + r) - g(d)) (P(d) + P(d + r)) = . 



(P(d + r) - P(d)) = (d - a)(d + r + a) P(d) 
(Q(d + r) - Qid)) = (d - a)(d + r + a) Q{d) , 
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we find that R and S are solutions of 



i?(d) + R{d + r) = (d + r/2)P(d) , S{d) + S{d + r) = (d + r/2)Q(d) 

Note tliat tliese equations satisfy both the degree and parity requirements. 
If a = —r/2, the situation is even simpler 

R{d) + R{d + r) = (P(d + r) - P(d))/(d + r/2) 
S{d) + S{d + r) = (Q(d + r) - Q(d))/(d + r/2) . 



Application: The unitary models. In subsec. we have seen that the differential 



operators P, Q may be written in the classical limit as 

P = 2Tp{d/2) , Q = 2T,(d/2) , 
where Tp is the p-th Chebychev's polynomial: 

Tp (cos v?) = cos p^p . 

As explained above, taking into account the degrees of the polynomials R and 5", we 
conclude that the polynomials Tq{{r + d)/2) — Tg(d/2) and Tp((r + d)/2) — Tp(d/2) must 
have a common root a = 2 cos lpq. Let also set a + r = 2 cos i/jq. We have 

cospVo = cosp^jQ and cosq^'^o = cosg^^o • 

The solution is 

2mTC 2mr 

V'o = ±^0 H = =FVo H • 

p q 

Since r = 2cos'i/'o — 2cosv9o5 excluding the solutions r = which is not acceptable, we 
have the different solutions: 

r = ±4 sinm7r/psinn7r/g, < 2m < p , < 2n < q . 

It is easy to verify that these results agree with the explicit solutions of the (2,3), 
(4, 3) and (4, 5) models. The results also show that, as expected, all unitary models lead 
to non-Borel summable topological expansions because all terms of the series have the 
same sign. These models thus suffer from the same disease as the pure gravity model. 
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Note finally that the number of different values of r is indeed the same as the number of 
operators in the minimal (p, q) conformal model. 

The general p = (2m + 1)^ ± 1 models. For m 7^ 0, r is a solution to more complicated 
algebraic equations. In the notation of previous subsection, we still have 

2nn 

tpo = ±ipo H . 

Q 

We set 

SO that, making a choice of signs, we have 

TITT 

j3 = — , r = 4 sin a sin /3 = 4 sin a sm( nn / q) , 
Q 

where 

^ (^^/^^ sm{{p - 2ql)a) sm{{p - 2ql)0) = . 

We note that sin((p — 2ql)(3) = sm{mvp/q), which can be factorized. We thus find an 
equation for a, 

= sin((p-2Q0«) = 0. 

The function A{a) satisfies the differential equation 

pA{a){cosqay — qA'{a) cos qa = K{j>, q) cos a , 

where X is a constant. This equation implies that A{'K/2q) A{?>tx /2q) < and thus A{a) 
vanishes at least once in the interval (0, tt). We conclude that for all these models the 
topological expansion is not Borel summable. 

7. 7. Instantons and large order behavior for one-matrix models 



Following we now show that the large order behavior analysis coincides with 
an instanton calculation of barrier penetration effects for the one-matrix models. For g 
negative and finite, instanton effects, corresponding to a single eigenvalue of the matrix 
model climbing to the top of the barrier, are responsible for the divergence of perturbation 
theory at large orders. In the large limit this effect is suppressed as e~^^^^^ and thus 
the tree-level free energy is analytic for > g > gc- Near gc, however, we shall show 
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that the function K{g) vanishes as {gc — so that in the double scaling limit 

the tunneling amplitude remains finite and moreover given identically by (|7.31| ). In what 



follows we use the notation of subsec. 2.1. 



In the large limit, the free energy is given by 

F = \nZ = (^J dXdfip{X)p{fi)ln\X- j dXp{X)V{X)^ , (7.37) 

where p(A) is the eigenvalue density. We wish to calculate the variation of the action in 
the large limit when one eigenvalue is displaced from position to A/. This variation 
of the distribution p(A) is of order 1/N and therefore can be obtained by the first order 
variation of ( [7.37| ), 

5F^N^ j dX dp{X) (2 j dii p{p) In |A - ;u| - ^^(A)^ (7.38) 

(which is non- vanishing since A is outside the support of p). For 

5p{X) = N-\5{X- Xf) - 5{X- Xi)) , 

eq. ( [7. 381 ) can be written (using eq. ( |2.13| )) 

5F = -2N j ' dX (uj{X) + — F'(A)) , 

= / dAv/(^'(A)) +i?(A) = -2iV / dAu;sing(A), 

9 Jxi JXi 



(7.39) 



where casing (A) is the singular part (which scales in the continuum limit) of a; (A) 
N-^{tr (M- A)-^). 

The case 1 = 2. In the special case / = 2, the explicit form ( |2.26| ) of a;sing(A) gives 



bF = -— \ dA(A2 + l + ia2)v/A2-a2 
9 JXi 



(7.40) 



with 0? = |(— 1 + v^l + 12(7) . It is easily verified (and intuitively clear) that the vari- 
ation of the action is minimized by moving the eigenvalue A^ = ±a at the edge of the 
distribution. The final position Aj is taken as the top of the barrier, A^ = —1 — 0? /2, at 
which [d/dXf)5F = 0. For g small and negative, it is easy to verify that dF is strictly 
negative and therefore tunneling is suppressed in the large limit. Let us now examine 
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what happens as g Qc = —1/12. Then A/ A^, and substituting = |(— 1 + y/x) into 
( [7.4q) , and evaluating in the hmit of small x = 1 — g/g^, gives 

5F = -^VQNx^/^ , 

in agreement with the result ( |7.14| ) (with z = N'^^^x). The conclusion is simple: instanton 
effects are suppressed in the large limit for fixed g > gc but remain finite in the double 
scaling limit. The scaling limit of the / = 2 one-matrix model thus corresponds to a complex 
solution of the Painleve I equation, presumably the complex Borel sum of the perturbative 
expansion. 

The case of general I. Returning now to the expression (|7.39|) for the case of general 
I, we take A^ to be the cut end-point Xi = a and A/ again the top of the barrier. In the 
scaling limit, we know that from the steepest descent analysis of subsec. \2.3i that co'sing(A) 
has the form (eq. (|2.36|) ) 



•-'A ^ 



where x = 1—g/gcis considered as a function of the cut end-point position a (as in (|2.36|) ). 
dF then becomes 

5F = -2Nh-^'^ I ^ I x{s) {s - A)-^/2 , (7.42) 
where A/ is determined by the condition 

= ^SF = -2NLj,,ng{Xf) . 
Changing variables in ( |7.41| ) to s = A + t{a — A), this condition can be rewritten 

l'^^{Xf + t{a-Xf)y-'=0. 

Comparing with eq. ( [7.36| ), we find the solution A/ = a(l — r^/8). Substituting in ( [7.42| ), 
we derive ^ 

■^6F = -2iV6-i/2 f ' dA (a - A)-^/^ = Nr(2a/b)^/^ . 

From ( |2.37| ) we finally conclude that 

^ = r^='l ^ 5F~lne(a:)~--^rz(2m)/2Z ^ (743) 
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yielding a result in agreement with the result ( 7.31 ) of the large order behavior analysis. 



Remarks. The interpretation of the large order behavior calculation in terms of in- 
stantons in the steepest descent method allows us to understand the direct correspondence 
between the property of Borel summability and the existence of the original integral. It 



has been noted |]71| that according to whether Hs odd or even, the original minimal matrix 
integral is well-defined or not because the integrand goes to zero in the first case while in 
the latter case it blows up for M large (see subsections \2. Sj , \2.4 }- The direct calculation 



given here of the instanton action fT^, using steepest descent, confirms that when the po- 



tential is unbounded from below the instanton action is indeed real and the series therefore 
non-Borel summable. Note moreover that replacing the minimum potential by a potential 
of higher degree which would be bounded from below does not solve the problem. Indeed 
the instanton result depends only on universal properties and thus instantons would still 
appear. In the latter case they would reflect the existence of another minimum of the 
potential, lower than the one in which one assumes the eigenvalues are contained. Such a 
minimum would also invalidate the direct calculation. 

7.8. I = m perturbed by I = m — 1 

The calculation of the previous subsection can be adapted to relate the loss of Borel 
summability when flowing from a model with / odd to another model with / even. For 
definiteness we consider here the case of flowing from / = m to / = m — 1, with m odd. 
The general case is treated identically. 

To describe a perturbation of an / = m model in the direction of an / = m — 1 model 
we write, in the notation of sec. |^, 

= -(1 - 6)"^ - ^(1 - 6)"^-! (7.44) 

(with W as in (|2.49|) ). A particular case would be the I = 3 Yang-Lee edge singularity 
perturbed by / = 2 pure gravity, and the reader might derive further intuition by substi- 
tuting this particular case in the more general formulae that follow. (The matrix model 
potential V{M) associated to ( [7. 441 ) is easily reconstructed by recalling that a potential 
V{M) = 'ZpOpM^^ in general leads to W(b) = 2J2p (§^TJ^ 9pb^^ but the exphcit form 
of V but will not be necessary in what follows here.) 
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The double scaling limit in this case involves taking 6 = 1 + N~'^/^-^'^~^^^u, and ^ = 
j\^-2/{2m+i)rp^ The resulting (all-genus) string equation that describes the flow of the 
/ = m model to the / = m — 1 model is 

Rm[u]+TR^_i[u]= z . (7.45) 

The scaling hmit corresponds to z ^ {g - ^jA^2m/(2m+i)^ T ~ m ~ z^/"^ ~ ^2/{2m+i)_ 
Thus, at leading order (genus zero), z and T are large with z/T"^ fixed in such a way that 
the leading order equation becomes 

u^ + Tu'^-^ = z. (7.46) 

To treat this (mixed) case more easily, we change variables to z = z/T"^, u{z) =u/T, and 
= 7"-(2"^+l)_ Eq_ ( [7.45] ) reads with these notations 



Rm[u]+Rm-l[u]^i (OH)' 



where the notation R indicates the substitution d — > Kd, and perturbation theory corre- 
sponds to an expansion in small k. Eq. ( [7.46|) becomes 



By the same linearization procedure (tt 1-^ w + e) as used before, we find at leading 
order (in the notation of ( [7.24] )) 



m — l 



Substituting the WKB ansatz e'/e = K~^u^/'^r gives an algebraic equation for r, 

m ^ m — l 

= ^ A^,, r^"^-^^- 4 ^ r^— 2-2. . (7.47) 



.=1 



depends on 2; = z/T^ implicitly through u = u/T. In the limits u large and small (and 
hence via (|7.46[ ) ' z large and small) , we recover the cases I = m and I = m — l respectively. 



Since we are assuming m odd, for some large enough value of z (and hence of w = u/T 
classical), the solution of eq. ( |7.47| ) goes from complex to real. Thus for T large enough, 
the equation for has real roots and the Borel summability is lost. 
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Now consider the value of ^ for which the interpolating potential V{X) corresponding 
to ( [7 .441 ) begins to admit real instantons. Employing the same steps leading from this 
potential to the result ( [7.43| ) for the instanton action, we find that the instanton action 
becomes real at precisely the value of T for which the Borel summability is lost, since both 



are determined by when the same equation (|7.47| ) has real solutions r. Thus above some 



value of ^, which corresponds to the critical value of = T/u in ( [7.47| ) for which the 
roots in become real, we see that the interpolating potential substituted in ( |7.37D allows 
real instantons, giving a physical interpretation for the loss of Borel summability.0 

7.9. Other properties of the equations. Moveable singularities 

We have investigated earlier in this section the large order behavior of perturbation 
theory and concluded that in many cases of interest, in particular pure gravity, the series 
is not Borel summable. In this subsection, we investigate a few other properties of these 
equations, again first considering the Painleve equation. 

The Painleve I equation. Eq. (|7.9|), 



u^{z)-^u"{z)=z, (7.48) 

is a second order differential equation. The general solution depends on two parameters 
fixed by the boundary conditions. It has as moveable singularities double poles of residue 
2, as one can immediately verify, and it can be proven that these are the only moveable 
singularities [^. Moreover it is easy to show that any real solution of the equation has 
an infinite number of double poles on the negative real axis. This leads to logarithmic 
singularities of the free energy with weight 2 and thus to double zeros of the partition 
function. (If we consider a general potential (not even), the free energy is divided by 2 
and thus the partition function has simple zeros.) 

A more quantitative analysis can be obtained by transforming the Painleve equation. 
We set 

u{z) = z^^^v{y) with y = \z^l^ , 



"^^ For D = 1 models coupled to gravity, "multicritical" models, determined by a potential V^(A) 
that has its first s — 1 derivatives vanishing at the critical point, were considered in [72|. These 
models have string susceptibility 7 = — (s — 2)/(s + 2), and their physical interpretation is in general 
unclear (but see ||Tl[| ). They do however provide a useful arena for studying how perturbations 
from one model to another can destroy Borel summability. We expect again a mod 2 grading that 
determines which models are Borel summable, and as well which models may flow to one another. 
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and then approximate the equation by keeping only the leading order terms for z large. 
The approximation is also valid when f ' ^ f , e.g. near the double poles. It is finally very 
easy to write the resulting equation: we simply replace u i— f , 2; i— ^ 1, in ( [7.48| ) to find 

v\y) - \v'\y) = 1 . 

The equation can be integrated once 

- ^^'^(2/) = "^{y) + const . 

The solution has in general the form of a Weierstrass elliptic function V{y/\/2), with a 
doubly periodic lattice of double poles in the complex plane. 

Boundary conditions. To understand the role of boundary conditions, we calculate 
the variation of a solution uq of eq. (|T]^) when the boundary conditions are infinitesimally 
changed, 

u{z) = uo{z) +e(z) . 
e satisfies the linearized equation ( [7.13a| ) 

e" — 6uoe = . 

It follows that e is a linear combination of two independent solutions e± which can be 
obtained, for z large, by a WKB ansatz (see eq. ( [7.14| )), 

e± oc ^-i/8e±(4V6/5)-^/\ 

z — >+oo 

This form shows that the equation is unstable for z large both when solved with z increasing 
or decreasing since a small change in the boundary conditions produces an exponentially 
increasing change in the solution. Moreover solutions with the asymptotic expansion ( |2.69| ) 
form only a one-parameter subset of solutions, because the coefficient of the growing expo- 
nential must vanish. The meaning of this remaining parameter is, at this point, obscure. 



This result is related to the property that the topological expansion (|2.69|) is non-Borel 
summable (see subsec. |7.^ ), and thus that its sum is ambiguous. 
From these properties the following conclusions can be drawn: 

(i) For the solutions with the correct asymptotic behavior, the existence of double poles 
implies that the operator H = — (d/dz)^ + f{z) has a discrete spectrum. This property 
is surprising since the operator if is a representation of the multiplication operator A (see 
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subsec. 3A). It seems to indicate that, in contrast to what is seen in perturbation theory, 
the hmiting distribution of eigenvalues of hermitian matrices is discrete. 

(ii) A careful analysis of the loop equations in the scaling limit [^, however, has 
shown that the spectrum of H cannot be discrete (see subsec. \2.8i ). 

(iii) Real zeros of the partition function, which is a sum of positive terms, are not 
expected. 

(iv) The initial matrix integral can only be defined by analytic continuation since 



the critical value Qc is negative. It can be argued that, as a consequence of this 



continuation, the function f{z) is actually a complex solution of Painleve equation, most 
likely the complex Borel sum of the series, which for z —>■ +00 behaves like ^/z and like 
for z —00. 

The conclusion is that the matrix integral does not define the sum over topologies 
beyond perturbation theory. 

The I = 3 model. In the case of the I = 3 model, the situation is more favorable. 



Since a real solution exists [^, the perturbation series is presumably Borel-summable. 
The initial matrix integral is convergent in the scaling region. A solution has been found 
numerically that behaves like z^^^ for z 00, —z^^^ for z — > —00 (for finite, the matrix 
model is defined both for g < Qc and g > Qc and therefore the solution should have an 
asymptotic expansion both for z +00 and z —00), and which has no real double 
poles. The spectrum of H is thus continuous. The solution thus satisfies all requirements 
and is presumably the correct solution to the initial problem. The non-positivity of the 
terms in perturbation theory, however, is itself a pathology of the model, and indeed it has 
been found [|T9[ that the / = 3 model corresponds to the Yang-Lee edge singularity, i.e. to 



the critical point of an Ising model in an imaginary magnetic field. 

Note that as in the preceding case, some information about the structure of the solu- 
tion in the / model, in the complex plane for z large, can be obtained by setting 

u{z)=z^/'v{y) , y = ^i+i/27(l + l/20. 

Again the effect of the transformation is to yield an equation for v which at leading order is 
the same as for u except that z in the r.h.s. has been replaced by 1. Because the equation 
follows from an action principle a first integration can be performed. 

General models. Although there is no general proof, there is very good evidence for a 
systematic difference between the odd and even cases in the one-matrix problem: Before 



115 



taking the large continuum limit, the operators A and B have a continuous spectrum 
and form a regular representation of the canonical commutation relations. This property 
still holds for the operator d^ — u when / is odd, while for / even the operator has a discrete 
spectrum.l?^ 



Moveable singularities. All the differential equations of the (p, q) models (subsections 
\3.3i , \4-2^ have as moveable singularities double poles. For all the equations investigated 
so far, the residues of these double poles are even integers. Although it was shown in 



that solutions to these equations with poles on the real axis are not physically relevant, 
for completeness we briefly recall here some previous results and add results for the Ising 
and tricritical Ising models. 

For the one- matrix models, we set 

u{x) ~ a/x^ , Ri[u] ~ pi{a)/x^^ , 

and using the recursion relation between Ris gives 

21 + 1 

Pi+i{a) = ^^y^-^ (a -/(/ + !) ■ 

The solutions of the equation Ri [u] = x have double poles with residues belonging to the 
set 2, . . . , j(j + 1), . . . , — 1) (all of which are even integers). 

For the critical and tricritical Ising models, a similar analysis shows that the possible 
residues are 2, 10 and 2, 6, 8, 14, 30, respectively — again even integers. In general for a 
[p, q) model, it can be shown that the residues are integers whose smallest and largest 
values are respectively 2 and -^{p'^ — i){q'^ — !)• These results are calculated using a 
normalization for the differential equation generated by an even matrix model potential, 
and the poles above all correspond to double zeroes of the partition function. For generic 
models with non-even potentials, on the other hand, the residues are all divided by two 
and the partition function has only simple zeroes. 
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In [74 1, it was further shown that it is impossible to flow by means of the perturbation of 



eq. ( |7.44 ) from the real solution of the / = 3 model found in [71| to a solution of the / = 2 pure 



gravity model. This result generalizes to show that flows from real solutions of arbitrary /-odd 
models to solutions of /-even are impossible |]75[| (see also |]28|] ), further distinguishing the /-odd 
and /-even cases. 
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8. Matrix canonical commutation relations, discrete action principle and dis- 
crete KdV flows 

As we shall show, many of the algebraic properties of the differential equations of the 
continuum limit considered in previous sections here are consequences of the properties 
of the recursion formulae when is finite. We shall examine the one-matrix case first 
because it is simpler and more explicit results are obtainable. In what follows an essential 
role is played by the matrix equivalent of differential operators, what we shall call local 
matrices: matrices X that have non-vanishing matrix elements only in a strip of finite 
width, i.e. satisfying X^n = for |m — n| > Ix- (For other consideration of the discrete 
case, and in particular the connection to Toda theory, see e.g. []76[| .) 

8.1. The one-matrix case 

We have seen that all perturbative properties of the models can be obtained from the 
commutation relations [P, Q] = 1 of the continuum differential operators which in turn are 
direct consequences of the matrix commutation relation ( |3.6| ) of two matrices representing 
the operations of differentiation and multiplication. 

We shall therefore assume that we are given two local matrices B and C, i.e. such that 
the matrix elements Bmn ? C'mn are no n- vanishing only for \m — n\ < 1 and | m — n | < ^ — 1 
respectively, which satisfy 

[B,C] = 1. (8.1) 

Moreover we assume B symmetric and C antisymmetric, and B^^n+i non-vanishing for all 
values of n. 

Reconstruction. With the hypothesis B^^n+i 7^ Vn, we can construct by induction 
an algebraic basis of polynomials Pn such that 

AP^(A) = BmnPn{X). Pq = coust , P_i = . (8.2) 

With the Pn defined, we can introduce a matrix A representing d/dA, 

p:^{X) = A^nPnW . (8.3) 

A is by definition a lower triangular matrix {A^n = for n > m) which also satisfies 



[B,A] = 1 
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(8.4) 



The matrix A, however, is not necessarily local. Setting X = ^{A — C), from 
and we see that 

[B,X] = 0. 

Because B is symmetric, both the symmetric and antisymmetric part of X commute with 
B. But an antisymmetric matrix commuting with B necessarily vanishes, as can be proven 
by a double induction on m and k on the matrix elements Xm,m+k (see app. D.l). Thus 
X is symmetric, and 

X = A + A^, C=l{A-A^). 

Because A is lower triangular with vanishing diagonal elements, it is uniquely defined 
by the second equation above and is thus local. Then X is also local, Xmn = for 
\m — n\ > / — 1. Again the commutation relation [S,^] = translates into recursion 
relations for the matrix elements. It is shown in app. D.2 that the commutation relation 
has a general solution of the form X = V'{B) where V'{B) is an arbitrary polynomial of 
degree I — 1 in B. We thus find 

A + A'^^V'iB), (8.5) 

a formula analogous to (p^). Note that C is then a representation of the operator d/dA 
on the functions e-^(^)/2p„(A). 

Equation ( p.5| ) however does not embody the full content of eq. ( p.l| ). Let us indeed 
assume the existence of a matrix A with Amn only for < m — n < I — 1, and satisfying 
(p75|). We consider the commutator 

J=[B,A] . 

Because A is lower triangular (Amn = for n > m), J satisfies Jmn = for n > m. 



Moreover using (|8.5|) we have 

J^={[B^A])^ = [A^^B] = J. 
Demanding [B, A] = J = 1 then leads to an additional equation 

Jmm -^m,m+l-4m+l,m -^m — l,7Ti-^m,m — 1 1 ; (^■^) 

with solution 

ITT' = Bm—l,m Am,Tn—l ■ (§•''') 
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The property that eq. ( |8.6| ) is the difference of the two successive equations (|8.7|) is un- 
doubtedly related to the property that in the continuum limit the relation [P, Q] = 1 has 
yielded the derivative of the string equation. 

Of course eq. (^^) can be directly obtained by comparing the terms of highest degree 
in eqs. Q and (^). Indeed setting P^(A) = p^A" + 0(A"-i), we obtain 

Pm Pm,m+lPm+l ; ^^Pm ^m,m — lPm — l : 



from which follows ( |8.7| ) . It is not proved however from this argument that this is the only 
missing equation. 

We have thus recovered the recursion relations ( |2.79| , |2.80| ) of the one-matrix case 



starting only from the commutation relation and some conditions on the matrices C, B. 
The recursion relations determine the matrix B, however, only as a function of its first 
/ — 2 matrix elements, whereas it is completely determined in the matrix model or when 
it is related to orthogonal polynomials. It follows that even when the integral J dAe~^*^^-* 
exists, the polynomials Pn(A) are not in general the orthogonal polynomials corresponding 
to the measure dAe"^'^'*'^. In particular the symmetric matrix 5', 



Smn = J dAe-^(^)p^(A)P„(A) 



is not in general diagonal. 



Remarks, (i) It is easy to verify that S commutes with B and A. Therefore, as shown 
in app. D.2, if S is local, it is a polynomial in P: S = S{B). Then expressing that it 
commutes with A we find S'{B) = 0, i.e. the polynomial S{B) is a constant and 5" is a 
multiple of the identity. Therefore, when the first / — 2 elements of B are not those given 
by the matrix model, the matrix S is not local. 

(ii) Based on considerations concerning the instability of the sequences generated by 
the recursion relations (analogous to the instability of the continuum differential equa- 
tions), we believe that when the / — 2 coefficients of B do not correspond to orthogonal 
polynomials, the continuum limit does not exist because the matrix elements Pn,^? B^^n+i 
do not have a smooth behavior for n large. In particular when the integral J dAe"^*-^-* 
does not converge, no real set of matrix elements has a smooth behavior and the continuum 
limit can only be defined when the P„ are taken orthogonal with respect e~^, with the 
integral calculated along a complex path on which it converges. This argument applies to 
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the even critical models. It emphasizes the importance of the boundary conditions for the 
differential equations representing the string equations in the continuum limit. 

(iii) We recognize in the decomposition of V'{B) into the sum of a lower and upper 
triangular matrix a structure very similar to the one encountered in the continuum case 
(sec. ^) in which we introduced 



An action principle. As in the continuum, the discrete recursion relations {^1^, ^^8^ ) 
can be obtained from a variational principle by varying the discrete action 

S{B) = -J2^^^B^-i,n + trV{B) . (8.8) 

n 

8.2. Discrete form of KdV flows, one-matrix case 

The solutions of the commutation relations depend on the coefficients of the poly- 
nomial V'{B). Let us explore what happens when these coefficients vary. We denote a 
parameter characterizing V by t. Since the polynomials Pn also depend on t, we set 

dP 



dt 



It follows that Timn = for n > m. Differentiating eqs. ( ^.2|) and (|8.3|) with respect to t 
and using the linear independence of the P^s, we find 



8 A 

— =[E,A], (8.10a) 

— =[E,S]. (8.106) 



dt 

For a set of parameters ti and corresponding matrices E^, (|8.9| ) implies 

a,E, -a,E,+ [E,,E,] = . 



Eqs. ( p. 10a, 5|) thus define a set of commuting discrete KdV flows analogous to the flows 



( ^.20| ) of the continuum differential equations. 
Since B is symmetric, ( |8.10fc| ) implies that 



[S, E + E'^]=0. (8.11) 



Differentiating ( |8.5| ) with respect to t and using (|8.10a| ), we find after transposing: 



[A, E + E-]=^, (8.12) 



dt 

120 



where dV/dt in the r.h.s. means derivative with respect to t at S fixed. From and 
[B, A] = 1 we see that a particular solution is 

E + = -^^^p^ + const . (8.13) 

The constant corresponds to a trivial rescaling of all polynomials Pn- Since E is lower 
triangular, it is entirely determined by (|8.13| ) — in particular we see that it is local. 

A general solution differs from ( p.l3| ) by a matrix which commutes both with B and 
A. A local matrix which commutes with B and A, however, must be a multiple of the 
identity (see remark (i) above). Therefore if E is local, eq. ( p.l3|) yields the most general 
solution of eqs. ( |8Tl| , |T2|). 

More general solutions are necessarily non-local. It can be shown that a general 
symmetric matrix X that commutes with both A and B vanishes only when its matrix 
elements Xmn with m + n < I — 3 vanish. The reason is easily understood: Variations of 
the matrix B can be separated into variations induced by a variation of the function V, 
corresponding to local matrices E of the form ( |8.13| ) , and variations induced by variations of 
the 1 — 2 first matrix elements of B not determined by the recursion relations, corresponding 
to non-local matrices E. This result indicates that KdV fiows generated by local and non- 
local matrices play a very different role. 

The result ( fj.l3D can be easily verified in the matrix model. Let ti be for example the 
coefficient of B^ in the potential V{B). Differentiating the orthogonality relations ( |3.1|) , 
we find 

+ Ef = ^ = B^ . (8.14) 

Note that the commutation relation [5, A] together with eq. ( ^.5|) constitute a particular 
case of this result: They correspond to an infinitesimal variation of V{M) resulting from 
a translation ofM: M M + tl, which induces an identical translation of B. 

Remarks 

(i) Since C is also a representation of a fixed operator, we expect a flow equation for 
it. One verifies indeed that ^(E — E^) (an antisymmetric matrix) generates a fiow for C 
and B. We denote by S!^ the antisymmetric matrix defined by 

Ei + EJ = B\ B\ = \{Ei - Ef ) . 
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As above we parametrize the potential V{B) as V{B) = Yli^i^^- With these definitions, 
we can write 

^ = [5;.^]. c = Y,iuB!,-\ 

i 

in complete analogy with ( [4.23|) . 



(ii) The fiow equations refiect the property that the matrices B and A are the repre- 
sentation in a basis which depends on some parameters of a fixed operator corresponding 
to a multiplication by A. Indeed two matrices B corresponding to two sets of parameters 
ti and t[ of the initial integrand are related by a linear transformation R, 

B{U) = R{U,t[)B{t[)R-^{U,t',) . 



Differentiating with respect to t^, one finds (|8.106|) with 



dR. 
Oti 

This definition of may differ from the previous one by the addition of a matrix com- 
muting with B. 

Conversely, let us assume ( |8.10^ ) and in addition that B is diagonalizable. In our prob- 
lem, this condition is satisfied when B is defined via polynomials orthogonal with respect 
to the measure dAe"^*-^-*, ensuring that the functions e~^'^^-'/^Pn(A) are the eigenvectors 
of B corresponding to the eigenvalue A. The other cases require a separate investigation 
of the behavior of -Bmn for m, n large. 

We can then express B in terms of its eigenvalues A^ (we write them as a discrete set 
although they could as well be continuous), 

B = R-^AR , 

which transforms ( p.l06 ) into 



sr = [^-^]- 

with 

E = RER-' + ^R-' . 

ot 

Since A is diagonal, the equation implies that E is also diagonal. It follows that 

122 



The spectrum of B is thus independent of t. This confirms indirectly that in the large N 
limit the spectrum cannot become discrete. 

The resolvent. It is easy to verify that, in direct analogy with the continuum case, the 
diagonal matrix elements of the resolvent G{z) = {B — z)~^ satisfy a four term recursion 
equation. Defining 

7n=[G(.)]^^, 

one finds (for an even potential) 

{in - 7n-l) = rn+lln+1 + Tn {in - In-l) - rn-lln-2 ■ 

(In the continuum limit, this equation becomes eq. (A. 7) of app. A.l.) 

8.3. Multi-matrix case 

The general case is not in the same satisfactory state. We shall therefore give a few 
general results valid for the multi-matrix case and then consider separately the two-matrix 
case for which more detailed results can be obtained. 

In the continuum limit, all the matrices reduce to two differential operators. It would 
be useful to derive, as in the one-matrix case, all properties entirely from consideration of 
two local matrices satisfying canonical commutation relations. Only partial results can be 
easily obtained. The major reason for this new difficulty is that while in the continuous 
case only two operators appear, the matrix model is originally defined in terms of several 
matrices forming equivalent representations of the commutation relation. As we shall see, 
the general problems we would like to solve are the following: (i) Characterize all pairs of 
local matrices B and A satisfying the commutation relations A] = 1. (ii) Understand 
whether all solutions correspond to a matrix model, or whether additional conditions are 
necessary to insure such a correspondence, (iii) Find the subclass of matrices B which lead 
to a continuum limit. Characterize all fiows generated by local matrices E: 

dB . . OA . 
dt L ' J ' L ' J ' 

such that B and A remain in the class defined above if they belong to this class for t = 0. 
Unfortunately we shall only give partial answers to these questions. 
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An action principle. It is, however, easy to verify that eqs. ( ^4.8| ) can also be directly 
derived from a more complicated variational principle than ( p.8|) , taking for action 

S {Ba) = J2^rVa (S^")) - 5^ c„ tr 

a=l a=l 

- ^n(ln[Si]n-i,n + ln[B^-i]n,n-l)- 
n 

Discrete canonical commutation relations 

Thus let B, A be two local matrices be such that Bmn 7^ only for 1 — r<n — m<l 
and A^n 7^ for m > n with A] = 1. (In the notation of sec. ^, we have in mind here 
the matrices Ai and Bi.) We again assume Bm,m+i 7^ for all m > 0. 

(i) It is proven in app. D.l that the matrix A is uniquely defined in terms of B. The 
matrix A can then be explicitly determined by first constructing inductively an algebraic 
basis of polynomials P„ such that 

\Pm{\) = B^r,P^{X) , (8.15) 

and thus 

P'^{\) = Ar^MX) . (8.16) 

(ii) The preceding results do not involve the locality of A and B. Expressing that A 
and B are local matrices yields r + 1 recursion relations for the matrix elements Bjn+k-i,m, 
< k < r. To prove this result, one solves [S, A] = 1 following the method of app. D.l. 
This leads to recursion relations for Ara+k,m which depend on the Ajn+k-n,m-, n = 1, ...,r 
which have already been determined. Expressing that Am+k,m vanishes for all k > s leads 
then to the corresponding equations for B. 

8.4- Generalized loop equations 

Let us write the equations which follow from the commutation relation A] = 1. 
We first specialize to m = n. Then 

-^n,n+l^n+l,n ^n,n— l,n 1 ; 

and thus 

Bn,n+lAn+l,n = U + 1 . 
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Writing the other equations and after some tedious algebraic manipulations, one finds 
generalized loop equations 

7V-1 /N-1 \^ 

[A{z - B)-' + {z- = - B)-^ 

n=0 \n=0 / (8.17) 

+ {z - B);^i{z - B)-^ . 

0<m<N-l<n 

This equation should be understood as a generating function for the set of equations 
obtained by expanding in powers of 1/z. 

Matrix model. Let us verify this result in the matrix model. The ingredients we need 
are the transformation law of the measure when ^ ^ Af ^ + £(Af We find 

- n n - ( 1 + ^ (E + E E >^>^~'"] 

i i<j \ i 1=0 i<j 

It is convenient to rewrite the last factor, using 

E "^r + E E ^^r-' = I (E >=^r + e e m^-^O ■ 

i 1=0 i<j i 1=0 i,j 

The corresponding variation of the integrand is the sum over eigenvalues, and we need 
only the variation for one eigenvalue. Let us call p(A) the integrand concentrating only on 
the dependence on A = A'^^^ then 

6p{X) =eAV(A) . 

We wish to calculate the integral of Sp multiplied by a polynomial Pn(A). Integrating 
by parts, we find 

J dAAV(A)Pn(A) = -J dAp(A) (A'=P.(A))' 

= -J dXp{X)B';i^A^iPi{X) . 

Gathering all terms, we recover the expansion in powers of I/2; of the loop equations ( |8.17|) . 
Note that the proof starting from the matrix model does not depend on the potentials 
being polynomials, in the same way as the proof starting from the commutation relations 
does not depend on the locality of B and A. On the other hand, the existence of the 
resolvent appearing in the loop equations implies some non-trivial topological properties 
of the matrix B. As we have seen in the one-matrix case, these have non-perturbative 
consequences and may be connected with the existence of a continuum limit. 
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8.5. Discrete form of KdV flows 

We assume now that the matrices B, A depend on a parameter t. It follows that the 
polynomials Pn depend also on t. We again set 

dPm 

rs^ = '^mnPn ■ (8.18) 
Ot 



It follows that T,mn = for n > m. Differentiating eq. ( |8.15| ) with respect to t and using 
the linear independence of the PnS, we find 

f ^[E,B], ^^E,^]. (8.19) 

For a set of parameters ti and corresponding matrices (|8.18|) implies that 

d,Ej-djE,+ = . 

Eq. ( |8.19|) thus defines a set of commuting discrete KdV flows analogous to the flows ( |3.20| ) 
of the continuum differential equations, and thus far the argument has been identical to 
the one-matrix case eq. ( ^.9|) and following. Unlike the one-matrix case, however, it is no 
longer obvious from these considerations when the matrices E are local. It is easy to verify 
that the flows associated with the matrix model are all local. We shall show below that in 
the multi-matrix model, as in the one-matrix case, non-local flows correspond instead to 
modiflcations of the boundary conditions of the recursion equations satisfled by the matrix 
elements of B, generating matrices B which no longer correspond to a matrix model. 
However, an interesting problem remains: Can all flows generated by local matrices which 
preserve the locality of the matrices A and B be associated with a matrix model? 

Explicit form of the flow generators in the multi-matrix models. If we assume again 
the set of equations (|4.8|) then we can flnd an explicit expression for the matrices E. Note 
that in what follows, to simplify notations, we shall set Ca = 1 (this corresponds just to 
a rescaling of the matrices M„ or Ba). Let us call E*^^) the matrix generator acting on 
Ai,Bi, 

^=[eW,BJ, ^=[E<'U,]. (8.20) 

In the same way, starting from the variation of the polynomials Pn associated with the 
matrix Bq-i, we can deflne a lower triangular matrix E^'?-!) which generates the flow of 

^g-l, Bq-i. 
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We now differentiate eqs. (^4.8|) with respect to the parameter t. The first equation 
yields 

We introduce a matrix E*^^\ 

ot 

Using the commutation relation [i?2, -Bi] = 1, we find 

^ = p(^\5i], ^ = p(^\5i]. (8.21) 

This argument can be repeated after introducing matrices E*^") a — 3, . . . , g. The 
result is 

^ = , ^ = . (8.22) 

Comparing with the definition of E^"?"^) above, we also have 

S(9) = (e('?-i))^ + X, 
where X is a matrix that commutes with Aq^i and Sg-i- When X vanishes, we can write 

eW + (E('^-i))^ = |5^K(S.), (8.23) 

a=l 

which determines the matrices and E^^-i) because they are both lower triangular, 
except for the diagonal elements for which only the sum for the two matrices is fixed. 
This was to be expected since an opposite change in the diagonal elements of these two 
matrices corresponds to a trivial change in the normalization of the polynomials P„ and 
P^. fJ'nPm Pn ^ fJ'n^Pn- K is easy to Verify equation (|8.23| ) in the matrix model. 

As in the one-matrix case, a flow generated by a matrix X commuting with A and B 
corresponds to a variation of the boundary conditions of the recursion relations. It remains 
to be shown that, as in the one-matrix case, such a matrix is necessarily non-local. This 
follows from the result proven in the app. D.2. 

Therefore the main difference from the one-matrix case is that we have not been able 
to characterize all local solutions of the commutation relations and all flows generated by 
local matrices as being related to a multi-matrix model. Note finally that the commutation 
relation A] = 1 does not yield additional constraints on E because the commutator 
[B, A] is a flow invariant. 
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9. The 0{n) matrix model 

In this section we employ the techniques developed in earlier sections to the case of 
the 0{n) matrix model coupled to gravity. Many of the results for the latter model provide 
interesting analogs for the other cases considered earlier here. (For early work on the 0{n) 
matrix model, see |1T7[] .) 

The partition function of the 0(n) matrix model is given by an integral over n + 1 
hermitian N x N matrices, n matrices Ai and a matrix M: 

dMd^....d^„e-»l^^<-^? + -^-*"' + '^(^'^'l, (9.1) 

with V{M) a general polynomial potential. 

The free energy of the 0{n) matrix model F = \n.Z can be interpreted as the partition 
function of a gas of loops, each indexed by an integer z, i = 1, ...,n, drawn on a random 
lattice of the form of a Feynman diagram |T^,0. In the special case n = 1 and for a 



specific class of cubic potentials V{M), the model can be shown to be equivalent to a two- 
matrix model of the form considered in sec. ^, representing an Ising model on a random 
triangulated lattice. The model cannot be solved exactly in the general case; we shall show 
however that it can be solved in the large limit, i.e. on the sphere, by steepest descent. 

The corresponding model on regular lattices can become critical only for values —2 < 
n < 2 of n mainly associated with non-physical symmetry groups (here the integral is only 
even defined for n < 2). It is thus convenient to set n = —2 cos 6. Although we could 
restrict to the interval < ^ < tt, it is convenient for book-keeping purposes to consider 
all positive values of 6. Note that the case n = 0, 6* = (2m-(- l)7r/2 reduces to the standard 
one-matrix model. 

The integral over the matrices Ai is gaussian and can be performed to give 

Z = J dM [det(M ® 1 + 1 ® M)]-i/2e-(^/3)*=^^(^) . (9.2) 

We can then parametrize M in terms of a unitary transformation and its eigenvalues A^. 
After integration over unitary matrices, the integral ( |9.2| ) becomes 



Z^jA'{A)l[{X, + X,)-/'l[dX,e (^/^^^^^^) (9.3a) 
f -iVE[A] 

= dXe ^ , (9.36) 
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with the effective action 

p. i. r/ie saddle-point equation 

In the planar limit N oo, Z can be calculated by the steepest descent method. The 
saddle point equation is 

We introduce the density of eigenvalues p(A) = jj J2i Hilbert's transform, 

P(A) 



dA- 



z-X 



the trace of the resolvent. In the large limit, p(A) becomes a continuous function and 
ujq becomes a function analytic except when z belongs to the spectrum of M, i.e. has a 
cut on a segment [a, b] of the real positive axis. 
Eq. (|9.4|) may be written in terms of ujq as 

uJo{X + iO)+uJo{X-iO)+nuJo{~X) = ^V'{X) , (AG [a, 6]). (9.5) 

This linear equation has a polynomial solution, 

^r{z) = --^{2V'{z)-nV'{-z)) . (9.6) 

Note that the cases n = ±2 are special and must be examined separately. The function 
a;(z), defined by 

UQ = Ur + uj/g , (9.7) 
then satisfies the homogeneous equation 

w(A + zO)+w(A-zO)+nw(-A) =0 . (9.8) 



Since Uoiz) behaves as for 2; large, uj{z) has the large z expansion 

1 

4 - n2 



iv{z) = --^{2V'iz)-nV'{-z)) + ^ + 0{z-') . (9.9) 
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A quadratic relation. We introduce the following function: 

r{z) = uj'^{z) + u'^{—z) + nu{z)u{—z) , 



(9.10) 



and verify that the discontinuity on the cut of r{z) vanishes as a consequence of equation 



r{z + iO) - r{z - iO) = [uj{z + iO) - uj{z - zO)] 

X [u{z + zO) + u{z - zO) + nuj{-z)) = . 

Therefore r is an even function, analytic in the whole complex plane. The behavior of u 
for z large implies that r is a polynomial. As in the usual one-matrix case, equation ( |9.10| ) 
can be directly derived from the saddle point equation or the loop equation. An expression 
for r{z) in terms of the potential follows. 

The one-cut solution: a useful representation. We introduce two auxiliary functions: 



2 sm U \ / 



2sin^ 

i 

2sTn^ 



uj-(z] = \e-'^/''uj(z)-e'^/^u(-z] 



(9.11) 



such that uj^{—z) = uj-{z), and we have 

u!+{z) uj- {z) = r{z) . 

Conversely, uj{z) is given in terms of W-i-, uj- by 

uj{z) = - (^e*^/V(^) + e-*^/2w_(^)) , (9.12) 

and eq. (^]^) is equivalent to the simple relations 

uj^(z - iO) = e^'^^u^iz + iO) , (9.13) 

which, by themselves, imply that uj{z) has cuts only on the positive axis. 

For generic values of n, the Riemann surface of u){z) has an infinite number of sheets. 
Only for the exceptional values n = —2 cos(7rp/q'), i.e. for 9 = 7cp/q, p, q being two relatively 
prime integers (and thus e*'^^ = ±1), is this number finite. Then the function u!{z) is the 
solution of an algebraic equation of degree q with polynomial coefficients. 

The case e**^^ = 1. This implies 9 = Tvp/q, where p is an even integer. The function 

s{z) = 1(^1 +u'L) (9.14) 
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has no discontinuity on the cut, and is analytic in the whole complex plane. It is therefore 
an even polynomial of a degree determined by the degree of the potential . We have thus 
the two following algebraic equations: 

LO+uo_ = r{z), ^1+1^1 = 2s{z) . (9.15) 

The solution of these equations is Uj^ = s ± V^, with 

A = s2-r'?, v^=i(w^-u;i). (9.16) 

Note that \/~K is thus an odd function. The function uj is given by eq. ( [9.121 ). 

The case e*'^^ = —1. This implies 6 = Tcp/q, where p is now an odd integer. The 
expressions are quite similar, but the role of s{z) and -\/A are formally exchanged. It is 
now the function 

«W = ^K-^-) (9-17) 
which has no discontinuity on the cut, and is analytic in the whole complex plane. It is 
therefore an odd polynomial. We have the two algebraic equations, 

uj+uj_=r{z), uj\ - uj'L = 2is{z) , (9.18) 

with solution Uj^ = v^A ± is, where 

Here -\/A is now even, and the function uj remains given by equation (|9.12| ) . 

One-cut solution. We still have to determine the coefficients of the polynomials r, s, A. 
They can be found from the additional condition that uj has only one cut [a, h] on the 
positive real axis. We can see from (|9.12| ) that the singularities of uj are the single roots 
of A. We demand that except for a and h (and —a,—b by parity), all the roots of A are 
double, such that A can be written: 

A = -{z^-a^){z''-b^)R-'{z) , 

where R{z) is an odd or even polynomial depending on the different cases. Due to the 
special form of the conditions (|9.13|) , all one-cut solutions u!± can be factor ized: 

z'^ — a?){z'^ 5 

where A and B are even functions, in general rational fractions because 0±(2;) may have 
zeros. The function O, which has only singularities at ±a, ±6, is a "minimal" solution of: 

0+(z) = 0_(-z) , Vt±{z- zO) = -e±^^0^(z + zO) . 

This factorization property is a consequence of the algebraic equation satisfied by uj{z). 
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9.2. Critical points 

A critical point is generated by the confluence of two different zeros of A or (and this 
is new with respect to the usual one-matrix model) when a cut endpoint (called a in what 
follows) approaches the origin where the integrand is singular. Let us examine the different 
possible situations: 

(i) a non-vanishing zero of A coalesces with the cut endpoint a: This is the case 
of an ordinary critical point of the one-matrix model, the determinant coming from the 
integration over Ai plays no special role, and just modifies the form of the potential. From 
the point of view of the statistical model, this is the low temperature phase in which all 
matter degrees of freedom are frozen. 

(ii) a = 0: this is a new critical point specific to the structure of the integral (|9.2| ), and 
the only case we shall consider from now on. The condition a = implies the divergence of 
N{tr Af), which characterizes the matter fluctuations. Indeed this quantity is proportional 
to Yli j + -^j) diverges only when some eigenvalue of M vanishes |T^. This 
argument is confirmed in the continuum limit by a determination of its scaling properties. 

Finally a general critical model in the continuum limit is obtained when both con- 
fluences occur simultaneously: a cut endpoint and some zeros of A approach the origin. 
Note, however, that in this limit the eigenvalue distribution approaches a singularity of the 
integrand. In such a situation the validity of the steepest descent method is questionable. 
The Ising model provides a useful test of the method. 

The resolvent. We thus consider only critical points for which a = 0. The function 
u:{z) at a critical point has a cut for < z < b. The general form of such a solution is 

uj±{z) = (^■^l-h'^/z^±ih/z^ ^^'^ (^A{z)^/l^^by^±ibB{z)/z^ , (9.19) 

where /, q are relatively prime integers with < I < q, and A, B are polynomials which 
can be chosen even without loss of generality. Indeed the situation A, B odd is equivalent 
to A, B even with the change I ^ q — I. 

It follows immediately that r(z), s{z) and ^{z) are polynomials of a form consistent 
with a one-cut solution, provided A and B vanish at ^ = 0. 

A minimal realization of a critical point with polynomial potentials is: 

-l/q 



uj±{z) =Ti{z/2bf'^+^ (^^/r^Wl^±ib/z^ \ (9.20a) 
u;± (z) = (726)2"^+2 (^^l-byz^± ib/z^ '"^^^ , (9.206) 
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where in both cases we have 



uj+{z - zO) = e^^(i-^/'?)tJ_ (z + zO) =^ e = 7r(l - l/q) . (9.21) 

The two cases q — I even and odd correspond to the two situations e**^^ = ±1. 
For z — > we find that 

uj{z) oc zP/'i , 

where in the two cases ( |9.20a, ^ ) we parametrize respectively 

case (a) : p = (2m + l)g — I 
case (6) : p = (2m + l)(j + / . 

A simple scaling argument shows that the same result would be obtained in the (p, q) 
string model of sec. ^ for the trace of the resolvent of the operator Q. 

Note that the values of p are such that m can also be defined as the integer part of 
p/2q since p/{2q) — 1 < m < p/{2q). Finally we see that for book-keeping purposes it is 
convenient to assign the angle 9 = np/q to the critical point characterized by the integers 

{p, q)- 

9.3. Scaling region 

We now wish to derive uj{z) in the scaling region, where the variable x = 1 — g/ Qc, 
which characterizes the deviation of the coupling constant from its critical value, is small. 



Functions that satisfy equation (|9.8|) and are singular only at z — ±a and z = oo have the 
general form: 

Usc,± = (^Va^ - ±iz^ ^''^ {C{z)^ - z'^ ±izD{z)^ , (9.22) 

where again C and D are even polynomials. A comparison between the large z behavior 
of and the small z behavior of u at the critical point yields the degrees of C and D. 
This determines them completely only for the minimal critical points m = 0, for which we 
find 

Wsc,± fx (^a/o^ - 2;2 ± iz^ 

To obtain the relation between x and a and completely determine the form of the poly- 
nomials C, D for multicritical points (m > 0), we calculate the deviation from the critical 
form at leading order for x small. 
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Deviation from the critical form at leading order. We now calculate the deviation 
from the critical form at leading order in the variable x = 1 — g/gc- We normalize the 
potential in such a way that for z large, the variation du!± is 



duj± ~ ^ixb/z . 

Introducing the function 

d 

^{z) = -^{gujoiz)) 



m = ^ . (9.23) 



we find that eq. (|9.7| ) implies 



It thus satisfies the homogeneous equation 
We also introduce the decomposition 

n±{z) = (e^'^/^n{z) - e^'^/^n{-z)') (9.24a) 

n{z) = - (^e'^/^n+{z) + e-^^/20_(2)) . (9.246) 

From the definition of ljq, we infer the behavior of fl for z large, Q{z) ^ 1/z. Moreover 
since O is the derivative of a function which has singularities of the form (z — zq)^^'^, 
zq = ±a, ±6, it can have a stronger singularity of the form (z — zq)~^I'^ . These conditions 
determine 0(z) uniquely as a function of the location of the singularities. The variation bu) 
is proportional to O at the critical point. At the critical point, has the form ( |9.19|) . To 
obtain its complete form, we need its small z behavior which must be consistent with the 
leading correction to the large z behavior of Wgc- Since O is independent of the potential, 
we can compare it to the form of Wgc in the case (a) (see preceding subsection) for m = 0. 
The leading correction to is then of order z'/'^"^, and the unique solution is 

z^l-h'^jz'^ ' ' J 

With this information we can now explicitly calculate the scaling functions for all critical 
points. 

Scaling function. From the preceding analysis, we conclude that for z large Wsc satisfies 

- const ^P/^ = O [xz^f'i-^^ . (9.25) 
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It is easy to verify that this fixes the polynomials C, D. We now show that the scaling 
function uj^c can be expressed in terms of the function d{z) given by the integral represen- 
tation 

^{z)= dt{t + izf/''-"'-\t-iz)'^, (9.26) 



the proof relying on a verification of condition (|9.25|) . Calculating the integral ( |9.26| ) we 



first verify that the function 'd^z) has a form consistent with expression (|9.22|) : 
^{z) = (Va2 -z^ + iz^ ^^^'^ (C{z) - + izD{z)^ , 

where C, D are two even polynomials of degree 2m and ±l/q = p/q — 2m — 1. 

Another representation of the function is useful. Setting z = a cos{q(p) and integrating 
over t, we find 

p/q- m f^^ \ r J 

Moreover if we introduce the parametrization of case (a) of the preceding subsection, 
p = (2m + l)q — I, we obtain 



d{l,z)= I dt(t + Z2)(P-'?-^)/(2'?)(t-Z2)(P-'?+0/(29) . 



It follows that 

z) - -z) = {2z)P/'' ^ 

where we have set 

apq - B{m + l,p/q-m) - — — ■ — . (9.28j 

T[{p + q)/q) 

Now we need to expand 'd{±z) for z large. 

Large z expansion. For z = —iX large, we find 

^(_,A) = ^_i)m^2xr/, nm+mp/Q-m) ^ ^2A)p/'^-2— 2^!!:i! 

T{p/q + l) m + 1 

+ O (ap/'?-^'"-^) , (9.29a) 

X'iSm-p/g 2p/(j-2m 

mX) ~ ^ — . (9.296) 

p/q — m 
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Case (a) . If p = (2m + 1)^ — /, and thus p/q < 2m + 1, 'd{—z) is asymptotically larger 
than the correction to 'i}{z). Moreover 2m — p/q = l/q — 1. The solution is then 



u,,,±{z) = ^o^{±z) . (9.30) 

Moreover comparing the expansion ( |9.29| ) with the expansions of the critical functions iv± 

and duj±, 

LO+i-iX) ~ (-l)'^+i(A/26)P/5 
6uj-{-iX) ~ -x2^-^/i{X/by/'i-^ , 
we obtain the normalization constant and the relation between a and x: 

«„ = -(4i.)->'/V-i, (9.31a) 



where we have used the definition (|9.28|) . 

Case (b). If p = (2m + 1)^ + /, and thus p/q > 2m + 1, the correction to 'd{z) is 
asymptotically larger than the correction to 'd{—z). Moreover p/q — 2m — 2 = l/q~l. The 
solution is then 

uJsc,±{^) = ^o^iTz) , (9.32) 

where, in our normalizations and in the set of variables p,qj, the same relations ( p.31j ) 
continue to determine 'do and the relation between a and x. 

9.4- The singular free energy 

We can find the singular part of the free energy. We have shown that 

^ ''g^^) = ^ (fdzV{z)n{z) , (9.33) 



dg \ dg J 2in 



where Q{z) is the function ( |9.23| ). Using the decomposition (|9.246P , we can rewrite equation 

as 

The critical function for a = is 



l-l/q 
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The scaling region. Let us now consider the case a ^ but small. For z small, the 
function ( |9.35| ) behaves as 2;^/'^"^. This, together with the other properties, determines the 
scaling form of Q{z): 



Va2 - ^2 V " ' y ' sin(^/2) 

Conversely, the next-to-leading term in the large z expansion of O provides the additional 
information needed to completely determine the first correction to the critical function for 
a small. This correction behaves like z~^~^l'^ , and therefore 



0+(a) - 0+(a = 0) oc U\ - z^ + ih/z 

z^\l\ — h-^ jz^ ^ ^ 

The leading correction to Osc,+ (^) is 

^sc,+ (-^A)~^?o2-'/V'/^A-l-'/^ , 

and thus 

0+(a) - O+(0) ~ —7^ [l) , ,^ , (Vl - hVz^ + ih/z) 

sm(6'/2) V6/ z^^l-lP'/z^^ ^ 

The identity 

d / , , \ .J 1 / ; TTT^ , \ 



- - h^/z^ + ih/z) = ib- y/1 - byz^ + ib/z 



allows us to cast this expression into the form 

With this expression we can integrate by parts in the integral (|9.34|) , giving the second 
derivative of the singular part of the free energy as 

2i7Tgclsme/2\4bJ J \^ ' ' 



We now substitute the identity 

e^e/2y'(^) +e-*^/V(-z) = 2z(7sin^ {^-'^ 1'^ uq{z) - l'^u^{-z) 



-K4sin^6'w_(z) 
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where we need U- only at leading order, as given by the expressions (|9.20| ). We see that 
the contribution to the integral coming from u}- vanishes. The contribution due to ujq can 
be calculated by taking the residue at infinity. Then only the leading behavior of ujq for z 
large is relevant, 

Ug sin6'cos(6'/2) 



2ig sin 9 {e 



In terms of the variable x = 1 — g / gc finally obtain 

From this expression, we derive the scaling of the free energy for all critical points: 

9/ 

F,';.(x) oc , ^ 7str = -— (9.36) 

^ p + q — I 

Discussion. We find a result in agreement with the (p, q) string models only in the 

case / = 1, i.e. when p is of the form p = (2m + ± 1. We note that this is the particular 

class of models for which we know the operators P, Q at leading order (see subsec. \4-^ - 

The trace of the resolvent of Q can be easily calculated. In the semiclassical limit we have 

indeed (eq. ([4.55| )) 

u{z, x) = P(d, x), Q(d, x) = z . 
In the case p = (2m + l)q'±l, Qisa Chebychev polynomial, 

Q{d,x) = u''/\x)T,{d/V^ , 

and P is given by (eq. ( [4.53|) ) 



P(d, X) = U^l^ {^'f\Tv-1rMI^) ■ 



Setting z = 2^*^/^ cos(g(/?), we see that cu, up to normalizations, agrees with the results in 
eqs. ( |9.26| , |9.30| , |9.32|) , where d is replaced by the form ( |9.27|) , obtained for the correspond- 
ing (p, (?) critical points of the 0{n) model. 

We conclude that on the sphere we have obtained in the case p = (2m + l)g ± 1 
complete agreement between the results of the (p, q) string models and those of the 0{ji) 
models. For other cases they seem to differ. One possible interpretation is that the most 
relevant operator of multimatrix models is not present here. And indeed one can find in 
the (p, q) model another relevant operator which has the proper dimension to yield the 
value of 7str found here. Introducing the parameter t coupled to this operator (which has 
the dimension x^'^^'^~^^ I one finds that the resolvent c<;(2, of the (p, string 
model coincides indeed for x = with the resolvent of the Oin) model [|7T 
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10. Open problems 

Many questions that require further investigation have been pointed out in the text. 
For example, the solutions of the differential equations which have been derived should be 
studied further. The perturbation expansion has been shown to be non-Borel summable 
in many cases of interest, and the solutions of the differential equations, having unphysical 
properties, do not lead to a solution beyond perturbation theory. In the pure gravity case, 
a candidate for a non-perturbative solution has been proposed based on the solution of a 
Fokker-Planck equation . 



It is possible to integrate over the relative unitary matrices for multi-matrix models 
in which the interactions among the matrices involves no closed loop. The method of 
orthogonal polynomials, however, has only been applied to the case of matrices interacting 
along a line. 

It will be useful to understand the underlying relation between the KdV and Liouville 
approaches to these theories in all generality, rather than the current successful comparison 
of partial results for low genus and specific correlation functions. 

Finally, the matrix models considered here representing d < 1 matter coupled to 
gravity are quantum mechanical models with a finite number of degrees of freedom. For 
application to more realistic string theories, it will be necessary to extend these techniques 
to the case of d > 1 matter coupled to gravity, which would correspond to a quantum 
field theory of matrices. The borderline case of d = 1 matter coupled to 2D gravity, 
realizable as a solvable quantum mechanical matrix model, admits an interpretation as 



a critical string theory with a two dimensional target space (for review, see |]Tl|). This 



model already provides many tantalizing hints of interesting physics that emerges as we 
approach the field theory case. The spacetime interpretation of these theories and the role 
of the Liouville "time" will greatly enhance our understanding of physical realizations of 
string theory. 
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Appendix A. KdV flows and KdV hierarchy 

We consider the Schrodinger operator — L 

L = d^ — u{x, t) , 

where in what follows d means d/dx, and t parametrizes the potential. We introduce a sec- 
ond (anti-hermitian) differential operator M of order 21 — 1, and look for the compatibility 
condition between the Schrodinger equation 

-mm = Em , 

where E is assumed to be t independent, and the condition 

dtip = -M{t)iP , (A.l) 

which then implies that the Schrodinger operators corresponding to different values of t 
are unitary equivalent. 

Note that M can be written 

I 

M = $^{m,(x,t), d^^-i} , 
k=i 

where the symbol {, } means anticommutator, and thus a priori depends on I independent 
functions. 

The compatibility condition for this linear system, 

[dt + M, L + E]^ = Q , 

implies that 

{atL-[L,M]}V = 
for all eigenfunctions ip. It follows that 

dtL=[L,M]. (A.2) 

A.l. The resolvent. Local conserved quantities 

Since all operators L{t) are unitary equivalent, any function of the eigenvalues is 
conserved in the flow ( |A.1| ) or ( |A.2| ). Let us consider the resolvent 

G{z) = (z + L)-^ . 
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Then we have 

dtGiz) = -(z + L)-^dtL(z + L)-i 

= {z + L)-i [M, L] {z + L)-^ - [G{z),M] , 

and thus taking the trace of both sides, we find 

dttTG{z)^0. (A.3) 

By expanding tr G{z) for z large, 

tvGiz) - -E / dxRk[u]{-z)-'-'/^ , i?o = I , (A.4) 

we generate an infinite number of conserved quantities which we exphcitly construct below 
and show to be the space integrals of local differential polynomials in u. 

Explicit construction of the conserved quantities. Let us consider the Schrodinger 
equation for the resolvent G{z) = {z + L)~^, 

(2; + d^ — u{x)^G{z] x,y) = 5{x — y) , (A. 5) 

where z does not belong to the spectrum of —L. 

Let us recall how G{z] x, y) can be expressed in terms of two independent solutions of 
the homogeneous equation 

{d? -u{x) + z)^i,2^Q . (A.6) 

If we partially normalize by 

and moreover impose the boundary conditions 

ip\{x) — > for a; — > —00, V2{x) — > for a; — > +00 , 

then it is easily verified that G{z; x, y) is given by 

G{z;x,y) = (pi{y)(p2{x) d{x - y) + (pi{x)(fi2{y) 9{y - x) . 
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The diagonal matrix elements r(x) = G{z; x, x) = (pi{x)(p2ix) satisfy 

r'(x) = ip[{x)ip2{x) + ipi{x)ip2{x) , 
r"{x) = 2[ip'i{x)ip2{x) + [u{x) — z)r{x)) , 
r'"(x) = 2(u'(x)r(x) + 2{u{x) - z)r'{x)) , 



where we have systematically used the Schrodinger equation (|A.6|) . 
The last equation, involving only r(x), is usefully written 

— zr'(x) = |r'"(x) — u{x)r'{x) — |w'(x)r(x) = j (d^ — 2{w, d})r . (A. 7) 

Let us now use this equation to expand r(x) for z large and negative. We can obtain r(x) 
at leading order for z large, directly starting from the definition ( [A .51) , and applying the 
WKB method (in this limit the non-commutation between d and x can be neglected). We 
find 

If dp 1 1 

r{x) - ' — 



271 J z-p'^-u{x) ^ ^-z + u{x) 

= + ^(-.)-^/^ + o ((-.)-^/^) . 

It follows that 

G{z-x,x) = -^Rk[u]{-z)-^-^/^ (A.8) 

for any smooth potential 'u(x), where Rk[u] is a local functional of u. Equation ( [A.7| ) leads 
to the recursion relation (|3.14| ) for the i?fc's. 

Remark. The conservation equation ( |A.3|) shows that as a consequence of the flow 



equation ( |A.2|) , dtRi can be written as a total derivative. 



Two useful properties. 

(i) G{z) defined by ( |A.5|) satisfies 



iTG{z) = {x\{L- z)-^\x) = -—G{z;x,x) . (A.9) 



5u{x) dz 
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Eq. ( |A.8| ) relates the quantities Ri to the expansion of G{z; x, x) for z large. An expansion 
of ( |A.9D for z large yields immediately 



du 



j dxRi[u] = -{I - \)Ri-i[u] . (A.IO) 



(ii) Let us denote by d/du the derivative corresponding to constant variations of u. 
Then the resolvent satisfies: 

Again an expansion of this identity for z large yields 

— R,[u] = -{l-\)R,_M. 

A. 2. The flow equation: discussion 

Let us examine the content of the flow equation ( [A.2|) . The l.h.s. is no longer a 
differential operator, while the r.h.s. is at most a differential operator of degree 21. Iden- 
tifying the coefficients of all independent operators yields / differential equations and one 
partial differential equation (hermiticity being taken into account) for the / coefficients 
nik and u. The operators M can then be constructed using the algebraic technique of 
pseudo-differential operators. 

Pseudo- differential operators. One considers non-integer powers of the operator 
L. Formally, the operator L'^ may be represented within an algebra of formal pseudo- 
differential operators as 

oo 

L'^ = d'"-|-^e,(a;)d'"-^ (A.ll) 

The choice of putting all pseudo-differential operators on the left is arbitrary. The equiv- 
alent expansion with operators on the right can be obtained by using the commutation 
relation: 

The successive terms of the expansion of L'^ can be obtained for instance by writing L'^ 
as a Laplace transform. 

Explicit construction of the flow equations. The formal expansion of L'~^/^ (an anti- 
hermitian operator) in powers of d is given by 

^1-1/2 ^ ^21-1 _ 2[_1 ^2Z-3| ^ . . . ^A.13) 
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(where only symmetrized odd powers of d appear in this case). We now decompose 
ij-^l'i' = _|_ ij_ 1/2^ where L^,_ = d^'~^ + ■ ■ ■ contains only non-negative pow- 

ers of d, and the remainder has the expansion 



l'-"" - E {^2.-1, d-(2^-^) } = {Ri, d-'}+ 0(d-3) + ■ ■ ■ . (A.14) 

i=l 



Here we have identified Ri = ei as the first term in the expansion of For LV2, for 

1 /2 

example, we find = d and Ri = —u/A. We justify below that Ri defined by equation 



( |A.14| ) is also the functional of u which appears in the large z expansion of tr G'(2;) in ( [A.4| ). 
Since L commutes with L^~^/'^, we have 

= [L,lL"'/'] . (A.15) 

But since L begins at d^, and since from the l.h.s. above the commutator can have only 
positive powers of d, only the leading (d^^) term from the r.h.s. can contribute, which 
results in 

[L^^^/^ L] = leading piece of [L, 2Ri d"^] = AR'i . (A.16) 

Z — 1 /2 

If we then take for M a linear combination of the operators L_|_ , 

i-i 



the equations for the coefficients nik are automatically satisfied and we find an equation 
for the potential u: 

d ^"^ 

dtu = -A—^ iikRk+i[u] . 

k=0 



A recursion relation. The quantities Ri in ( A.16 ) are easily seen to satisfy a simple 



recursion relation. From L^+^Z^ = LL' = ^/^L, we find 
Commuting both sides with L and using ( [A.16|) , simple algebra gives 



R'l+i - -R'l" - uR'i - -ju'Ri , 
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where we recognize equation (|3.14 ). While this recursion formula only determines R[, by 



demanding that the Ri {I ^ 0) vanish at m = 0, as implied by the definition, we obtain 

„ 1 „ 1 „ 3 2 1 // 

^0 = ' ^1 = -7^ ' = - —u , 

2 4 lb lb 

i?3 = + - iuu" + \u'^) - , (A.17) 

= _ ^ . ^ ^2^//X ^ (2^,^,(4) + ^u'u'" + 3W"2) _ ^^(6) . 

25b 128^ ^ 25b^ ' 25b 

/ — 1 /2 

We summarize as well the first few , 

rl/2 J r3/2 i3 3. 

^+ = d , = d^ - -{w, d} , 



L'^^ = d'-\{uA'} + Y^{{^u' + u")A] 
7 

bi 



(A.18) 



^ I (I3w(^) + lOW + lOu^ + 25w'2) , d 
34 I 



Examples of flow equations. The simplest equation only involves -Ri and is thus linear 

dtu — niu . 

the unitary transformation simply corresponds to the translation of the coordinate x. The 
second equation is non-linear instead and is actually the original KdV equation 

d^u = ^ [u'" - &uu') . 

A. 3. Large z expansion of the resolvent: residue and trace of pseudo-differential operators 

We now relate the large z expansion of the resolvent to two operations called "residue" 
and "trace" of pseudo-differential operators. To calculate the coefficients of the expansion 
of G{z]x,x), we assume for convenience that the spectrum of —L is strictly positive and 
consider the integral 

I{s,x) = -\ j dzz'~^G{z]x,x) , (A.19) 

where the integration contour can be taken for example going from —zoo -\- iO to -|-zoo -f- zO. 
The residue of I{s,x) at its poles located at s = k + 1/2 is the coefficient of z~^~^/'^ in 
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the expansion of G{z]x,x). After the change of variables z ^ —Lz and integrating over 
z, (|A.19| ) becomes 

I{s,z) = in{x\{-Ly-^\x) , (A.20) 

where (a;|(— L)'^~^|x) is defined by analytic continuation from s < |. The poles of the 
expression are related to the short distance, large momentum singularities. 

If we expand L^~^ as a formal series in d, poles in s will appear when one power of 
d approaches —1, the usual logarithmic divergent term. A short calculation then shows 
that the coefficient of the expansion of G{z; x, x) is given by the residue 

of L'^"^/^, where the residue of the pseudo-differential operator A = XliL-oo '^«('^) 
defined by 

Res A = a-i . (A.21) 

Then we have 

' ' k=0 



Similarly the trace of A is defined by 

trA = y da; Res A = j dxa_i . (A.22) 

The trace trL^+^Z^ is thus the coefficient of {—z)~^~^^'^ in the expansion of trG(z), 

tr G{z) = -i^(-^)-^-V2trL^-i/2 , 

k=0 

justifying the nomenclature and in particular the cyclicity property!! tr AB = trBA for 
any two differential operators A, B. 

We now give two simple examples illustrating the usefulness of these operations. 

(i) Since 



The cyclicity property can also be verified directly by considering basis elements a{x)d"^ and 
6(x)d". We find 

trad""M" = /EJlo (7) Resa5(^)d™+"-^' = / {^^^+^) a5(""+"+^). Integration by parts then 
gives 

trad""6d" = / 6a(""+"+^) = tr M^ad"', where we have used the identity 

(m+"+i)' following from the definition ( j) = k{k — 1) . . . {k — j + 
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we can take the trace of both sides. The cychcity of the trace imphes 

SttvL'-^/^ = , 

and we thus recover the conserved quantities ( |A.4| ). 
(ii) Since Ri+i = \ ResL^+^/^, we see that 

= -(i + i)B,|u] , 

recovering the identity (|A.10[) . FinaUy, introducing a polynomial n(L) in the variable L^/^, 
such that M = n_|_(L), we can write the flow equation (|A.2| ) as 

d 

dtL = \L,U+(L)] = -2— Resn(L) . (A.23) 

ox 

A. 4. Commutation of higher KdV flows. 

Let us consider the set of flows associated with the operators M = L^^^^ and call ti 
the corresponding parameters. We show now that all these flows commute. We have to 
prove that the linear system 

D,L.g-lL,L'WV0 
is compatible. This is again a zero curvature condition 

^ [D„ D,] = d,L^:'/' - a.4+^/^ + [Ll+'/\ L^:'/'] . 
The equation can be rewritten 

and thus, using the flow equations, we can write 

([L^+V2, 4+^/2] _ [L'^+V^, 4+^/2] + [4+V2^4+V2]j^ ^0 _ 

Replacing 4^^^^ ^he above by L-'+^Z^ — L'L^^^^ for j = /c, we flnd that the above 
equation is identically satisfied, and hence the flows commute. 
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A. 5. The canonical commutation relations 



It is now very easy to relate the previous considerations to the solutions of the string 
equation [P, Q] = 1. Let us examine the compatibility of the equation, where we take P 
of the form 

K{Q) being a polynomial in Q^^^, with the KdV flow 

g=[Q,Qr/1 • (A.24) 

Differentiating the string equation with respect to tk and using the flow equation ( [A.24|) , 
we flnd 

+ = 0. 

We rewrite the second term by using the Jacobi identity and the string equation and flnd, 
as expected, the compatibility condition 

Let us now calculate explicitly the l.h.s.: 

dP dK+ fdQdK\ 

I ^77 Pin I ^ ' 



dtk dtk \dtk dQ 



where, as in the verification of commutation of KdV flows, we have used the property that 
differentiation with respect to tk and truncation of the differential part commute. Thus 
we have 

and it follows that 



In the r.h.s., q'^^^'^ can be replaced by Q'^"'"^/^. Moreover, 
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The compatibility condition between the KdV flows and the string equation is thus satisfied 
with the choice 



OK 
dtk 



k-l/2 
+ 



justifying the equations which appear in subsec. \3.3^ Note also the relation between the 
action ( p.l8| ) and the form of P 



S{Q)=tYS{Q), P(x[S'{Q)] + 



(A.25) 



Appendix B. Generalized KdV flows 

Most of the preceding algebraic results can be generalized to higher order differential 
operators 

9-1 

i=o 

The explicit form of the equations is of course slightly different. In this appendix we 
provide a few of these expressions. 



The large z expansion of the resolvent. From the analysis of app. [A. j[ , we know that 



the large z expansion of the resolvent is related to the residues of L'^ ^. The residues do 
not vanish only when s is of the form s = 1 + m/g. We thus find 



G{z-x,x)= J] ResL""/'? 



(B.2) 

trG'(z) = (-z)-i-"^/nrL"^/'J , 

m= — 1 

showing in particular that the traces tr L"^/*^ form a complete set of local conserved quan- 
tities in the fiow. Note that relations ( [B.2| ) can be verified in the WKB approximation in 
which 

dp 



.x\[z 



+ L) ^\x) = G{z\x,x) ~ — [ 



Rescaling p, we see that x, x) has an expansion of the form (|B.2|) 
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B.l. Explicit construction 

The new feature is that dtL is now a differential operator of order q — 1. The operators 
M however can again be expanded on a basis formed by the set of differential parts of 
rational powers of L. Indeed the relation 

[L,L7/'^] = -[L,L'I^/''] 

shows that the commutator on the r.h.s. is a hermitian differential operator of order q — 1 
as required. 

Note that since the proof of commutation of KdV flows did not depend on the specific 
form of the L operator, the property remains true in this more general case. Also since the 
problem has become purely algebraic through the use of the pseudo-differential operator 
formalism (at least as long as we do not discuss the existence of solutions of the equation), 
we can immediately generalize previous considerations to general differential operators L 
of arbitrary (not necessarily even) degree q. We then expand the second operator M in 
the basis L^^^. 

(p, q) string equations. We now consider the string equations [P, Q] = 1 with 

Q = d-^ + {^;g_2(x), d'?-2} + ■ ■ ■ + 2vo{x) . (B.3) 

(By a change of basis of the form Q —>■ f~^{x)Qf{x), the coefficient of d*^"^ may always 
be set to zero.) We then exf 
variable Q^/i. We can write 



be set to zero.) We then expand P on the basis Q^_l'^ . Let K[Q) be a polynomial in the 



P=[K{Q)\^ 

(in K{Q), integer powers of Q can of course be omitted). Again we can study the compat- 
ibility of the string equation with general KdV flows generated by 

The calculation is the same as in the case q = 2. One finds the solution 

^-{k/q)Q'/'^-' ^ P = J2(k/q)t,Q'J'^-\ 



consistent, up to sign conventions, with ( [4.22| , [4.23| ). Comparing with eq. ( |4.24| ), we also 
see that the relation (|A.25| ) is valid for all (p, q) models. 
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B.2. {p^q) and {q,p) actions 



In subsec. ^Tj, we have indicated that the string equations could be derived from an 
action principle (eq. ( [4.24| )). This action seems to break the initial symmetry between the 
two operators P and Q. Let us first examine this problem here in the simplest case, i.e. 
with tm = only for m^l or rriy^p + q. 

The action tr Qp/'J+i was constructed from a q order differential operator Q. We 
assume p> q, and discuss what happens if we consider instead the action trP'^/^"'"^, where 
P is a p^^ order differential operator. We shall show that the two actions tr QP'/'?+^ and 
^j-pg/p+i give rise to the same equations of motion (modulo extra integration constants 
that arise because P contains higher derivatives). 

We expand 

P = dP + {V2,d^"'}+ ■■■ +2^0 • (B.4) 

Since p > q, we can take Q = PlJ" , where the coefficients in an expansion of Q of the 
form ( [4.18| ) may be expressed in terms of the v^s of ( [B.4| ). Note however that Q will only 
depend on the q — 1 quantities Vp-2i . . . , Vp-q. 

First we eliminate Sq, . . . , Vp-q-i by imposing the equations of motion 

.trP'?/p+i = o (a = 0,...,p-g-l) . (B.5) 



These imply that Res {d", P^/p j = for ct = 0, . . . ,p — g — 1, and consequently that 

P^/^ = {ep-,+i,d-(^-'^+i)}+ ... , (B.6) 

where the additional terms involve only derivatives d~'^ with k > p — q -\- 1. 

We now show that the actions trP'^/^^^ and tr Q^/'^"'"^ are proportional after imposing 

trQP/^+i = tr [P^J^Y'"^^ = tr [p'^'^ - P^p)''^'^' 

= tr P^/^'+i - b/? + l)tr (P Pi'") (B.7) 
+ i(p/Q + l)^tr(pi-^/P(P!/^)')+ ... , 

where the neglected terms are of lower order in d either because they involve at least one 
commutator [P'^/^, P^^^] or a higher power of Pl/^ (see the appendix of ||3^]). From (|B.6|) , 



we see that only the first two terms on the r.h.s above can contribute a nonvanishing d ^ 
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term. But since ti PP'^^ = (both contain only positive powers of d), we also have 



l^PP<i/p = trPPi/P = trP9/P+^ Substituting in (^) gives the result 



trgp/g+i = -EtrP-j/P+i , (B.8) 
Q 

which shows that the two actions generate the same equations of motion. 

A general action. We shall see that if we had restricted ourselves to operators P of 
the form P = Q^'^ , the argument would have been even simpler. Let us consider a more 
general case where the action S{Q) is the trace of polynomial S{Q) of degree p + q in Q^^'^, 
normalized by its term of highest degree: 

S{Q) = tr S{Q), P = [K{Q)]+ with K{Q) = S'{Q) = Q^'^^ + ■ ■ ■ . 

A first order calculation. To explain the idea we first restrict to 

S{Q) = ^-tr (q^+p/i + t^tvQ'^/A , (B.9) 
p + q \ / 

and calculate to first order in tm- We find 

= tr {K-K_)^+''/P 

= tiK^+'i/P - (1 + q/p)trK'^^PK_ + i(l + q/p){q/p)tr K'^/p'^KI + ■■■ , 
where again only the two first terms contribute because q < p. Then we have 

and thus, since tr QK_ = tr Q{K — K^) = tr QK {QK^ has only positive powers) there 
results 

trpi+g/p = (^trQi+P/'? + (m/p)t^trQ('"-P)/'?Q^/^) . 

Let us now calculate also trP'^^P at leading order: 

trP""/P = tr (Qp/'^ -QL^'^j 

= trQ'"/'? - (m/p)trQ^/'^Q("^-P)/^ . 

We can eliminate the term tr Qp/'^Q^^~p^/'^ between the two equations. Introducing the 
notation 

S(P) = (trP^+'ifP - (q/p)tmtrP'^/P^ 

p + q \ V 

152 



we obtain 

S{P) = -S{Q)+0{tl) . (B.IO) 

This shows, at first order in tm, that there exists an action S{P) equivalent to S{Q). 

General proof. Let us now consider a general polynomial S{P) of degree p + q m P^^p. 
For convenience, we normalize it by S'{P) ~ pp/i for P large. The general proof relies on 
the identity 

tr^(P) = tr S{K) - tr S'{K)K- + ^tr S"{K)K^ + ■■■ 
= tr S{K) -trS' {K)K_ , 
because p > q. Again since trQK_ = trQ{K — K^) = tr QK, we have 

tr S'{K)K_ =trQK + tr{S'{K) -Q)K_ . 

We then choose the polynomial 5" such that 

This equation is satisfied by inverting Q = S'[K{Q)'j in the functional sense, expanding 
for Q large, and taking the polynomial part in K^^p of the solution. It follows that 

tr^(P) = tr S{K) - tr QK . 

Differentiating the r.h.s. with respect to Q, we obtain —K = S'{Q), and therefore 

tr^(P) = -tr^(Q) . 

The quantity trS'(P) is thus the action expressed in terms of the P operator. In the 
functional sense, S{P) is the Legendre transform of S{Q). 

Appendix C. Matrix models and jacobians 

In sec. H, we gave a justification for the appearance of the Vandermonde determinant 
( p.2| ) in the one-matrix integral. As usual a more powerful method to calculate Jaco- 
bians relies on determination of a metric tensor. Let M be a hermitian matrix which we 
parametrize in terms of a unitary matrix U and eigenvalues Xk- M = AU. Let us then 
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calculate the square of the line element in the new variables. We express the variation dU 
of U in terms of a hermitian matrix dT, dU = idT U. Then 

tr {dMf = tr (U^ (dA + z[A, dT])[/)^ = ^^(dAfc)^ + ^ (A, - A,)^ |dT,,f . (C.l) 

k i,j 

Note that the independent variables are the variation of the eigenvalues dAfc and RedTi-,-, 
ImdTij, for i < j. From eq. ( |C.1| ), we immediately obtain the Jacobian M A, t/ which 



is proportional to v G, where G is the determinant of the metric tensor. In the notation 
of sec. 0, we have explicitly 

G = ll{X^- Xjf = n - ^^O' = ^'(^) ' 

i^j i<j 



and we thus recover the identity ( |2.1|) . 

The expression ( |(J.1| ) has another useful application here: it yields the laplacian 

applied to symmetric functions of the eigenvalues of the matrix M. It is easy to verify that 
this laplacian is equivalent to a sum of free terms. Setting A(A)'i/'(A) = ^(A), we have 

A(A)£A-i(A) = A-1(A)Aa2(A)^A-1(A) . 

After a series of commutations, this expression reduces to 



k L k ^ ^ 

where the square brackets are meant to indicate that the second term is no longer a 
differential operator, the derivatives acting only on A. The quantity J2k (^/^-^fc) ^(^) is 
a totally antisymmetric polynomial in the A^'s. Any antisymmetric polynomial, however, 
is necessarily proportional to A. Since the quantity has a degree smaller than A, it must 
vanish identically so that 

2 



2 



A(A)£A-'(A) = 5:(^) . (C.2) 



This result has several applications. It provides a partial solution to the D = 1 
matrix problem. As long as unitary excitations can be neglected, the quantum mechanics 
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of hermitian matrices can be reduced to a free fermion problem: The ground state tj^ is a 
completely symmetric function of the eigenvalues. Because the Vandermonde determinant 
is completely antisymmetric, the transformed function (p must also be antisymmetric. 

The Itzykson-Zuber integral. The same result can be used to prove the basic identity 
for the multi-matrix model. 

The free hamiltonian —C has the plane-wave solutions: 

_^^rtr KM ^^^J^2^^tr KM ^ ^^3^ 

where K is an arbitrary hermitian matrix. We diagonalize K, 

K = V^TV , 

and integrate both sides of ( |C.3| ) over the unitary matrix V. The eigenfunction 

then remains a function only of the eigenvalues of the matrix M. Using the result ( |(J.2| ), 
we can set A(r)A(A)'i/j(A) = (p{A) and the function (p satisfies 

k ^ ' 

where (^(A) in an antisymmetric function of the A^'s and thus is superposition of Slater 
determinants of free solutions. We can therefore write 

<^(A)= /dp(^fc)dete^^'^^- with Y.i^\ = Y^^1. 

i i 

We now note that ^(A) is symmetric in the exchange A F. This fixes the solution 
up to a global normalization, with the result 

J dUe'''^'^^^ oc A-\T)A-\A) dete^^^^^' . 

If, instead of integrating over a unitary matrix, we integrate e^^^ over K, then, because 
A~^(r) is antisymmetric in the 7i's, we can replace the determinant by only one term. 
Moreover the jacobian for K \—>- V,r generates as we have seen a factor A^(r). We thereby 
recover the identity (|4.1| ). 
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Appendix D. Discrete canonical commutation relations: A few additional re- 
sults 

In this appendix, we consider the matrix solutions to certain problems posed as matrix 
commutation relations. 

D.l. A uniqueness theorem 

Let ^ be a lower triangular matrix with vanishing diagonal elements, B a matrix such 
that Bn,n+i 7^ for all n > and B^n = for n > m + 1, satisfying [B, A] = 1. Then we 
shall show that 

given B , A is unique . 

Indeed let us assume that we have found two solutions Ai, A^ of this linear equation, 
then X = Ax — A-i commutes with B: \B^X\ = 0. We consider the diagonal elements of 
this commutator, 

[-S;^]mm — — Bm,m+lXm+l,m -^m,m— l-^m— l,m • 

An induction over m thus shows that Xm+i,m vanishes. Let us then assume that Xm,m-k 
vanishes for all k < n. It follows that 

[-^) -^]m,m— n — -^m,m+l-^m+l,m— n -^m,m— n— n— l,m— n — • 

Again an induction over m shows that Xm,m-n-i vanishes, and thus finally X = 0. The 
equation [B,A] = 1 for A has a unique solution. 

Application. If B is symmetric, as in the one-matrix problem, X an antisymmetric 
matrix, then [B, X] — implies X = 0. The proof follows simply from the previous result. 
We set X = A — A^ where A is a lower triangular matrix with vanishing diagonal elements. 
Then 

[B,A] = [B,A^] . 

The first commutator has vanishing matrix elements [B,A]mn for n > m, the second 
[B, A'^]mn for n < m. Thus J — [S, A] is a diagonal matrix and 

J=[B,A] = [B,A'^] = = [A,B] = . 

According to the preceding result, A thus vanishes as well as X. 
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D.2. Another useful result 

Let X be a matrix such that Xmn = for n > m + s {s > 0), and let X commute 
with B, where B has the form of the matrix Bi of the multi-matrix model (i.e. B is local 
with Bmn = for n > m + 1, and Bm,m+i 0)- We shall show that 

X = V'{B) , (D.l) 

where V{B) is some polynomial in B. 

Let us write the set of equations [-B,^] = 0. We start from the matrix element 
n = m + s of the commutator. Only one term contributes to each product of matrices 

■Bm,m+lXrn+l,m+s ^m,m+s — l-^m+s — l,m+s ~ ; 

which we rewrite 

Xm+l,m+s -^m+s — l,m+s 



Xm,m+s — l ^m,m+l 



Multiplying the numerator and the denominator in the r.h.s. by the product Bm+i,m+2 • • • Bm+8-2,m+8-i, 
we write 

^m+l,m+s Tn+l,m+s 



and therefore 



Xm,m+8-l 



^m,m+s-l — -^s-l-°m,m+s-l • 



Next we write the matrix element n = m + s — 1. Two terms contribute and we 
consider the equation as an inhomogeneous recursion relation for Xm,m+s-2, 

Bm,m+lXm+l,m+s — l ~l" BjYi^mXm,m+s — l ^m,m+s— 2-Bm+s— 2,m+s — 1 
Xm,m+s—l-Bfn+s—l,m+s—l — • 

Because [B^~^,B] = 0, Xm,m+s-2 = satisfies this equation. The general 

solution is thus the sum of this special solution and the general solution of the homogeneous 
equation 

Bm,m+lXm+l,m+s—l -^m,m+s— 2-Srjx-|-s— 2,m+s— 1 — • 

We now use the same argument as for the earlier equation and conclude 

Xm,m+s-2 = Ks-2B^^_^_g_2 . 
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Iterating the argument, we find that X can be written 



X = V'{B)+Y, (D.2) 

where V{B) is an arbitrary polynomial of degree s and y is a lower triangular matrix with 
vanishing diagonal elements that commutes with S, 

[s,y]=o. 



In app. D.i, however, we have shown that this equation has the unique solution y = 0, 



thus establishing (|D.l 



Remarks. 

(i) It follows in particular that if in addition [^, A] = 0, with = 1, then 
V"{B) = 0, and thus X is proportional to the identity. 

(ii) If with the same conditions on the matrices B and X we look for a solution of the 
equation [X, B] = A, with Amn = for n > m, we find from the previous arguments that 
X can be written in the form (|D.2| ) 

X = V'{B)+Y , 

where Y is an lower triangular matrix with vanishing diagonal elements satisfying [Y, B] = 
A. 

(iii) Of course all results remain true if we transpose all matrices. 

D.3. The two-matrix model 

Some problems unsolved in the case of the general g-matrix model can be solved in 
the two-matrix case. For example, let i? be a local matrix such that 

Bmn 7^0 form — r-|-l<?T.<m-t-l(r>2). 

We look for B such that there exists another local lower triangular matrix A such that 

Amn 7^ only for m > n , 

and which satisfies 

[B,A] = 1. 
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Note that if B and A are solutions to this problem, 5" ^BS and S ^AS are also solutions 
if 5" is a diagonal matrix, Smn = SmSmn- Therefore we can fix some of the matrix elements 

Bmn, for example Bm,m+1 or Bm,m-r+l- 

We shall construct a set of solutions to this problem by adding one assumption and 
show that the solutions then correspond to the general two-matrix model: we shall assume 
that we can find another local matrix B such that X = A + B commutes with B, 

X = A + B, [X,B] = ^ [B,B] = 1 , 

with 

Amn = for m < n, -Bmn 7^ only form — l<n<m + s — 1 . 
The decomposition is ambiguous for the elements n = m — 1, so we leave A^^rn-i arbitrary. 



The matrix X is also local. Then according to the result of app. \D.3i we have 

A + B = V'iB) , 

where V is an arbitrary polynomial of degree s. It follows that Bmn = for n < m + s, 
and therefore 

[B, X] = ^ [B, A] = [B,B] = J . (D.3) 

The commutators [B, A] and [B, B] have non-vanishing matrix elements Jmn only for 
m > n, n > m — r, respectively. Therefore the matrix J has the same width as B. The 
equation J = 1 then yields just the right number of equations, i.e. r recursion relations, 
which determine the matrix B. For m = n, we use the commutator [B, A] to find 

-^n,n+l-^n+l,n ^n,n— l,n 1 ; 

and thus 

If we write the other equations in terms of A, we get a subset of the loop equations 
( ^.17| ). If instead we use the commutator [B, B] = 1, we can use the result (ii) of app. \L).^ 
the role of B played by B. It follows that we can introduce a polynomial V{B) of degree 
r such that 

Bmn=[V'iB)]^^ forn<m, 
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which can be rewritten 

B + A = V'{B), Amn = for n<m. 
We recognize exactly the structure of the two matrix model. In particular 

[B^,A^] = 1, 

and 

Note that we could have stated the initial problem more symmetrically by starting directly 
from the commutator [B, B] = 1. 

Remarks. 

(i) In the two-matrix case we have proven with one additional assumption eqs. 
and thus we know the form of the generators of the KdV flows of the two-matrix model. 

(ii) The ^2 invariant two-matrix problem is even simpler because then the matrices 
satisfy B2 = Bf, and the model can described in terms of only two matrices B and A 
satisfying 

A + B^ = V'{B), (DA) 

and thus 

E + = ^ {V{B) + V{B^) - B^B) . (D.5) 

D.4. The general ^2 invariant model 

In the case of the multi-matrix model with ^2 invariance, the general problem of 
finding solutions to the canonical commutation relations can be reduced to that of finding 
the solution of commutation relations C] = 1 where B and C are two local matrices, 
respectively symmetric and antisymmetric, satisfying 

Bmn = for |m — n| > g , Cmn = for |m — ?t,| > r . 

Let us introduce a matrix D, in general non-local, such that 

B = D^-^ and Dmn = forn > m + 1 . 

Note that although D is not local, the calculation of matrix elements of powers of D 
involves only finite sums. 



160 



It is easy to verify that the matrix elements of D can be systematically calculated 
from the matrix elements of B. The first equation is 

m+q—2 

Bm,m+q-l = Y\. ^h3 + '^ ■ 

j=m 

We assume that the elements Bm,m+q-i never vanish. This equation determines all matrix 
elements Dm,m+i except for the first q — 2 which remain arbitrary. The next equation is 

m+g — 3 m+g — 2 

Bm,m+q-2 — JJ^ + l . 

j=m i=m 

This equation determines Dmm again when the q — 2 first elements are given. 
More generally for k >1, 

m+k — 2 m+q—2 i—1 

,m+k—l 

= n ^^'^+1 Yl A,i+fc-g+i n ^M+i + known terms . 

j=m i=m l=i-\-k — q+l 

We always find q — 1 terms recursion equations; however for the first values of m the 
situation is different and only the first k — 1 elements are undetermined. For /c < 1, the 
equations are similar but all elements are determined. 

Let us now consider the equation C] = 1. We first specialize to the matrix element 
n = m + q-\-r — 2 because only two terms contribute, 

m+r+g— 3 m+q—2 

Cm,m+r-l JJ^ ~ W Djj^iCm+q-l,m+q+r-2 = . 

j=m+r—l 3=^1 

Setting 
we obtain 

c' — c' = n 

^m,m+i — 1 ^m+q—l,m+q+r—2 ^ • 

This equation determines all C'^^rn+r-i except the first q — 1- However we remember that 
the q — 2 elements of D are not determined. We can use them to reduce C^^ to a 

constant. Also if D has a continuum limit, then we expect that only the constant solution 
corresponds to a matrix C which also has a continuum limit. 

Generalizing this argument ton — m>g — 1, we conclude that C can be written 

C = W{D) - L , 
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where is a polynomial of degree r — 1 and L is a lower triangular matrix with vanishing 
diagonal elements, Lmn = for n > m. Since C is antisymmetric, we have also 



= -C => L + = W + , 
which determines C and L in terms of W. Moreover we have 

[B,C] = [L,B] = 1 . 

Now we consider the equations corresponding torn < n < m + q — 1. The other 
equations are satisfied by symmetry. A first equation is 

L'm,m — l ■Bm—l,m+q—2 -^mjTn+ij— 1 -^m+q — l,m+q— 2 — • 

The difference between this equation and the previous ones is that the first term vanishes 
for m = 0. Therefore it seems likely that the only solution compatible with a continuum 
limit is Lm,m-i = 0. The argument extends to all non-diagonal matrix elements of the 
commutator and thus 

Lm,m-k = for A; < - 1 . 
Finally the last equation is 

Lm,m—q+l ^m—q+l,m -^m,m+q— 1 ^m—q+ljm — 1 ) 

which leads with the same arguments to 

r R 

^m..m. — (7-1-1 ^m. — c 



q-l 

In terms of W , these equations read 

Wm,u + Wn,m = for |n - m| < g - 1 , 

2m (D.6) 

,m— q+1 — 9+1,771 — / ^ X 

{q — i)nrn-q+l,m 

Note that the number of equations we have obtained is equal to the number of coefficients 
of S. 
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An action principle. Eqs. ( p.6|) , which determine B, can be derived from an action 
principle with action 

S{B) =tvV{B)-J2 In Bm-q+l,m , 

r-1 (D.7) 

V{B) = ^14trSi+^/(«-^) . 

fc=0 

KdV flows. KdV flows are generated by antisymmetric matrices E which can be 
reconstructed from analogous considerations. One finds 

d 

^mn^ -Q^[V{B)]^^ for7i>m, 
where V{B) is defined in eq. ( p.7| ). 
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